$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
1. Experimentele workflow
Figuur 1 illustreert een georganiseerde en modulaire workflowpijplijn die is ontworpen om moleculaire biomarkers te categoriseren in categorieën Risico- of Diagnostische indicaties met behulp van machine learning- en deep learning-technieken. Hier volgt een gedetailleerde stapsgewijze uitleg van het proces:
- Dataset-acquisitie
De dataset is verzameld uit MarkerDB 2.0, een zorgvuldig samengestelde database van moleculaire biomarkers.
- Datavoorverwerking
Dit zorgt voor datakwaliteit en consistentie door nulwaarden en ontbrekende vermeldingen te verwerken, categorische features te coderen met LabelEncoder en numerieke features (bijv. entrez_gene_id) te schalen met StandardScaler.
- Feature-engineering
Het selecteren van relevante kenmerken die de biomarkerclassificatie beïnvloeden, datalekken voorkomen en de modelprestaties verbeteren, wordt aanbevolen.
- Modelontwikkeling
De gebruikte benaderingen zijn Logistic Regression, Random Forest, XGBoost, terwijl in DL: MLP, MLP met LR Scheduler, AutoInt, TabNet. De uitgevoerde activiteiten zijn hyperparameterafstemming, kruisvalidatie en optimalisatie voor classificatietaken.
- Datasplitsing
Hier wordt de trainingsset (80%) gebruikt om modelparameters te leren, en de testset (20%) wordt gebruikt om het model te evalueren op onzichtbare data. In de laatste fase worden evaluatie-metrics gebruikt om ervoor te zorgen dat prestaties holistisch worden gemeten, vooral in onevenwichtige klassenscenario's.

Figuur 1: Stroomdiagram van het voorgestelde werk. De workflow van het voorgestelde probleem wordt in deze figuur weergegeven. De stappen zijn data-preprocessing, feature engineering, modelontwikkeling, splitsing, modelevaluatie met behulp van metrics en biomarker-analyse. Klik hier om een grotere versie van deze figuur te bekijken.
2. Datasetbeschrijving
De naam van de dataset is all_sequence_variants met de csv-extensie. Er zijn kolommen met namen als 'id', 'variation', 'position', 'external_link', 'gene_symbol', 'entrez_gene_id', 'reference', conditions', 'indication_types', waarbij allemaal objecttype zijn behalve id als int-type en entrez_gene_id als float-type. Het totale aantal rijen is 42818 met 9 kolommen. Bovendien bevat de dataset nulwaarden voor variatie, positie, external_link, entrez_gene_id, referentie, condities en indication_types. De dataset is afkomstig van de link gegeven als: https://markerdb.ca/downloads17
3. Dataset-preprocessing
In deze fase begint de voorverwerking van data met het verzamelen van dataset-informatie, beschrijvingen, het identificeren van duplicaten en het lokaliseren van ontbrekende/nul waarden. Er is geen impact van het external_link kenmerk wanneer het wordt gemeten met de correlatiecoëfficiënt, en daarom wordt het verwijderd. De nulwaarden van categorische en numerieke kolommen worden respectievelijk gevuld met de mediaan en het gemiddelde. Ten slotte is de vorm van input en output in de dataset (42787, 6)(42787,). De invoerkolommen zijn 'id', 'variatie', 'positie', 'gene_symbol', 'entrez_gene_id', 'voorwaarden', terwijl de doelkolom indication_types is.
4. Biomarkers: Categorieën, data-uitdagingen en de rol van moleculair risico
Een biomarker is een kwantificeerbare eigenschap die ofwel normale of abnormale biologische activiteiten of reacties op verschillende blootstellingen of interventies aangeeft. Biomarkers kunnen moleculair, histologisch, radiografisch of fysiologisch van aard zijn en zijn onderverdeeld in zeven primaire categorieën: (1) Een diagnostische biomarker bepaalt of een individu een specifieke ziekte of aandoening heeft en of deze binnen een subtype valt, (2) Een monitoringbiomarker geeft de progressie van een ziekte en de respons op behandeling aan, (3) Een responsbiomarker geeft aan of een patiënt reageert op een medicijn of therapie, zoals een verandering in hartslag, (4) Een voorspellende biomarker kan vaststellen of iemand risico loopt op het ontwikkelen van een bepaalde ziekte en hoe hij of zij kan reageren op een bepaalde behandeling, (5) Een prognostische biomarker kan inzicht geven in de kans op ziekteprogressie of terugkeer als een patiënt vatbaar is voor een specifieke gezondheidsuitkomst, (6) Een risicobiomarker beoordeelt het potentiële risiconiveau voor een aandoening voordat een patiënt een officiële diagnose krijgt, en (7) Een veiligheidsbiomarker beoordeelt het potentieel voor toxische of bijwerking die kunnen ontstaan door behandeling. Deze studie richt zich voornamelijk op de analyse van moleculaire biomarkers die verband houden met risicobeoordeling en diagnostiek, zoals weergegeven in Figuur 2. Gezien de enorme hoeveelheid klinische data die wereldwijd beschikbaar is, zijn biomarkers essentieel voor het sturen van klinische beslissingen, het bevorderen van onderzoek en het faciliteren van gepersonaliseerde geneeskunde.

