3  Material und Methoden

Die deutsche Multizentrische Allergiestudie (MAS-90) wurde initiiert, um den Voraussagewert verschiedener klinischer und immunologischer Parameter sowie den Einfluss von Umweltfaktoren auf die Entwicklung von atopischen Krankheiten in früher Kindheit zu untersuchen.

Fünf deutsche Städte (Berlin, Düsseldorf, Freiburg, Mainz und München) waren das Einzugsgebiet für die Rekrutierung von Neugeborenen für diese prospektive Studie. 1314 von 7609 Neugeborenen, die im Zeitraum zwischen dem 1.1.1990 und dem 31.12.1990 geboren wurden, wurden in die Geburtskohorte aufgenommen. 499 Kinder (38%) wurden als Patienten mit einem hohen Risiko für Atopie eingestuft. Sie hatten mindestens zweiatopische Familienmitglieder und/oder wiesen eine IgE-Konzentration im Nabelschnurblut größer als 0,9 kU/l auf. Die Kinder der Kontrollgruppe (815 Kinder = 62%) wurden unter den verbliebenen Neugeborenen zufällig ausgewählt. Sie hatten entweder einen oder keinen Atopiker in der nahen Verwandtschaft und einen Nabelschnurblut-IgE-Wert kleiner als 0,9 kU/l. Eltern wurden als Atopiker definiert, wenn sie über eines der folgenden Symptome beziehungsweise eine der folgenden Krankheiten berichteten: Atopisches Ekzem, allergische Rhinitis, Asthma, Nahrungsmittelunverträglichkeit, Pollen-, Hausstaubmilben-, Katzenhaar-, Hundehaar- oder Schimmelpilzallergie.

	Abb. 10: Studiendesign und Auswahl der Neugeborenen in der Multizentrischen Allergie Studie.

Die Kinder der Kohorte kamen regelmäßig mit ihren Eltern zu Nachuntersuchungen jeweils im Lebensalter von einem Monat, drei Monaten, sechs Monaten, 12 Monaten, 18 Monaten, 24 Monaten und anschließend einmal jährlich. Die Untersuchungen beinhalteten:

Fragebogen und Interview: Die Eltern füllten einen Fragenkatalog aus und antworteten auf gezielte Fragen zu Atopie-spezifischen Symptomen und Erkrankungen seit der letzten Untersuchung. Außerdem wurden Fragen zur Ernährung, Entwicklung des Kindes, zu Umweltfaktoren (Zigarettenrauch, Haustiere), zur häuslichen Situation und zu psychologischen Problemen gestellt. Die MAS-Fragebögen und Interviews wurden auf der Grundlage der ISAAC (International Study on Asthma and Atopy in Childhood)-Fragebögen entwickelt [71].

Körperliche Untersuchung: Bei jedem Termin wurde eine standardisierte körperliche Untersuchung durchgeführt. Dabei war ein Schwerpunkt die standardisierte Dokumentation von Hautläsionen.

Blutentnahme: Im Serum der Kinder wurden im Alter von 1, 2, 3, 5, 7 und 10 Jahren Gesamt-IgE- und spezifische IgE-Antikörperkonzentrationen gegen neun gängige Nahrungsmittel- und Inhalationallergene bestimmt. Die spezifischen IgE-Antikörper gegen Hühnerei, Kuhmilch, Sojabohnen, Weizen, Hausstaubmilben, Katzenhaare, Gräsergemisch und Birkenpollen wurden mit einem Radioallergosorbent Fluorenscence Immunoassay (CAP-RAST-FEIA, Pharmacia and Upjohn, Freiburg) bestimmt. Die IgE-Werte wurden für die statistischen Analysen log-transformiert, um eine Gauss’sche Verteilung zu erreichen.

Aus EDTA-Blut wurde die genomische DNA extrahiert. Insgesamt steht die genomische DNA von 888 Kindern sowie 1231 Eltern zur Verfügung. Insgesamt ergaben sich 382 komplette Trios (Mutter, Vater, Kind).

Lungenfunktionsuntersuchung: Lungenfunktionstests bei den Kindern wurden an allen deutschen Zentren mit dem gleichen Ganzkörperplethysmographen (Master-Lab, E Jäger, Würzburg) durchgeführt. Bestimmt wurden das forcierte expiratorische Volumen der ersten Sekunde (FEV₁) und die forcierte Vitalkapazität (FVC).

