|
|
Genetic Analysis of Variation in Neuron Number Richelle Cutler
Strom
Chapter 3: Mapping Genes
Controlling Brain Weight Variation Using Recombinant Inbred strains and F2
Intercross Progeny
Introduction
The high heritability of brain weight together with the estimate of a
relatively small number of genes affecting this trait, indicate that it
should be possible to map at least several responsible for variation in
brain weight. In this study, I map quantitative trait loci that are
responsible for variation in brain weight between mouse strains using the
recombinant inbred (RI) strains, BXD and AXB/BXA, and an F2
intercross.
Belknap and colleagues (1992) used the BXD RI strains to map QTLs
responsible for brain weight differences between DBA/2J and C57BL/6J
(Belknap, 1992). In the Belknap study, only 20 RI strains were used, whereas
my BXD analysis includes 26 strains. Also, Belknap’s study uses 360 markers,
whereas I use a data set that includes 529 markers. The increase in both
strain number and marker density greatly improves the power of detecting
significant associations. In addition, this study combines two sets of RI
strains and an F2 intercross for the purpose of
QTL confirmation.
The same genetic interval could be linked to variation in brain weight in
more than one cross. Identifying the same QTLs is especially likely in this
study since there are some common parental origins between the crosses,
ABXD5F2, AXB/BXA, and BXD. For example, ABXD5F2
share the A/J parental with the AXB/BXA RI strains, DBA/2J is a parental for
the BXD strains, including BXD5, which is a parental in ABXD5F2.
In addition, all of the crosses have C57BL/6J origins in common. The
detection of the same linkage interval in multiple crosses will provide QTL
verification. I also expect to identify unique QTLs within each cross since
there will be different allele segregating between the crosses.
There are several advantages in using RI strains for QTL analysis. First,
since the individuals of each RI strain are isogenic an infinite number of
mice can be phenotyped to obtain a strain average. Obtaining a strain
average serves to reduce the environmental and technical variance, that can
obscure the association between a phenotype and a genotype. Second, with RI
strains there is no genotyping necessary since there are hundreds of marker
loci already typed for each RI strain. Another advantage with RI strains is
that during the generation of each RI strain, there were multiple
opportunities for recombination, resulting in a 4-fold increase in resolving
power for determining linkage distance.
The F2 intercross has the advantage of being able to detect
genes with a lower effect size, which are not detected in an RI analysis due
to the limited number of strains in a set (Johnson et al., 1992). However,
detecting and positioning QTLs with reasonable accuracy requires genotyping
and phenotyping a large number of intercross progeny. The number of progeny
necessary to genotype can be reduced by first phenotyping and then selecting
the animals from the phenotypic extremes for genotyping (Darvasi, 1997;
Lander and Botstein, 1989). This technique is called selective genotyping
and in addition to saving time, the progeny with phenotypic extremes will
have parental genotypes, which are the most informative for linkage
detection (Tanksley, 1993).
Materials and methods
Mapping Progeny
The BXD set consisted of 26 strains that were generated from parental
strains, DBA/2J and C57BL/6J, which have a 73 mg difference in brain weight.
All the BXD mice were obtained directly from the Jackson Laboratory (Bar
Harbor, ME). The AXB, BXA RI strains were produced from reciprocal crosses
between the inbred strains A/J and C57BL/6J. The strains A/J and C57BL/6J
have an 80 mg difference in brain weight. Brain weights were obtained for 24
AXB/BXA strains. Mating pairs for each of the 24 AXB/BXA RI strains were
obtained from the Jackson Laboratory and the majority of the mice were
generated in the mouse colony at the University of Tennessee. The ABXD5F2
progeny were generated from a cross between A/J and BXD5 (n = 542).
The strains A/J and BXD5 have a 106 mg difference in brain weight. The BXD5
strain has one of the heaviest brain weights of all inbred strains—528 mg,
while A/J has one of the smallest—422.5 mg.
Corrected brain weights
Fixed brain weights were obtained as described in Chapter 2. The effects
of the variables sex, age, body weight, litter size and parity on brain
weight were determined by multiple regression as in Chapter 2. Standardizing
brain weights to the average body weight, age, litter size, parity, and
female sex can minimize the effects of these variables on brain weight.
Mapping with standardized or corrected brain weights, increases the
probability of mapping only brain weight QTLs. The BXD brain weights were
standardized to that of a 21.5 gm, 79-day-old female mouse (the body weight
and age average across BXD strains), using the regression coefficients for
body, age, and sex shown in Table 2.10. Litter size and parity data were not
available for the RI strains. The sex code for a female mouse is 1 and is 0
for a male mouse. For example, the equation to standardize brain weight
would be:
Corrected brain weight = fixed brain weight + 4.43*(21.5 – 20) +
10.44*(1–0) + –0.18*(79 – 40).
In the equation, the coefficients for body weight, sex, and age, are
italicized. Corrected brain weight averages for the 26 BXD strains can be
found in Table 3.1. Similarly, AXB/BXA brain weights were standardized using
their respective coefficients, to a 19.5 gm, 78.8-day-old female. Corrected
brain weight averages for the 24 AXB/BXD strains can be found in Table 3.2.
The ABXD5F2 brain weights were standardized to
that of a 40-day-old, 19 gm, female mouse, coming from a litter size of 10,
and born in a litter with a parity of three.
Table 3.1. Corrected brain weight and body weight
data for BXD strains and parentals, C57BL/6J and DBA/2J.
Strain |
Brain wt. (mg) ± SE |
n |
Body wt.(gm) |
Mean age(days) |
C57BL/6J |
487.0 ± 2.8 |
72 |
24 |
108 |
DBA/2J |
424.7 ± 3.7 |
35 |
23 |
149 |
BXD1 |
465.7 ± 3.7 |
16 |
21 |
71 |
BXD2 |
424.3 ± 4.0 |
17 |
23 |
62 |
BXD5 |
527.6 ± 5.5 |
20 |
22 |
123 |
BXD6 |
419.7 ± 5.9 |
12 |
22 |
193 |
BXD8 |
460.7 ± 5.9 |
29 |
23 |
170 |
BXD9 |
432.2 ± 3.3 |
21 |
24 |
82 |
BXD11 |
442.2 ± 4.4 |
10 |
21 |
66 |
BXD12 |
436.2 ± 5.6 |
14 |
23 |
116 |
BXD13 |
406.1 ± 3.4 |
21 |
22 |
57 |
BXD14 |
442.5 ± 3.3 |
14 |
24 |
85 |
BXD15 |
436.0 ± 3.6 |
10 |
26 |
129 |
BXD16 |
462.5 ± 6.4 |
11 |
21 |
61 |
BXD18 |
421.6 ± 2.4 |
15 |
23 |
73 |
BXD19 |
434.1 ± 4.0 |
22 |
19 |
88 |
BXD20 |
407.4 ± 3.2 |
20 |
19 |
74 |
BXD21 |
457.6 ± 5.8 |
10 |
25 |
114 |
BXD22 |
450.4 ± 3.3 |
16 |
23 |
84 |
BXD23 |
439.8 ± 3.5 |
16 |
20 |
83 |
BXD24 |
399.7 ± 4.7 |
16 |
24 |
189 |
BXD25 |
433.9 ± 2.6 |
19 |
18 |
77 |
BXD27 |
382.6 ± 4.2 |
32 |
21 |
99 |
BXD28 |
408.9 ± 2.6 |
33 |
22 |
68 |
BXD29 |
398.7 ± 5.6 |
18 |
20 |
73 |
BXD30 |
392.8 ± 3.4 |
18 |
18 |
80 |
BXD31 |
417.4 ± 5.9 |
11 |
23 |
75 |
BXD32 |
432.7 ± 3.0 |
55 |
26 |
129 |
Table 3.2. Corrected brain weight and body weight
data for AXB/BXA strains and parentals, A/J and C57BL/6J.
Strains |
Brain wt. (mg) ± SE |
n |
Body wt. (gm) |
Mean age (days) |
C57BL/6J |
487.0 ± 2.8 |
72 |
24 |
108 |
A/J |
422.5 ± 3.9 |
25 |
22 |
190 |
BXA1 |
450.3 ± 5.3 |
19 |
20 |
56 |
BXA2 |
455.2 ± 3.0 |
19 |
20 |
67 |
BXA4 |
417.4 ± 3.5 |
18 |
19 |
80 |
BXA7 |
457.9 ± 4.2 |
15 |
22 |
79 |
BXA8 |
436.5 ± 3.1 |
14 |
21 |
80 |
BXA11 |
465.2 ± 6.5 |
14 |
21 |
72 |
BXA12 |
464.1 ± 4.4 |
19 |
17 |
81 |
BXA13 |
463.3 ± 3.1 |
38 |
17 |
56 |
BXA14 |
414.7 ± 3.4 |
19 |
23 |
66 |
BXA16 |
437.7 ± 4.1 |
17 |
19 |
73 |
BXA17 |
462.8 ± 6.4 |
5 |
25 |
134 |
BXA24 |
436.5 ± 3.9 |
15 |
21 |
80 |
BXA25 |
474.7 ± 16.1 |
2 |
24 |
176 |
BXA26 |
430.7 ± 3.6 |
21 |
19 |
67 |
AXB1 |
428.4 ± 4.0 |
22 |
20 |
56 |
AXB2 |
475.9 ± 6.3 |
10 |
21 |
92 |
AXB4 |
450.9 ± 3.0 |
4 |
19 |
73 |
AXB5 |
460.7 ± 6.1 |
19 |
20 |
71 |
AXB6 |
467.9 ± 4.9 |
11 |
23 |
93 |
AXB8 |
441.8 ± 2.4 |
16 |
18 |
67 |
AXB10 |
475.0 ± 4.6 |
18 |
20 |
79 |
AXB12 |
414.1 ± 3.1 |
14 |
17 |
81 |
AXB13 |
480.9 ± 6.5 |
4 |
21 |
129 |
AXB15 |
430.1 ± 2.4 |
23 |
19 |
58 |
AXB18 |
468.4 ± 4.4 |
9 |
18 |
89 |
AXB19 |
466.6 ± 3.0 |
22 |
21 |
69 |
AXB20 |
463.2 ± 4.4 |
15 |
18 |
85 |
AXB24 |
440.4 ± 4.5 |
19 |
19 |
77 |
Genotyping
The BXD and AXB/BXA RI strains have already been genotyped at
approximately 600 and 700 loci, respectively. (Marshall et al., 1992;
Sampson et al., 1998). For the ABXD5F2 progeny, I
performed a genome-wide screen by typing 126 microsatellites dispersed at an
average interval length of 30 cM across the entire genome, Chr Y excepted
(Fig. 3.1). Genomic DNA was isolated from spleens of the intercross progeny
using the simplified DNA isolation procedure described by Laird et al.
(1991), with minor modifications. Details of the genotyping procedure are
available on the WWW at
(http://nervenet.org/papers/PCR.html).
Microsatellites were typed using the polymerase chain reaction (PCR). The
PCR protocol was set up as described by Dietrich et al. (1992). A PCR
reaction contained (10 µl): 1X PCR buffer, 1.92 mM MgCl2, 025 U/µl Taq
polymerase (0.5 U), 0.2 mM dNTPs, 132 nM Mit map pairs and 50 ng genomic
DNA. The map pairs were purchased from Research Genetics, Huntsville, AL. A
loading dye (60% sucrose, 1.0 mM cresol red, 11% in-assay) was added to the
reaction before the PCR. Polymerase chain reactions were carried out in
96-well microtiterplates, which allowed the use of a 12-channel multichannel
pipettor to dispense reagents and load gels. I used the high stringency
thermocycling protocol called "touchdown" described by Don et al. (1991).
The protocol is 94° for 3 min., followed by 5 cycles of (94° for 30 sec,
60°C and then –2°C with each cycle, and 72° for 2 min). The protocol then
continues with 30 cycles of 94° for 30 sec, 50° for 2 min. and 72° for 2
min., with a final cycle at 72° for 7 min.
Figure 3.1. Microsatellite marker positions for ABXD5F2
intercross screening.
The PCR products were separated by electrophoresis in 3% Metaphor agarose
gels (FMC Bioproducts, Rockland ME.) in 1X TBE buffer. Gels were generally
run for 90 min at 170 V, and maintained below 30 °C with cooled circulating
buffer. PCR products were visualized by in-gel staining with ethidium
bromide (Fig 3.2). Most loci could be scored as homozygous or heterozygous
with complete confidence. To detect QTLs in the ABXD5F2
progeny, I genotyped ~4% of with the lowest corrected brain weights and ~4%
with the highest corrected brain weights (the lowest 22 cases and the
highest 23 cases). These 45 extreme cases were chosen because with the
addition of the two parentals and F1 DNA the
number of DNA samples equaled 48, which used exactly half a microtiter
plate, and allowed the testing of two markers per plate. To determine the
effect size of a particular QTL, all the F2
progeny were genotyped.
Figure 3.2.
PCR amplified microsatellite polymorphic between A/J and BXD5 and
separated by electrophoresis on a 3% agarose gel.
Linkage analysis for a quantitative trait involves comparing the trait
means of the different marker classes, i.e., parental and recombinant
genotypes AA, AB, BB. A significance difference between trait values between
the different classes indicates a linkage between the genetic marker and the
QTL and this can be tested by an analysis of variance (Kearsey and Pooni,
1996). I used the computer program Map Manager QT to perform linkage testing
(Manly, 1996). Map manager QT reports the measure of linkage between a test
interval and the phenotype with a likelihood ratio statistic (LRS). The LRS
is computed by taking the n loge of the ratio of the sum of
squares (SS) of the groups over the SS of the total regression. A linkage
statistic is commonly reported with a LOD score (logarithm of the odds). The
LOD score represents the log10 of the ratio of two probabilities,
that of linkage, and that of no linkage. The LRS statistic is essentially
equivalent to a LOD score and can be converted by dividing it by 4.61.
The ability to detect and positionally map a QTL depends on how close the
marker is to the QTL. When regions of the map are not adequately covered
with markers a procedure called interval mapping can increase the ability to
detect and accurately position a QTL. Interval mapping uses simple linear
regression equations to predict additive phenotypic deviation at intervening
positions between flanking markers. The predicted additive or dominant
deviation coefficients are computed across the interval and compared to the
observed additive deviation coefficients and at each position a residual sum
of squares (RSS) is calculated (Haley and Knott, 1992). A QTL effect and
positional assignment can be seriously skewed by the presence secondary QTLs
outside the QTL interval. To eliminate this interference I performed
composite interval mapping, in which QTLs mapping outside the interval were
fitted as cofactors in a multiple regression analysis (Jansen, 1993; Zeng,
1994). I used Map Manager QT to perform composite interval mapping.
Significance levels
In testing multiple markers in a genome-wide scan one will obtain
significance associations just by chance. Therefore, it is necessary to
determine the probability of a significant linkage in genome-wide test. The
genome-wide probability thresholds can be determined with a permutation
test. A permutation test shuffles the trait data, reassigning the data to
new individuals and then compares this "null distribution" with the
genotypic pattern (Churchill and Doerge, 1994). The trait data is then
reshuffled from 1,000 to 10,000 times and after each permutation the single
best linkage statistic is recorded. The genome-wide significance thresholds
are computed from the incidence of test statistics in the complete
permutation test. For example, if the highest linkage statistic was reached
only once in 1000 permutations, then the probability of obtaining that high
a linkage statistic by chance is 1 in 1000. The statistics obtained with the
experimental data are compared to the distribution of permutation scores to
determine their genome-wide probability. The false positive rates of 0.63,
0.05, and 0.001 correspond to a suggestive, significant, and highly
significant linkage. The advantage of using permutation tests is that the
threshold values are specifically tailored to the parameters of a given data
set, e.g., chromosome number, genetic length, genetic marker density, and
trait distribution. I used the program Map Manager QT to perform the
permutation tests (Fig. 3.3). I computed conditional empirical thresholds (CET)
to control for the effects of secondary QTLs (Doerge and Churchill, 1996).
Figure 3.3. Permutation of BXD brain weight data for conditional
genome-wide significance testing. The BXD
data set was permuted 20,000 times at 1 cM intervals across the genome with
conditional loci, R1Rik100 and Iapls2-13. The suggestive LOD,
(p =0.63), is 1.2, the significant LOD, (p =0.05), is 3.1 and
the highly significant LOD, (p =0.001), is 5.7.
Results
Recombinant inbreds
The best match with variation in brain weight was found at D11Mit53.
The correlation coefficient between brain weight and alleles, assigned
either a 1 and 0 at D11Mit53, is +0.90 (Table 3.3). The underlined
alleles indicate discordance with the brain weight phenotype. In Table 3.3,
the allele assignments are based on the strain’s average brain weight
resemblance to that of C57BL/6J (B) or DBA/2J (D). The U stands for unknown
genotype and is assigned for phenotypes between that of the parentals. There
are five strains that have intermediate genotypes. Refer back to Fig. 2.3A
for the parental brain weight distribution.
A linkage analysis with corrected BXD brain weights resulted in a LOD
score of 4.7 at D11Ncvs74 at 0 cM on proximal Chr 11 and at 12 cM a
LOD score of 4.5 at D11Mit53. Interval mapping on Chr 11 increases
the LOD score to 5.2 between Glns-ps1 and D11Mit53.
Glns-ps1 and D11Mit53 map to 8.6 cM and 12.1 cM, respectively.
The concordance of brain weight variation and alleles at Glns-ps1 and
D11Mit53 is shown in Table 3.3. To map secondary QTLs that may not be
detected due to the interference of D11Mit53, I tested for linkage
while controlling for the variance associated with D11Mit53. This
resulted in detecting secondary QTLs on distal Chr 1 at D1Rik100 with
an LOD of 2.8 and on mid Chr X at Iapls2-13 with LOD 2.8.
Table 3.3. BXD strain distribution pattern for
brain weight and genotypes at five loci on Chr 11
Locus |
cM |
BXD strains |
|
|
1 |
2 |
5 |
6 |
8 |
9 |
11 |
12 |
13 |
14 |
15 |
16 |
18 |
19 |
20 |
21 |
22 |
23 |
24 |
25 |
27 |
28 |
29 |
30 |
31 |
32 |
D11Mit2 |
5.4 |
B |
B |
B |
D |
B |
D |
B |
B |
D |
D |
B |
B |
D |
B |
D |
B |
B |
D |
D |
B |
D |
D |
D |
B |
D |
D |
D11Ncvs73 |
6.4 |
B |
B |
B |
D |
B |
D |
B |
B |
D |
D |
D |
B |
D |
B |
D |
B |
B |
D |
D |
B |
D |
D |
D |
B |
D |
D |
Glns-ps1 |
8.6 |
B |
D |
B |
D |
B |
D |
B |
B |
D |
D |
D |
B |
D |
B |
D |
B |
B |
D |
D |
B |
D |
D |
D |
D |
D |
D |
Brain
weight |
|
B |
D |
B |
D |
B |
D |
U |
U |
D |
U |
U |
B |
D |
D |
D |
B |
B |
B |
D |
U |
D |
D |
D |
D |
D |
D |
D11Mit53 |
12.1 |
B |
D |
B |
D |
B |
B |
B |
B |
D |
D |
D |
B |
D |
D |
D |
B |
B |
B |
D |
B |
D |
D |
D |
D |
D |
D |
D11Mit51 |
15.5 |
D |
D |
B |
D |
B |
B |
B |
B |
D |
D |
D |
B |
D |
B |
D |
B |
B |
B |
B |
B |
D |
D |
D |
D |
D |
D |
Composite interval mapping on Chr 11 with the two secondary loci as
cofactors resulted in an LOD score of 5.1 at D11Ncvs74 and a peak LOD
score of 8.4 between Glns-ps1 and D11Mit53, at 12 cM (Fig.
3.4). A 2 LOD confidence interval (CI) will have a 99% probability of
containing the QTL. The 2 LOD CI for the QTL linked to D11Mit53 is
between 11–16 cM. Correcting for the variance of the secondary loci
increases the ability to detect QTLs and to accurately map and estimate QTL
effects. A permutation performed on this data set defined the conditional
empirical threshold for a highly significant linkage to be 5.7 (Fig. 3.3).
The brain weight QTL mapped to proximal Chr 11 at 12 cM, exceeds the highly
significant threshold. Thus, there is less than 0.001 probability that the
linkage is spurious in a genome-wide scan! This QTL has been named brain
size control 1 (Bsc1).
Figure 3.4. Linkage map demonstrates the QTL Bsc1 on Chr 11
in the BXD data set.Linkage statistics were computed at 1 cM intervals while
controlling for D1Rik100 on distal Chr 1, and Iapls2-13 in
middle of the X Chr.
Linkage analysis in the AXB/BXA strains indicated that a gene modulating
brain weight is close to a marker on distal Chr 7. Using simple linear
regression, the peak LOD on distal Chr 7 at D7Mit71 was 3.8. The
correlation coefficient for alleles at D7Mit71 and brain weight is
+0.72. Controlling for this QTL interval on Chr 7 indicated a linkage on
proximal Chr 14 (D14Mit14 and D14Mit48), with an LOD of 3.0
using interval mapping. Composite interval mapping on Chr 7, while
controlling for the proximal Chr 14 region increased the LOD score on distal
Chr 7 at D7Mit71, to 6.7 (Fig. 3.5). Permuting the trait data 10,000
times with the conditional locus on Chr 14, resulted in only 1 in 1000 tests
as significant as the strain distribution pattern at D7Mit71. The QTL
mapped to proximal Chr 7 at 65.2 cM, with a 2 LOD CI between 64-72 cM,
exceeds the highly significant threshold! This QTL is named Brain size
control 2, (Bsc2).
Figure 3.5. Linkage map demonstrates the QTL Bsc2 on Chr 7
in the AXB/BXA data set. Linkage statistics were computed at 1 cM intervals
while controlling for loci D14Mit14 and D14Mit48.
ABXD5F2 Intercross progeny
An association was found between brain weight variation in the 45 extreme
ABXD5F2 progeny and their genotypes at two loci on Chr 14.
Linkage testing resulted in LOD scores on Chrs 14 of 5.8 at D14Mit260
and 5.4 at D14Mit69. Three secondary QTLs were identified, one on mid
Chr 8 at D8Mit242 with an LOD score of 4.3, one on distal Chr 5 at
D5Mit291 with an LOD score of 3.6, and one on proximal Chr 18 at
D18Mit36 with a LOD score of 3.0. Controlling for the Chr 14 interval
increased only the LOD score at Chr 18 at D18Mit36 to 3.7. A
permutation test defined the LOD score required for a significant linkage to
be 4.7. Thus, the linkage at D18Mit36 is only suggested. Composite
interval mapping while controlling for D18Mit36 boosts the LOD score
at D14Mit260 and D14Mit69 to 6.8 and 5.9, respectively (Fig.
3.6). The QTLs on Chr 14 map to 25 cM with a 2 LOD CI between 20–30 cM and
to 59.1 cM, with a 2 LOD CI between 50–69 cM. Permuting the ABXDF2
trait data 5000 times with the conditional locus on Chr 18, resulted in only
0.7% and 1.6% of the tests with LOD scores as high as D14Mit260 and
D14Mit69, respectively. Thus, these two QTLs exceed the significant
threshold, p < 0.05 and have been named Brain size control 3 (Bsc3)
and Brain size control 4 (Bsc4).
In the AXB/BXA cross, a secondary QTL was detected on proximal Chr 14
between markers D14Mit14 and D14Mit48, loci that map to 0 cM
and 13 cM, respectively. The peak linkage statistic was at 9 cM. The Chr 14
committee report places D14Mit48 a little more distal at 16.5 cM. The
Chr 14 QTLs in the ABXD5F2 cross, D14Mit260
and D14Mit69, map to 25 cM and 59.1 cM, respectively. These are close
to the positions reported by the Chr 14 committee, 29 cM and 43 cM,
respectively. The secondary QTL identified in AXB/BXA strains is too
proximal to be the Bsc2 in ABXD5F2.
However, the secondary QTL in AXB/BXA does give credence to the small
increase in the linkage statistic proximal to the D14Mit260, between
markers D14Mit98 and D14Mit144 (Fig. 3.6).
Figure 3.6. Linkage map demonstrates the QTLs Bsc3 and
Bsc4 on Chr 14 in the ABXD5F2 data set.
Linkage statistics were computed at 1 cM intervals while controlling for the
locus D18Mit17.
QTL effect size
To examine the effect size of QTLs on Chr 14 and Chr 18 in the ABXD5F2
cross I performed a regression analysis on the brain weight z-scores for
each genotype class. A single allelic substitution at the QTL linked to
D14Mit260 is responsible for an effect size of 0.3 SD units, with a 95%
confidence interval equal to ± 0.21. An analysis of variance determined that
the regression was highly significant ( p < 0.0001) and that variance
at D14Mit260 accounts for 5.9% of the variance in brain weight. The
standard deviation of the corrected brain weights for ABXD5F2
is 21.7. Thus, each allelic substitution at Bsc3 translates into ~6.5
mg difference in brain weight. The effect of a single allelic substitution
at the QTL near D14Mit69 is also equal to 0.3 SD ± 0.23, explains
5.1% of the variance, and each allelic substitution translates into ~ 6.5 mg
difference in brain weight. The effect of a single allelic substitution at
the QTL near D18Mit20 is equal to 0.25 SD ± 0.23, explains 3.5% of
the variance, and produces a ~5.4 mg difference in brain weight. The
heritability estimate for ABXD5F2, reported in
Chapter 2, is 24%. Adding the explained variance of the three QTLs equals
14.5%, therefore, together they explain 60% of the genetic variance and 15%
of the total phenotypic variance in brain weight among ABXD5F2
progeny.
Discussion
I have mapped four QTLs that have major effects on brain weight in mouse
to Chrs 7, 11, 14. These brain weight QTL have been named Bsc1, 2, 3,
and 4. Bsc1 maps proximally to Chr 11 at 12 cM, Bsc2
maps distally to Chr 7 at 65.2, and Bsc3 and Bsc4 map to Chr
14 at 25 cM and 59.1 cM, respectively. Secondary brain weight QTLs were
mapped to Chrs 1, 5, 8, 11, 14, 18, and X.
In Chapter 2, a correlation between albino coat color and brain weight in
the AXB/BXA strains indicated that a QTL controlling brain weight mapped
near the albino allele on distal Chr 7. The significant linkage between
brain weight and D7Mit71 on distal Chr 7 in the AXB/BXA cross
confirms the linkage between coat color and brain weight.
A secondary QTL was mapped to the X Chr in the BXD set. I looked for
evidence of a sex-difference in brain weight among reciprocal C57BL/6J and
DBA/2J F1 males, since they inherit their X Chr
from different strains. There was no difference in brain weight between
reciprocal hybrid males. However, this comparison was not paired; there were
6 B6D2F1 males and 18 D2B6F1
males. Regardless the similarity of reciprocal F1
males does not preclude a sex association because detection could be
confounded by autosomal variation if the effect size is too small.
A study by Belknap et al. (1992) mapped suggestive QTLs controlling brain
weight in BXD strains to Chrs 2, 6, 7, 11, and 13. The locus Belknap
identified on Chr 7 was proximal around 14 cM, while my AXB/BXA analysis
detected a QTL on distal Chr 7. The loci identified on Chr 11 by Belknap are
proximal, including Pmv-22, Glns, and Hba, which map
between 8 cM and 16 cM. In my BXD mapping analysis, the best linkage to
brain weight variance was between 11 cM and 16 cM, at Glns-ps1 and
D11Mit53, respectively, with the peak LOD score at 12 cM. Therefore, I
have confirmed the Chr 11 QTL reported by Belknap et al. (1992). There are
several explanations why I did not confirm more of the QTLs detected in the
Belknap et al. (1992) study. First, the Belknap study used only 20 strains,
whereas I use 26. Second, in the Belknap study the brain weights were
regressed with respect to body weights by computing brain to body weight
ratios, rather than a regression analysis. And third, the significance of
Belknap's QTLs are only suggestive since they are nominal p values, whereas
mine are genome-wide significant.
Effect size
Hundreds of genes are probably involved in controlling neuron number,
however, only a dozen or so contribute to the variation in cell number. Many
of these QTLs will have small effects on global neuron number while some
will are large effects. It is the genes with large effects that will be
detected in a linkage analysis, and thus, the four brain weight QTLs are
expected to have relatively large effect size. The two QTLs mapped on Chr 14
in the ABXD5F2 cross, account for 5.9% and 5.1% of
the phenotypic variance in brain weight. In addition, the suggestive QTL
mapped in the ABXD5F2 cross on Chr 18, accounts
for 3.5% of the phenotypic variance in brain weight. The explained variance
estimates are not likely to be overestimates since the sample size of the
ABXD5F2 cross (n =522) is large enough for
reliable estimates of effect size. Thus, in the ABXDF2
cross, three major effect QTLs account for 60% of the genetic variance and
contribute to ~15% of the phenotypic variance in brain weight.
Surprisingly, significant QTLs were not identified in multiple crosses.
Thus, each cross has unique major effect brain weight QTLs segregating,
despite some common origins of the parental strains. In addition to the two
major brain weight QTLs mapped to Chr 14 in the ABXD5F2
cross, a secondary QTL on proximal Chr 14 was mapped in the AXB/BXA strains.
The detection of three QTLs for brain weight on Chr 14 indicates that brain
weight QTLs cluster on Chr 14.
Candidate Genes
The complexity of a quantitative trait can limit the precision to which a
QTL can be mapped, making the direct cloning of the QTL intractable.
However, due to the large number of genes that are already cloned, in many
cases, positional cloning will not be necessary. One can simply search for
candidate genes that map within the QTL interval. Candidates would have a
compatible function or expression pattern, or, when mutated, disturb the
phenotype under study. For instance, the tubby gene was a candidate
for a QTL controlling diet-induced obesity susceptibility in mice because a
mutation in tubby, which mapped to the same region as the QTL, caused
a severe mouse obesity phenotype (Kleyn et al., 1996). Most QTLs have been
identified fortuitously in this way, whereby a gene is recognized by a major
morphological mutation and genetic variants among strains are subsequently
found.
A search for brain weight QTL candidates was performed using the Portable
Dictionary of the Mouse Genome developed by Robert Williams and attainable
through the web (http://mickey.utmem.edu/front.html)
(Williams, 1994), and from the Mouse Genome Informatics (MGI) database,
located at the Jackson Laboratory (http://www.informatics.jax.org/).
The human genome was also inspected for candidate genes in mouse conserved-synteny
regions using a homology search in the MGI database. QTL candidates were
selected if they mapped within the QTL interval, are expressed in the brain,
and are known to be involved in cell proliferation, or when absent or
mutated disrupt brain development. Two candidates for the secondary QTL that
mapped to the proximal end of Chr 11, between 0 and 1 cM, are the mouse
homolog of the human gene neurofibromatosis type 2 (NF2) and the beta
subunit of brain-specific Ca2+/calmodulin-dependent protein kinase II (Camk2b).
The human NF2 gene encodes a tumor suppressor protein called
schwannomin. A mutation in NF2 causes an autosomal dominant disorder
that increases the chance of developing benign tumors in the nervous system.
The expression of NF2 is found on embryonic days 15 and 16 in cells
that are migrating from the ventricular zone to the cortical plate. NF2
is likely to play a role in controlling cell proliferation, perhaps in onset
of differentiation. The candidate Camk2b is first expressed in the
brain on E12.5 and continues to be expressed in adult brain. The specific
role that Camk2b performs in the brain is unknown.
A number of potential Bsc1 candidates map within the 2 LOD CI
between 11 cM and 16 cM on Chr 11. These include RAB1, a member of
the RAS oncogene family, Krct, a member of a novel subfamily of
serine/threonine kinases, Otx1, creatine kinase brain, related
sequence 5 (Ckb-rs5), epidermal growth factor receptor (Egfr),
and epidermal growth factor receptor-related sequence, (Egfr-rs).
All of these genes are expressed in the brain. Otx1 is an attractive
candidate because of its role in cortical neurogenesis. Otx1 is
expression in mouse begins at E10 in the rostal brain and continues until
late gestation. Mice homozygous null for Otx1 have defective
sense-organ development and neurogenesis (Suda et al., 1999). Egfr
maps to 9.0 cM, a little too proximally for Bsc2, but was thought to
be an attractive candidate because of EGF’s effect on cell proliferation
(see Introduction–developmental control of neuron number). However,
EGFR has been dismissed as a candidate because the average brain
weights from mice carrying a mutation in EGFR called waved-2 (wv2)
are the same as wildtype brains (wv2– 488 ± 20.7, n =4 and
wt– 485 ± 9.0, n =12). The mutant wv2 has a single point
mutation in EGFR that produces a 5–10 fold decrease in
autophosphorylation and over a 90% reduction in substrate phosphorylation (Leutteke
et al., 1994).
The best Bsc2 candidate is Msx3. The gene Msx3 maps
to 69 cM on Chr 7 within the 2 LOD CI for Bsc2 and is expressed in
the ventricular zone of the dorsal neural tube during early development
(Wang et al., 1996). The best Bsc3 candidate is orthodenticle,
Otx2. Otx2 maps to 20.5 cM on Chr 14, is expressed in the
developing central nervous system, and plays an important role in regional
morphogenesis of the anterior brain structures (Suda et al., 1999).
Unfortunately, I did not find any candidate genes for Bsc4 on distal
Chr 14.
This brain weight mapping study has identified the location of four brain
weight QTLs. In addition a survey of candidate has suggested some exciting
candidates for brain weight QTLs. In particular three homeobox containing
genes, Msx, Otx1 and Otx2. Homeobox genes have been
suggested to function in some instances, as a regional proliferation index.
Furture studies could involve testing some of these candidate genes.
However, narrowing the QTL interval will probably be necessary to determine
the best candidate and proceed with cloning.
It will be exciting to discover the identity of the genetic variants and
the mechanism by which they influence brain weight. In addition, the
influence of these genetic variants in the variation of brain weight in
other species can be tested to determine if there are common evolutionary
mechanisms, as expected.
Return
to Table of Contents |
|
|