Figuur 2: Biomarkercategorie voor klinische besluitvorming. De beslissingen van de biomarker worden in dit diagram weergegeven. Op basis van de analyses die op de biomarkers worden toegepast, worden klinische beslissingen genomen en worden ML-algoritmen geïmplementeerd. Klik hier om een grotere versie van deze figuur te bekijken.
Het vinden en begrijpen van klinisch bruikbare biomarkerinformatie is een uitdaging vanwege het grote scala aan typen, kenmerken en kwaliteit van biomarkers. De snelle groei van biomarkerstudies in de afgelopen twee decennia heeft de toegang tot uitgebreide en actuele gegevens, vooral voor moleculaire biomarkers, verder bemoeilijkt. Dit heeft geleid tot problemen zoals inconsistente bevindingen, het herontdekken van bestaande biomarkers, tegenstrijdige resultaten en over het hoofd geziene kansen omgevestigde biomarkers te benutten. De opkomst van "omics"-metingen heeft dit probleem verergerd, wat heeft geleid tot een proliferatie van nieuwe moleculaire biomarkers in klinische tests en de literatuur, waardoor het steeds moeilijker wordt om biomarkerkennis effectief te navigeren en toe te passen19.
MarkerDB 2.0 biedt toegang tot diverse moleculaire biomarkers voor experimentele doeleinden20. Verschillende biomarkers worden gebruikt in klinisch onderzoek en gezondheidszorgtoepassingen om ziekten te diagnosticeren, voorspellen en monitoren. Onder deze eiwitbiomarkers gebruiken biofluiden zoals serum, bloed, cerebrospinaal vocht (CSF), peritoneale vloeistof, vruchtwater, plasma en urine om aandoeningen zoals diabetes mellitus, tuberculose, het syndroom van Down en myopathie-gerelateerde aandoeningen te analyseren. Evenzo gebruiken genetische biomarkers gensymbolen zoals LRP1, LPP, MAPT en SPI1 om aandoeningen te analyseren zoals leeftijdsgebonden maculadegeneratie, allergische ziekten, de ziekte van Alzheimer, kanker en astma. Bovendien helpen multimodale moleculaire biomarkers het patiëntrisico en diagnostische indicatoren te meten21. Dit werk benadrukt de noodzaak van gezondheidszorganalyse met ML-benaderingen om ziektediagnose, voorspelling en gepersonaliseerde behandeling te verbeteren.
5. Modellering van statistische ML-benaderingen voor biomarkerontdekking
De toepassing van ML-benaderingen heeft de gezondheidszorg op tal van manieren getransformeerd. Specifiek zijn moleculaire biomarkers geanalyseerd met behulp van diverse ML-technieken, waaronder Principal Component Analysis (PCA), K-means clustering (KMA), Nearest Neighbor Analysis (NNA) en neurale netwerkalgoritmen22. Deze modellen identificeren ingewikkelde patronen uit gekozen kenmerken, beheren onvolledige of inconsistente data en ondersteunen realtime tracking van ziekteverloop. Naast het diagnosticeren van ziekten bieden ML-algoritmen belangrijke inzichten in patiëntkwesties door data-analyse en visualisatie, waardoor snelle en gerichte zorg mogelijk is23. Hierdoor verbetert dit de gezondheidsuitkomsten van patiënten in bijzondere omstandigheden. Bovendien wordt het vermogen van ML om biomarkerontdekking en therapeutische resultaten over verschillende ziekten, waaronder de analyse van moleculaire biomarkers, te revolutioneren.
ML-algoritmen worden ingedeeld in vijf hoofdtypen, zoals weergegeven in Figuur 3. Supervised learning is het proces van trainen op gelabelde datasets om input-outputrelaties te begrijpen, waardoor het ideaal is voor taken zoals classificatie en regressie, zoals beslissingsbomen en ondersteunende vectormachines. Ongecontroleerd leren identificeert patronen binnen niet-gelabelde data en wordt vaak toegepast voor clustering en dimensionaliteitsreductie, met technieken zoals k-means clustering en principal component analysis. Semi-supervised learning voegt zowel gelabelde als niet-gelabelde data samen om voorspellingen te genereren. Reinforcement learning richt zich op besluitvorming gebaseerd op feedback uit de omgeving, en wordt vaak gebruikt in gaming en robotica, met methoden zoals Q-learning en deep reinforcement learning-netwerken. Deep learning maakt gebruik van kunstmatige neurale netwerken met meerdere lagen om ingewikkelde patronen in data te identificeren. Deze techniek helpt bij het voorspellen van biomarkers voor veelvoorkomende ziekten, waaronder kanker, beroerte, de ziekte van Alzheimer en de ziekte van Parkinson. Talrijke klinisch goedgekeurde toepassingen maken gebruik van deze technologie24,25. De precieze diagnose van patiënten met kwaadaardige tumoren hangt doorgaans af van de verworven en pathologische analyse van weefselmonsters, die cruciale informatie opleveren over histologische graad, subtype, stadium en diverse andere tumorbiomarkers.