Ein pharmakologischer Provokationstest wurde mit Histamin-dihydrochlorid in aufsteigenden Dosen durchgeführt. Die Ausgangsdosis betrug 0,125 g/l Histamindihydrochlorid, anschließend wurden immer höhere Dosen (0,5 g/l; 2,0 g/l; 8,0 g/l) verabreicht, bis der Patient dies nicht mehr tolerierte oder die FEV₁ um mehr als 20% absank. Die Histaminkonzentration, die einen Abfall der FEV₁ von 20% bewirkte (PC₂₀FEV₁), wurde durch lineare Interpolation der beiden letzten Messwerte berechnet. Wurde schon nach der Anfangskonzentration ein Abfall von mindestens 20% verzeichnet, wie es bei 5% der Studienteilnehmer der Fall war, so wurde der Nullwert als erster Datenpunkt verwendet. Wenn die Maximaldosis von 8,0 g/l nicht ausreichte, um einen 20% Abfall zu bewirken, wurde diese Dosis für die Berechnung verwendet.

Staubanalysen: Teppichstaubproben der häuslichen Umgebung wurden im Alter von 6 und 18 Monaten, 3 Jahren, 4 Jahren, 5 Jahren auf die wichtigsten Innenraumallergene überprüft: Hausstaubmilben, Katzen- und Hundehaar.

	Abb. 11: Untersuchungen von Geburt bis zum siebten Lebensjahr bei den Kindern der MAS-Kohorte.

Bei der Untersuchung am siebten Geburtstag (± sechs Wochen), nahmen noch 939 (71,5%) der Kinder teil. Die Phänotypen der allergischen Erkrankungen wurden wie folgt definiert:

Als Asthmatikergalt ein Kind bei mindestens zwei aufgetretenen Episoden von pfeifender Atmung mit Luftnot.

Das atopisches Ekzem wurde bei einer Kombination von mindestens drei typischen Hautveränderungen an drei typischen Körperstellen diagnostiziert.

Allergische Rhinitis im Alter von 7 Jahren wurde definiert als Auftreten von [Seite 36↓]Niesattacken und/oder einer laufenden, verstopften, juckenden Nase ohne Erkältung während der letzten 12 Monate. Zusätzlich mussten spezifische IgE-Antikörper (≥0,7 kU/l) gegen Gras- und/oder Birken-Pollen vorliegen.

Ein Kind galt als sensibilisiert, wenn der spezifische IgE-Wert bei mindestens einem von neun Hauptallergenen (Milbe, Hund, Katze, Kuhmilch, Ei, Soja, Weizen, Birke, Gras) ≥0,7 kU/l betrug.

Mütterliches Rauchen wurde als Tabakrauchexposition des Kindes gewertet, da Mutter und Kind sich normalerweise in enger räumlicher Nähe zueinander befinden. Pränatale Exposition wurde als mütterliches Rauchen während der Schwangerschaft definiert und von postnataler Exposition ausgegangen, wenn die Mutter bei ≥ 2 der ersten drei Nachfolgetermine angab zu rauchen oder zu ≥ 4 Zeitpunkten der insgesamt 7 Untersuchungen. Uneindeutige Angaben wurden als fehlend definiert.

K. Barnes rekrutierte zwischen 1993 bis 2002 783 Individuen auf Barbados (Karibische Inseln). Die Population setzt sich aus Asthmatikern und deren unmittelbaren Verwandten zusammen. 13-20% der gesamten Bevölkerung von Barbados leiden unter Asthma. Der Anteil an Atopikern wird noch höher geschätzt. Barnes und andere haben bei der Bevölkerung von Barbados eine extrem hohe mittlere Gesamt-IgE-Konzentration gefunden, die die hohe Allergenexposititon in dieser tropischen Umgebung widerspiegelt. Das Hauptallergen der Hausstaubmilbe Dermatophagoides pteronyssinus beispielsweise liegt mit 60 µg/g Hausstaub [72] rund 30 mal höher als die „Risikokonzentration“ in den USA bzw. Europa [73]. Die exzessive Exposition erhöht die Wahrscheinlichkeit, dass genetisch prädisponierte Individuen atopische Erkrankungen ausbilden. Es gibt keine parasitären Erkrankungen auf Barbados, die die hohen IgE-Werte erklären könnten. 7% der Barbados Population rauchen, unter den Asthmatikern sind nur 5% Raucher. Alle Kernfamilien wurden über einen Asthma-Probanden ausgewählt. Die Asthmatiker wurden systematisch aus der Notaufnahme des Queen Elizabeth [Seite 37↓]Hospitals oder aus Arztpraxen auf der Basis von folgenden Einschlusskriterien ausgewählt: Erstens sollten sie älter als sechs Jahre alt sein, zweitens sollte eine ärztliche Diagnose Asthma gestellt worden sein und drittens sollten beide Elternteile verfügbar sein oder bei einem fehlenden Elternteil verfügbare Geschwister, um den fehlenden Genotyp abzuleiten. Es wurden alle Kernfamilienmitglieder eingeschlossen und wenn möglich, der Stammbaum so weit wie möglich ausgedehnt.

