The BTL Procedure (Experimental)

Suppose you have genotyped 20 members of an experimental backcross population at five markers spanning two chromosomes, and you have also recorded the disease resistance status of each subject as resistant or not resistant. You are interested in finding whether there are BTL in the genetic region spanned by your marker set, and if so, where those BTL are and how strongly linked they are to the disease resistance locus. The first step is to input the data, and then to use PROC BTL with the appropriate options to request all single-marker models of the data to be calculated, as in the following program:

data MarkerDat;
   input (m1-m5) ($) trait;
   datalines;
A B B A B 1
B A B B A 1
A B A B B 1
A A B B A 1
A B A B B 1
B A A A B 1
A A B A A 1
B B A B A 1
A A B B A 1
A A A B A 1
A B A A B 0
B A B A A 0
B A A A B 0
B B A B B 0
B A A B B 0
B B A A B 0
A B B A A 0
B A B A A 0
A B A A B 0
B A B A B 0
;

proc btl data=MarkerDat;
   marker m1-m5 /all=1;
   model trait;
run;

The results from the analysis are shown in Figure 4.1.

Similarly, all two-marker models can be calculated as follows, with the results shown in Figure 4.2.

proc btl data=MarkerDat;
   marker m1-m5 /all=2;
   model trait;
run;

Since m1m4 appears to be the best two-marker effect, you can then estimate the recombination and penetrance parameters for this BTL model. First you have to enter the mapping information for the markers as follows:

data MarkerMap;
   input marker $ chromosome position location;
   datalines;
m1 1 1 0
m2 1 2 4.3
m3 1 3 16
m4 2 1 0
m5 2 2 5.5
;

Now you can use the PARMEST statement to request the parameter estimates to be calculated, as in the following code. PROC BTL estimates penetrance values with each recombination parameter set to 0.5.

proc btl data=MarkerDat map=MarkerMap;
   marker m1 m4 /group=chromosome;
   model trait;
   parmest cross=b gen=1 r=0.5;
run;

Figure 4.3 displays information about the model that includes the two-marker effect m1m4. The "Parameter Estimates" table shows that penetrance values are not in the valid range (between 0 and 1) for this model with the given values of .