Figuur 3: Categorisering van machine learning-algoritmen voor klinische besluitvorming. De typen ML-algoritmen in het diagram worden geïllustreerd om het toegepaste methodologium beter te begrijpen. Klik hier om een grotere versie van deze figuur te bekijken.
Statistiek en ML zijn twee aparte disciplines die vaak samenkomen. Statistiek omvat het verzamelen, analyseren en interpreteren van data, meestal met kleinere, goed gedefinieerde datasets, en richt zich voornamelijk op het vergelijken en samenvatten van informatie. ML daarentegen is een subset van kunstmatige intelligentie die voorspellende modellen creëert, vaak gericht op grote en complexe datasets. Hoewel statistiek een belangrijke rol speelt bij bevestigend onderzoek en het begrijpen van onderliggende mechanismen, is ML effectiever in verkennend onderzoek, het ontdekken van complexe patronen en het beheren van ontbrekende data. Methoden zoals logistische regressie kunnen worden beschouwd als zowel statistische modellen als supervised ML-modellen, wat de vage verschillen tussen de twee domeinen illustreert. Uiteindelijk moeten onderzoekers hun onderzoeksdoelstellingen en de aard van hun data evalueren om de meest geschikte methode te kiezen.
6. Modellering van MLP en varianten van MLP
MLP, ook bekend als een multilayer perceptron, is een type neurale netwerkarchitectuur met meerdere neuronlagen tussen de invoer- en uitvoerlagen, zoals geïllustreerd in Figuur 4. MLP werd geïntroduceerd in de jaren 1950 en kreeg brede toepassing door de vooruitgang en vooruitgang in backpropagation in de jaren 1980, wat hielp het verlies te minimaliseren. De basisprincipes van MLP zijn neurale netwerken, gewichten en biaswaarden voor het berekenen van de output. Ten slotte om de activatiefunctie en drempelwaarde te kennen om de output in de outputlaag te krijgen. In de voorgestelde methodologie komen de knooppunten in de inputlaag, verborgen laag en outputlaag overeen met klinische en genkenmerken, volledig verbonden lagen met ReLU-activatie, en een enkele neuron met sigmoidactivatie om respectievelijk de binaire uitkomst van indicatietype (Risico, Diagnostisch) te voorspellen. De belangrijkste voordelen zijn dat ze eenvoudig te trainen en te interpreteren zijn. De eenvoudige MLP kan echter geen grote datasets verwerken en is gevoelig voor de leersnelheidvan 26. Daarom wordt MLP met Learning Rate Scheduler overwogen om deze nadelen te overwinnen met de architectuur die in Figuur 3B wordt gegeven. MLP met LR is een krachtig dynamisch mechanisme met steile verval, exponentiële verval, inverse tijdsverval, ReduceLROnPlateau en cosinus-annealing. In plaats van een vaste LR-waarde verandert de leersnelheid op basis van de modelprestaties en training.