250 Kinder mit Verdacht auf Asthma im Alter von 5 bis 18 Jahren (Mittelwert=10,8 Jahre) aus dem Einzugsgebiet von Freiburg wurden in der Kinderklinik des Universitätsklinikums Freiburg untersucht. Es wurde eine standardisierte Untersuchung durchgeführt, die Fragebögen zu klinischen Symptomen und Medikation, die Bestimmung von totalem und spezifischem IgE, Pricktests, Lungenfunktionstests, bronchiale Histaminprovokationen sowie die Bestimmung des belastungsabhängigen Abfalls des FEV₁-Wertes mittels Laufbandbelastung beinhaltete. Die Teilnehmer wurden angewiesen, frühzeitig die Asthma- bzw. Allergiemedikation vor der Untersuchung abzusetzen. Die Probanden wurden mit Pricktests durch intrakutane Gaben von 17 häufigen Allergenen (Hausstaubmilben, verschiedene Gräser- und Baumpollen, Aspergillus fumigatus, Alternaria alternata, Cladosporium herbarum, Hund, Katze, Hase, Ente und Pferd) zusammen mit Positiv- (Histamin) und Negativkontrollen getestet. Der Durchmesser der resultierenden Schwellung wurde nach 15 Minuten bestimmt. Der Pricktest wurde als positiv gewertet, wenn die Schwellung mindestens halb so groß wie die der Positivkontrolle war.

Spezifische IgE-Antikörperkonzentrationen wurden mit einem enzymgekoppelten Immunosorbent Assay gegen zwei Mischungen von Gräserpollen, Milbenallergenen (Dermatophagoides pteronyssinus und Dermatophagoides farinae), Katze, Hund, Haselnuss, und Birkenpollen bestimmt (Magic Lite, Chiron Diagnostics, Fernwald). Ein Test galt bei ≥1,43ML [Seite 38↓]Units als positiv. Die Gesamt-IgE-Konzentration wurde mit Hilfe eines Enzym Allergosorbent Test (Phadezym, Pharmacia, Uppsala, Schweden) gemessen.

Die Lungenfunktion und Histaminprovokation wurden nach einem dem der MAS-Studie vergleichbaren Standardprotokoll durchgeführt. Um belastungsinduzierte FEV₁-Abfälle zu evaluieren, wurden die Probanden 6 Minuten lang standardisiert körperlich belastet (Laufband). Die erste spirometrische Messung und Peak flow-Bestimmung wurden 2-3 Minuten und die zweite Messung 5-6 Minuten nach der Belastung durchgeführt.

Die Diagnose eines Asthma bronchiale stützte sich auf asthmatische Symptome, den Gebrauch von antiasthmatischen Medikamenten sowie einer nachgewiesenen BHR (Histaminprovokation und/oder Laufbandbelastung).

[Seite 39↓]

Die durch Venenpunktion gewonnenen Blutproben wurden in EDTA-Röhrchen abgenommen und bei –20°C gelagert. Zur Extraktion der genomischen DNA aus Vollblut wurde der QIAamp DNA Midi Kit von QIAGEN nach dem Arbeitsprotokoll des Herstellers verwendet: 2 ml EDTA-Blut wurden in ein Reagenzröhrchen überführt, mit 200 μl QIAGEN Proteaselösung und 3 ml Puffer AL gründlich gemischt und für 20 Minuten bei einer Temperatur von 70°C in einem Wasserbad inkubiert. Nach dem Hinzufügen von 2 ml Ethanol (100%) wurde das Gemisch über eineQIAamp Midi Säule für 3 Minuten bei 3000 rpm (~1850g) in zwei Arbeitsschritten zentrifugiert. Das Filtrat wurde verworfen und die Säule mit 2 ml Puffer AW1 beladen und für 1 Minute bei 5000rpm zentrifugiert. Entsprechend wurde die Säule mit Puffer AW2 beladen und für 15 Minuten zentrifugiert. Mit 300µl Puffer AE wurde die an die Säule gebundene DNA durch Zentrifugation für 5 Minuten mit 5000 rpm eluiert und in Reaktionsgefäße überführt. Alle Zentrifugationsschritte fanden bei Raumtemperatur statt. Die Lagerung der DNA erfolgte bei –20°C und –80°C.

