Select      
 Site search   
  Home    Publications

Browse Publications
 
List of Contents

Quantitative Neurogenetics & QTL Mapping

Genetics of Myopia

Control of Neuron Number and Stereology

Growth Cones and Dying Axons

Retina Development and Visual System Mutants

Grant Application

U.S. Patent

Abstracts


Need Help?
Help with Publications
Help with Nervenet
Contact Us

     

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

   


Neurogenetics at University of Tennessee Health Science Center

Text Only Top of Page

Home Page  |  Genome DBs  |  Phenome DBs  |  Publications  |  People & Associates
Mouse Brain Library  |  Related Sites  |  Complextrait.org

Nervenet.org  |   MBL.ORG

Robert W. Williams | Alex Williams © 2002, Nervenet.org