Figuur 4: Modellering van meerlagige Perceptron en MLP met een LR-scheduler. (A) Multi-Layer Perceptron: De structuur van MLP wordt geïllustreerd met Input Layer, Hidden Layer en Output Layer. De eerste laag ontvangt de input en verwerkt deze in de tweede laag met activaties, en de output wordt ontvangen in de laatste laag. MLP is eenvoudig te implementeren. (B) MLP met LR Scheduler: Dit lijkt op MLP; er wordt echter een Learning Rate Scheduler toegevoegd om de trainingssnelheid te regelen voor een betere convergentiesnelheid. Het leertempo is op een plateau afgenomen. Dit vermindert de minima van overschrijdingen. Klik hier om een grotere versie van deze figuur te bekijken.
Deze aanpak is efficiënter, betrouwbaarder en effectiever voor tabelgegevens. Daarnaast zijn er voordelen zoals stabiele convergentieverbetering op ruisachtige data, verbeterde generalisatie en minder inspanning bij hyperparameterafstemming. Voor sommige voorspellingstaken kunnen er echter moeilijkheden ontstaan, zoals vast komen te zitten in lokale minima. Bovendien worden in de brede en diepe MLP-architectuur zoals gepresenteerd in Figuur 5 brede en diepe componenten geïntroduceerd in de MLP om uit te blinken in correlatiemapping, het onthouden van regels en het leren van nieuwe patronen, terwijl hoge standaarden van generalisatie behouden blijven. Door deze aspecten te combineren tot de uitvoerlaag, wordt de binaire uitvoervoorspeld met 27. Er zijn echter beperkingen aan het generaliseren van onzichtbare kenmerken en overgeneralisatie. Samenvattend wordt MLP met batchnorm en uitval gebruikt om het trainingsproces te normaliseren, waarbij uitval wordt geïntegreerd voor het leren en wordt generaliseerd van onbekende attributen. Deze techniek kan met name worden toegepast op hoogdimensionale data en helpt overfitting te voorkomen. Desalniettemin heeft deze techniek last van beperkingen in batchgrootte, een verhoogd geheugenverbruik (door Batch Normalization-lagen) en is het geen universele oplossing28,29.

Figuur 5: Modellering van MLP-breed en diep netwerk. Het diagram illustreert het diepe pad van de MLP, dat niet-lineaire relaties kan vastleggen en de output combineert met een sigmoid-eenheid voor classificatie. Deze architectuur legt patroonleer en kennis vast voor heterogene biomarkergegevens. Klik hier om een grotere versie van deze figuur te bekijken.
7. TabNet
De TabNet is een deep learning-model dat bekend staat als Tabular Network, ontwikkeld door het Google Research-team30. De belangrijkste motivatie van het TabNet-model is het overbruggen van verschillen tussen de DL-benaderingen voor beeld-, tekst-, spraak- (CNN, Autoencoder, Transformers) en boomgebaseerde ML-benaderingen (XGBoost, Random Forest, LightGBM). Bovendien hebben de beperkingen van DL-benaderingen een enorme impact op het TabNet-model, waarbij MLP een geparametriseerd model is. Bovenal ontbreekt het aan interpretatie van DL-benaderingen bij echte problemen, zelfbegeleid leren en sequentiële leercapaciteit31. De TabNet-modellen zijn beperkt tot toepassingen zoals verzekeringsrisicovoorspelling32, schatting van steenkoolasgehalte33, voorspelling van chemische toxiciteit34 en detectie van fraude in educatieve tests35. De architectuur van TabNet is weergegeven in Figuur 6.