Die DNA-Konzentration der aufgereinigten Proben wurde mittels UV-Photometrie bestimmt. Der Leerwert wurde mit dem gleichen Mischungsverhältnis wie die Probe, d.h. 5 µl AE-Puffer plus 75 µl H₂O ermittelt. Die DNA-Lösung wurde mit destilliertem Wasser 1:15 in Greiner 96-well Platten verdünnt (5 µl DNA-Lösung, 75 µl Wasser).

Um eine definierte Ausgangskonzentration für die folgenden Arbeitsschritte zu erhalten, wurde aus den Ergebnissen der photometrischen Messung ein Verdünnungsschema berechnet. Die Proben wurden mit einer 1:10-Verdünnung des Elutionspuffers von QIAGEN auf eine einheitliche Konzentration von [Seite 42↓]20 ng/µl gebracht.

Für die Polymerase-Ketten-Reaktion wurden folgende Primer benutzt:

HNMT -(C314T)-SNP

C5aR -(C245T)-SNP

IL4- C-(590T)-SNP

IL13- (Arg130Gln)-SNP

IL13- (C-1055T)-SNP

IL4RA- (Ser478Pro)- und IL4RA- (Gln551Arg)-SNP

Diese Primer schließen beide SNPs ein.

IL4RA- (Ile50Val)-SNP

[Seite 44↓]

Tab. 3: Protokolle für die Single- bzw. Duplex-PCR-Ansätze. Bei allen PCR-Ansätzen wurde mit der Annealingtemperatur = 52°C gearbeitet.

Die PCR lief unter folgenden Konditionen ab:

Initialer Denaturierungsschritt bei 95°C für 10min

PCR-Zyklus:

Es wurde mit 50 Zyklen gearbeitet und nach zweiminütiger Finalelongation auf 4°C abgekühlt.

Um den SNP des Kandidatengens C5aR nachzuweisen, wurde das PCR-Produkt mit folgendem Reaktionsansatz bei 37°C für 90 Minuten restriktionsenzymatisch gespalten.

Die Reaktion wurde mit 1,5µl Probenauftragspuffer gemischt und in Taschen eines 3%-Agarosegels aufgetragen. Die Fragmente wurden gelelektrophoretisch getrennt. Eine Lagerung der Proben zur späteren Verarbeitung wurde bei –20°C vorgenommen.

	Wobei N für eine beliebige Base steht

Abb. 12: Genotypanalyse des C5aR-(C-245T)-SNP verschiedener Individuen mittels PCR-Amplifikation und Restriktionsverdau mit BPU 1103. Das (T245)-Allel wird nicht verdaut. Es resultiert eine Bandenlänge von 301 bp (obere Bande). Liegt das (C245)-Wildtypallel vor, ergeben sich Fragmente mit einer Größe von 258 bp (untere Bande) und 43 bp (nicht zu sehen). Das Enzym BPU 1103 schneidet die Sequenz folgendermaßen:

	Wobei N für eine beliebige Base steht

Jeweils 10 µl des PCR-Produktes beziehungsweise des Restriktionsansatzes wurden auf einer Mikrotiterplatte mit 1 µl Proben-Auftragspuffer gemischt und in die Taschen eines 3%-Agarosegels pipettiert und in TBE-Puffer elektrophoretisch bei 90 V und 400 mA für 60 min aufgetrennt. Als Marker wurde eine 100 bp-DNA-Leiter verwendet. Die Färbung erfolgte in einer Ethidiumbromidlösung (1 µg/ml) für 10 Minuten, anschließend wurde in einem Wasserbad für 10 Minuten entfärbt. Unter einem UV-Transilluminator mit einer Wellenlänge von 312 nm konnten die aufgetrennten Banden durch Fluoreszenz detektiert und photographiert werden.

Der LightCycler besteht aus zwei funktionellen Komponenten, der Thermocycler-Einheit und der Fluoreszenz-Meßeinheit. Das Gerät kann mit Hilfe fluoreszenzmarkierter Oligonukleotidproben sequenzspezifische Fluoreszenz messen. Die Reaktion findet in 32 Glaskapillaren statt, die gleichzeitig als Küvetten für die Fluoreszenzmessung dienen. Eine Hochleistungsleuchtdiode emittiert blaues Licht, das nach Passage eines Filters und eines Kollimators eine Wellenlänge von 470 nm hat und die Proben anregt. Drei unabhängige Photoelemente können Licht der Wellenlängen 530 nm, 640 nm und 710 nm gleichzeitig detektieren.

Um eine Punktmutation nachzuweisen, werden zwei Oligonukleotidsonden eingesetzt. Die erste Sonde, der sogenannte „Sensor“, bindet an den 5’-Abschnitt der zu untersuchenden Region und ist an seinem eigenen 3’-Ende mit dem Fluorophor Fluoreszin gekoppelt. Die zweite Sonde, der „Anchor“, bindet an den 3’-Abschnitt der spezifischen genomischen DNA und ist entweder mit [Seite 47↓]dem Farbstoff LightCyclerRed-640 oder LightCyclerRed-705 an seinem eigenen 5’-Ende gekoppelt. Binden nun beide Sonden an den Zielabschnitt der genomischen DNA, dann liegen die Farbstoffe in unmittelbarer Nähe zueinander. Wird der „Sensor“-Farbstoff durch die Bestrahlung mit der Leuchtdiode angeregt, kann er die aufgenommene Energie an den „Anchor“-Farbstoff weitergeben, der seinerseits angeregt wird und Licht einer anderen Frequenz emittiert, das vom LightCycler detektiert werden kann. Dieser sogenannte Fluorescence Resonance Energy Transfer (FRET) tritt nur auf, wenn beide Sonden an ihre Zielsequenz binden und die Farbstoffe Kopf-an-Kopf zu liegen kommen.

	Abb. 13: Prinzip des FRET bei Hybridisierung der Sonden mit der Zielsequenz und Anregung (nach LightCycler-Manual, Roche).

Das Thermoelement kühlt das Reaktionsgemisch nach einem initialen Denaturierungsschritt bei 95°C auf 45°C ab. Anschliessend wird die Temperatur schrittweise erhöht und zeitgleich die Fluoreszenz gemessen. So kann eine Schmelzkurve aufgezeichnet werden, da die Intensität der Fluoreszenz proportional mit der Ablösung der Sonden von der untersuchten DNA abnimmt. Eine einzige unterschiedliche Base in der untersuchten Zielsequenz (Punktmutation) bewirkt eine schlechtere Hybridisierung mit der Sonde und so eine um einige Grad Celsius verschobene Schmelzkurve.

	Abb. 14: Oben: Die Abbildung zeigt schematisch das Prinzip der Mutationsanalyse mit einer Wildtyp-spezifischen Sonde. Unten: Die HNMT-(C314T)-Genotypisierung zeigt beispielhaft wie das Hybrid aus Sonde und untersuchtem Amplikon durch die „fehlerhafte“ Anlagerung nur einer Base instabiler wird und seine Schmelztemperatur dadurch um einige Grad Celsius niedriger liegt.

Für die Genotypisierung wurden folgende Oligonukleotid-Sonden eingesetzt:

HNMT -(C314T)-SNP

IL4-C -(590T)-SNP

IL13- (Arg130Gln)-SNP

IL13- (C-1055T)-SNP

IL4RA-( Ser478Pro)-SNP

IL4RA- (Gln551Arg)-SNP

IL4RA- (Ile50Val)-SNP

LightCycler Sondenansatz

Zu je 8 µl des PCR-Produktes wurden 8 µl Sondenansatz gegeben, in der Kapillare gemischt und für 30 Sekunden bei 3000 rpm zentrifugiert.

3.3.8  Statistische Methoden

Populationsbasierte Assoziationsstudien