Figuur 6: Modellering van een tabulair netwerk. Het Tabulaire netwerk neemt invoerfuncties in een tabelformaat en past de transformatortechniek op een sequentiele manier toe in elk blok. Het GU-blok, ook wel het Gated Linear Unit-blok genoemd, regelt de informatiestroom door het netwerk, wat de selectie van functies en stabiel leren vergemakkelijkt. De output wordt in de outputlaag gepresenteerd uit de geaccumuleerde output. Klik hier om een grotere versie van deze figuur te bekijken.
Uit het bovenstaande architectuurdiagram is het systeem verdeeld in vier componenten: Inputblok, Shared Feature Transformer, Sequentiële Beslissingsstappen en de verliesberekening met behulp van optimalisatie. In de categorische kenmerken worden de kenmerken omgezet naar gehele getallen na de inbeddinglaag. Tegelijkertijd worden de continue kenmerken genomen zoals ze zijn. Ten slotte worden de features omgezet in een uniforme featurevector waarbij x ∈ Rd. In de volgende stap worden volledig verbonden lagen met activatiefuncties geïmplementeerd, gevolgd door een transformatie om de data om te zetten in een vector, en wordt een beslissing genomen. In de derde stap worden de aandachtstransformator, maskering, beslissingstransformator en voorspellingsaccumulatie achtereenvolgens gebruikt. De doelstellingen van deze stap zijn selectieve aandacht, sequentiële redenering en beslissingspaden. In de laatste stap worden binaire kruisentropie of vergelijkbare methoden gebruikt om het verlies te vinden en de waarde te optimaliseren. Verder wordt regularisatie uitgevoerd met spaarzame aandacht om overfitting te voorkomen.
8. AutoInt
Automatisch Feature Interaction Learning via Self-Attendentive Neural Networks (AutoInt) is de volledige vorm van AutoInt, geïntroduceerd in Click-through rate (CTR) voorspelling36 en later toegepast op diverse toepassingen, waaronder softwaredefectvoorspelling37, aanbevelingssystemen38 en genoomdetectie39. Er zijn beperkingen in de bestaande ML-benaderingen, waaronder het omgaan met hoogdimensionale, spaarzame data, combinatorische kenmerken van hogere orde, het begrijpen van de voorspellingsmogelijkheden van de eigenschap en de noodzaak van handmatige feature-engineering. Al deze genoemde uitdagingen worden door het AutoInt-model aangepakt met efficiëntie, effectiviteit en uitlegbaarheid. Het doel van AutoInt in de voorgestelde methodologie is te voorspellen of de gegeven variant van een risico- of diagnostisch indicatietype is. Er zijn zeven stappen voorgesteld voor de gegeven dataset, zoals weergegeven in Figuur 6, samen met de broncode. Ten eerste worden het lezen van data en het voorverwerken geïmplementeerd om alle categorische kenmerken als gehele getallen te coderen en de labels te mappen naar 0 en 1 voor respectievelijk Risico en Diagnostisch. In de tweede stap wordt feature-representatie uitgevoerd om de embedding-matrix te creëren; op zijn beurt wordt de matrix omgezet om de numerieke kenmerken te normaliseren. Er zijn interactieve lagen waarbij de neutrale automatische interactie tussen de kenmerken wordt geleerd zonder handmatige interactie met de transformatoren. Bovendien leren meerdere hoofden verschillende soorten interacties, en uiteindelijk worden de kenmerken samengevoegd, wat de representatie versterkt. Residuele verbindingen zijn aanwezig om zowel lagere als hogere orde interacties te leren, zoals 1 kenmerk, 2 kenmerken, 3 kenmerken, enzovoort. Ten slotte worden vanuit het residuele netwerk de uitvoer doorgegeven aan de sigmoidfunctie voor binaire waardeberekeningen, zoals weergegeven in Figuur 7.

Figuur 7: Modellering van het AutoInt-netwerk. Het AutoInt-netwerk leert automatisch features, waarbij de embedding-laag de features neemt en ze in de volgende laag omzet met behulp van een zelf-aandacht mechanisme en residuele verbindingen. De laatste component bevat de MLP-laag en activatiefuncties om de output te produceren. Klik hier om een grotere versie van deze figuur te bekijken.