Die Allelfrequenzen für jeden SNP wurden errechnet. Der χ²- Test wurde angewendet, um die Abweichung der Genotypfrequenzen vom Hardy-Weinberg-Gleichgewicht zu bestimmen. Der Expectation-Maximization-Algorithmus des Programms Arlequin (Schneider S, Excoffier L. Arlequin (Version 2.0), Genetics and Biometry Laboratory, 2000)) wurde angewandt, um die Haplotypfrequenzen für jedes SNP-Paar zu bestimmen. Die Kopplungsungleichgewicht-(LD)-Statistik D’ wurde auf der Basis der so [Seite 51↓]ermittelten Frequenzen errechnet [74]. Ein Wahrscheinlichkeitstest des Programms Arlequin wurde durchgeführt, um die Signifikanz von Abweichungen vom Kopplungsungleichgewicht zu untersuchen.

Das Grenzwerteffektmodell Generalized Estimation Equotations (GEE) wurde angewendet, um die Assoziation zwischen Genotyp und den longitudinal gemessenen IgE-Werten zu ermitteln. GEE sind statistische Methoden um Einflussgrößen auf korrelierte Daten mittels Regression zu untersuchen.

Um die Veränderung der IgE-Werte während der siebenJahre für die verschiedenen Genotypen zu testen, wurden Interaktionsterme zwischen Alter, Alter², Alter³ und Genotyp in das Modell aufgenommen.

Für jeden SNP wurden je ein kodominantes Modell und das in früheren Studien verwendete genetische Modell angewendet [75]. Gen-Gen- und Gen-Umweltinteraktionen wurden durch Einfügen von Interaktionstermen in das Regressionsmodell untersucht (z.B. mütterliches Rauchen). P-Werte <0,05 wurden als statistisch signifikant gewertet.

Um den möglichen Einfluss von systematischen Selektionsfehlern auf die Ergebnisse zu testen, wurde ein gewichtetes lineares Modell (weighted linear model (WLM)) angewendet, um die beobachteten Daten bezüglich der ursprünglichen Kohorte zu gewichten. Gewinnt man ähnliche statistische Koeffizienten bzw. Aussagen mit Hilfe des gewichteten linearen Modells und des Grenzwertmodells, ist ein systematischer Selektionsfehler unwahrscheinlich.

Für die Gene mit mehr als zwei Varianten wurden Haplotypanalysen durchgeführt. Die Drei-Lokus-Haplotypanalysen für das IL4RA-Gens wurden mit PHASE (Version 1.0) durchgeführt (Stephens M, 2001). Gewichtete lineare [Seite 52↓]Regression wurde dann für die IgE-Werte zu jedem Zeitpunkt getrennt durchgeführt.

Familienbasierte Assoziationsstudie

Wie bei der populationsbasierten Assoziationsstudie wurde auch hier das Grenzeffektmodell Generalized Estimation Equotations (GEE) verwendet. Es wurde jedoch eine andere Kovariante definiert, die statt des Genotyps des Kindes die Transmission des mutierten Allels von beiden Eltern berücksichtigt.

Der TDT-Test, der auf der Auswertung kompletter Trios (Vater, Mutter, Kind) beruht, wurde im Abschnitt 1.2.2.4 bereits genauer beschrieben. Eine Modifikation des TDT für quantitative Merkmale wurde wie folgt vorgenommen und jeweils auf die longitudinalen Daten der Kinder der MAS-Kohorte angewandt [76]:

Der Grenzeffekt auf das arithmetische Mittel wird durch Kovarianten bestimmt und wird wie folgt errechnet: Log (IgE)ij = β₀ + β₁x_i1 + β₂ (Geschlecht)_i +β₃ (Alter)_ij +β₄ (Alter²)_ij +β₅ (Alter³)_ij + ε_ij. Die Kovarianz wird definiert als: Corr (ε_ij,ε_ik) =R_i (α). Log (IgE)_ij bedeutet der j. IgE-Wert des i. Kindes, und x₁ ist eine Kovariante, die die Transmission von Allel 1 beider Eltern eines Kindes darstellt. Die Transmission nimmt Werte von x_i1 = H_i ^M (T_i ^M – ½) + H_i ^P (T_i ^P – ½) an, wobei H_i ^M und H_i ^P jeweils einen Wert von 1 oder 0 annehmen, abhängig davon, ob die Mutter und der Vater des Kindes heterozygot für Allel 1 ist. Die Transmissionsindikatoren T_i ^M and T_i ^P werden Werte von 1 annehmen, wenn Allel 1 von jeweils Mutter oder Vater auf das Kind transmittiert wurde und 0 wenn nicht.

Alle Analysen wurden mit dem statistischen Paket STATA (Version 7) durchgeführt und die Analysen der Longitudinaldaten sowie die familienbasierten Analysen mit SAS (Version 8.2).