
Cancer Research 67, 2508, March 15, 2007. doi: 10.1158/0008-5472.CAN-06-3157
© 2007 American Association for Cancer Research
Molecular Biology, Pathobiology, and Genetics |
Fine Mapping and Candidate Gene Analyses of Pulmonary Adenoma Resistance 1, a Major Genetic Determinant of Mouse Lung Adenoma Resistance
Min Wang,
Zhongqiu Zhang,
Zhuo Zhang,
Haris Vikis,
Ying Yan,
Yian Wang and
Ming You
Department of Surgery and The Alvin J. Siteman Cancer Center, Washington University School of Medicine, St. Louis, Missouri
Requests for reprints: Ming You, Department of Surgery and The Alvin J. Siteman Cancer Center, Washington University School of Medicine, Box 8109, 660 Euclid Avenue, St. Louis, MO 63110. Phone: 314-362-9294; Fax: 314-362-9366; E-mail: youm{at}msnotes.wustl.edu.
 |
Abstract
|
|---|
Pulmonary adenoma resistance 1 (Par1) is a major genetic determinant of mouse lung adenoma resistance. Although Par1 was previously mapped to mouse chromosome 11 by conventional linkage analyses, its candidate region was broad and undefined. In our present study, we generated Par1 congenic mice using two mouse strains A/J (Par1/) and Mus spretus (Par1/+). Analyzing these congenic mice enabled us to fine map the Par1 quantitative trait loci (QTL) into a 2.0-cM (2.2 Mb) chromosomal region between genetic marker D11Mit70 and the gene Hoxb9. We then conducted systematic candidate gene screening through nucleotide polymorphism and expression analyses. Genes showing differential lung tissue expression or carrying nonsynonymous single nucleotide polymorphisms were identified and discussed. In particular, we evaluated tumor suppressor gene Tob1 for its Par1 candidacy. Our findings have narrowed the Par1 QTL region and will greatly facilitate the identification of the major genetic determinant of mouse lung adenoma resistance. [Cancer Res 2007;67(6):250816]
 |
Introduction
|
|---|
Lung cancer is the leading cause of cancer death in men and women in the United States (1). Although tobacco smoking is the major risk factor, there is strong evidence for genetic susceptibility and gene-environment interactions in lung cancer development (25). However, genetic heterogeneity and enormous variation in exposure levels to environmental agents make it difficult to identify lung cancer susceptibility loci in humans, especially for low-penetrance cancer susceptibility genes. Inbred strains of mice exhibit substantial difference in their susceptibilities to both spontaneous and carcinogen-induced lung adenoma development, thus offering a powerful model system to study human lung cancer susceptibility (6). To date, dozens of genetic loci have been linked to mouse lung adenoma susceptibility through linkage analyses (7). The majority of these genetic loci, or quantitative trait loci (QTL), exert small dose of functions, and some of them could only be detected in specific mouse lines though advanced statistical methods (8). However, major mouse lung adenoma susceptibility loci were also detected. For instance, pulmonary adenoma susceptibility (Pas) 1 was mapped to distal chromosome 6 and accounts for
40% of variances in mouse lung adenoma susceptibility (9). Other QTL, such as Pas2 to Pas4 and pulmonary adenoma resistance (Par) 1 to Par4, are less potent than Pas1 but also make significant contribution to mouse lung adenoma susceptibility in individual experimental systems (10). Specifically, Par1 functions as a major modifier of Pas1 and suppresses mouse lung adenoma development (11, 12).
Par1 was previously mapped to mouse chromosome 11 in two different mouse populations (A/J x Mus spretus) x C57BL/6J and reciprocal (SMXA24 x A/J) x A/J (11, 12), but the derived QTL regions from both studies were broad and the precise location of Par1 has not been delineated. In the present study, we generated Par1 congenic strains of mice from Par1/+ M. spretus and Par1/ A/J mice by targeting a 29-cM chromosomal region between microsatellite markers D11Mit15 and D11Mit301. Subsequent analyses on these congenic mice enabled us to narrow down the Par1 QTL into a 2.0-cM (
2.2 Mb) chromosomal region between Hoxb9 and D11Mit70. Further analyses ensued for genes located in this region. Consequently, multiple genes were identified with nonsynonymous single nucleotide polymorphisms (SNP) or differential lung tissue expression, and tumor suppressor gene Tob1 was studied in particular.
 |
Materials and Methods
|
|---|
Mouse breeding. Inbred A/J and M. spretus mice were purchased from The Jackson Laboratory (Bar Harbor, ME). Animals were housed in plastic cages with hardwood bedding and dust covers, in a HEPA filtered, environmentally controlled room (24 ± 1°C, 12/12-h light/dark cycle). Animals were given Rodent Lab Chow, #5001 (Purina, St. Louis, MO) and water ad libitum. The breeding scheme was to transfer a
29-cM genomic fragment of lung tumorresistant strain M. spretus, encompassed by microsatellite markers D11Mit15 and D11Mit301, onto the genetic background of lung adenoma susceptible strain A/J through sequential multigenerations of backcrossing. Selection of this chromosome 11 region governed the breeding procedure until N7. At N8, 50 male congenic strains containing different chromosomal segments of Par1 QTL were generated. These individual strains were then each crossed with four A/J females to produce the N9 generation. After genotyping, N9 congenic mice with same genotypes of their parental N8 strain were selected for lung tumor bioassay.
Genotyping and phenotyping. Genomic DNA was isolated from mouse tail. Briefly, tail clips were homogenized and incubated overnight at 55°C in nucleic acid lysis solution (Pronase 0.4 mg/mL, 10% sodium dodecylsulfate (w/v), 10 mmol/L Tris, 400 mmol/L NaCl, and 2 mmol/L EDTA) followed by saturated NaCl extraction, precipitation with ice-cold 100% alcohol, washing twice with 70% alcohol, and dissolving in 1x TE buffer. The following markers have been purchased from Research Genetics, Inc. (Huntsville, AL) and used for congenic mouse genotyping: D11Mit15, D11Mit117, D11Mit212, D11Mit70, D11Mit14, D11Mit223, and D11Mit301. Hoxb9 marker was based on a SNP identified in the 3'-untranslated region. The forward primer was end labeled with [
-32P]ATP, and 30 cycles of PCR were done at 94°C for denaturation, 55°C for annealing, and 72°C for extension. Denaturing polyacrylamide gels (8%) were used for resolution of the radiolabeled PCR products followed by autoradiography.
For phenotype collection, 5-week-old N9 mice were given a single i.p. injection of urethane (1 mg/g body weight) in 0.2 mL PBS. All animals were euthanized by CO2 asphyxiation 20 weeks after urethane initiation. A portion of lung tumors and normal tissue were removed and flash frozen in liquid nitrogen. The rest of lungs were fixed in Tellyesniczky's solution and examined with the aid of a dissecting microscope for tumor multiplicity.
Data mining. Position, name, and description of each gene within the Par1 candidate region were retrieved from the latest version of Ensembl mouse genome database (version 39.36, released on June 2006).1 Information of gene expression in mouse lung tissue was retrieved from EST Profile Viewer of National Center for Biotechnology Information (NCBI) Unigene database.2
Reverse transcription-PCR expression and nucleotide polymorphism analyses. Four 6-week-old A/J and M. spretus female mice (two for each strain) were euthanized and their lung tissues were collected for RNA isolation. Briefly, 100 mg lung tissue was pulverized and total RNA was extracted using Trizol reagent according to the manufacturer's protocol (Life Technologies, Gaithersburg, MD). The quality of the isolated RNA was assessed by absorbance at 260 nm, the A260/A280 ratio (1.71.9), and electrophoresis on 1% agarose/formaldehyde gels that indicated the intensity and integrity of the 28S and 18S bands. Two micrograms of total RNA were used in reverse transcription reaction to synthesize the first-strand cDNA using oligo(dT) primer. To improve reverse transcription-PCR (RT-PCR) accuracy, we have done assays at least twice for each gene. We also used DNase Itreated RNA to prepare cDNA samples and set up a genomic DNA control for possible pseudogene amplification. RT-PCR primers were designed for each gene based on their published sequences and will be available on request. For semiquantitative RT-PCR assays, the linear amplification region for each gene was predetermined by plot experiments. The reaction profile generally consisted of one cycle at 95°C for 2 min followed by 30 to 35 cycles at 94°C, 55°C, and 72°C each for 45 s with minor modifications. Twenty-two cycles of amplification was used for the ß-actin control. For the Tob1 gene, a 5% DMSO was used in both PCR and sequencing reaction. Electrophoreses on 1.5% agarose gels were done to resolve PCR products. Amplified PCR products were purified with QIAquick gel extraction kit (Qiagen, Valencia, CA) and subjected to direct sequencing.
Cells and cloning. Complete Tob1 open reading frames (ORF) of A/J and M. spretus mice were amplified from respective strain cDNAs by PCR with primers designed based on published 129/SvJ strain sequence (Genbank accession no. NM_009427). Sequences were confirmed in both directions. The cDNAs were cloned into vector pcDNA3.1, which carries epitopes V5 and His (Invitrogen, Carlsbad, CA). The constructs were sequenced in both directions. Protein expressions were confirmed by transient transfection in 293 cells and Western blotting with V5 antibody (Invitrogen).
NIH3T3 cells (CRL1658; American Type Culture Collection, Manassas, VA) were grown in DMEM with 10% fetal bovine serum. The stable transfected cells were obtained through selection with G418 (Life Technologies, Invitrogen). The expression of Tob1 was examined by Western blotting with V5 antibody.
Cell growth analysis. NIH3T3 cells expressing Tob1 transfectants (0.5 x 104 and 0.5 x 103) were seeded in 24-well plates for the growth curve assay and 60-mm dishes for colony formation assay, respectively. For the growth curve assay, cells were counted from the 2nd day (day 1) after splitting until day 5. The cells were harvested and fixed with 4% paraformaldehyde followed by staining with 0.1% crystal violet for 15 min. The cells were then washed with PBS twice, and 0.2% of Triton X-100 were added to make cells permeable. The absorption value was measured by using 96-well plate reader at 595 nm wavelength. For the colony formation assay, after 10 days, the cells were fixed with ethanol and stained with 0.1% crystal violet. The number of colonies was counted under microscopy.
Western blot. Mouse lungs from A/J and M. spretus were homogenized in lysis buffer [PBS containing 5 mmol/L MgCl2, 1 mmol/L DTT, 0.2 mmol/L phenylmethylsulfonyl fluoride, 2 µg/mL leupeptin, 5 µg/mL aprotinin, 1 mmol/L benzamidine (pH7.5)] and cleared by centrifugation. NP40 was then added to 1% after centrifugation. Protein concentrations were determined by Bradford assay (Bio-Rad, Hercules, CA). Western blotting followed regular procedures. Phosphorylation of Tob1 was detected using 1:100 antihuman phospho-Tob1 antibody (code no. 18911, IBL Co. Ltd., Takasaki-Shi, Gunma, Japan). This antibody has shown crossed-reaction with mouse phospho-Tob1 protein and thus was used. Transient transfections of NIH3T3 cells with V5-tagged Tob1 plasmids were conducted using the manufacturer's provided protocol (Invitrogen).
Immunohistochemistry. Tissue specimens were collected from the lung tumor bioassay. Immunohistochemical analysis with antiphospho-Tob antibody was done by the standard method. In brief, deparaffinized 5-µm thin tissue sections were placed in microwave oven in 0.01 mol/L sodium citrate buffer (pH 6.0), for two periods of 5 min each. After tissue sections were blocked with 0.3% hydrogen peroxide in PBS for 30 min at room temperature, each section was incubated overnight at 4°C with antiphospho-Tob antibody at 1:50. The antirabbit IgG secondary antibody was added next morning. Counterstaining was done with hematoxylin.
Statistical analysis. Two-tailed Student's t test was used to determine P values for the lung tumor multiplicity difference between the congenic strains and control strains. The same test has also been used to compare cell growth-suppressing function among the A/J and M. spretus Tob1 isoforms.
 |
Results
|
|---|
Fine map Par1 QTL into a 2.2-Mb candidate region through congenic mice. Par1 was detected previously on mouse chromosome 11 by two groups using conventional linkage analyses. Manenti et al. (11) studied (A/J x M. spretus) x C57BL/6J mice and detected a peak logarithm of odds (LOD) score of 5.3 at the Rara gene locus (Fig. 1A
). The genetic difference at this locus explained 23% of the phenotypic variance when only mice carrying the highly penetrant Pas1 allele (i.e., Pas1/+) were analyzed, whereas no QTL could be detected in Pas1/ mice. Pataer et al. (12) later worked on reciprocal (SMXA24 x A/J) x A/J mice and detected a peak LOD score of 4.35 at microsatellite marker D11Mit70 (Fig. 1B). In their system, the genetic difference at the Par1 locus explained
15% of total phenotype variance. Both analyses delineated that Par1 acts to negatively modify the Pas1 locus. However, neither study was able to provide a defined candidate region for the Par1 QTL because of limited marker density and population size. We also noticed that based on the present mouse genome map, marker order in the peak region was mistakenly assigned by the linkage software in the Manenti et al. study (Fig. 1A). A correctly defined Par1 candidate region is imperative for positional cloning of the gene. In our previous studies, we have successfully used congenic mice to fine map two major mouse lung adenoma susceptibility loci Pas1 and Par2. Based on the initial linkage mapping results, we selected a 29-cM chromosomal region encompassed by D11Mit15 and D11Mit301 markers to generate Par1 congenic mice. Par1 congenic strains were constructed by transferring this 29-cM fragment from the lung tumorresistant strain M. spretus (Par1/+, donor strain) onto the genetic background of lung tumorsusceptible strain A/J (Par1/, recipient strain) through a total of nine generations of backcrossing. Theoretically, nine generations of backcross mating will render
99.81% of recipient (i.e., A/J strain) genome in congenic mice (13). During the breeding, we genotyped other mouse lung adenoma susceptibility loci (e.g., Pas1Pas4) to assure a clean A/J background.
N9 generations of mice were carrying variant Par1 subcongenic fragments and subjected to a standard urethane-induced lung tumor bioassay (i.e., single i.p. injection at dose of 1 mg/g body weight). As shown in Fig. 1C, we observed a significantly lower lung tumor multiplicity (18.7 ± 5.0 tumors per mouse) in congenic strain AM (heterozygous A/J x M. spretus alleles for the entire 29-cM Par1 region) than in control strain AA (30.5 ± 4.5 tumors per mouse, homozygous A/J alleles for the region; P < 0.00001). Thus, the substitution of one copy of A/J allele by M. spretus in the Par1 QTL decreased the mouse lung tumor multiplicity by
1.6-fold. Consistent with the previous study, we also observed a
2.4-fold decrease in total tumor volume in the AM mice (data not shown). We analyzed the N9 subcongenic mice on their distinguished Par1 segment genotypes and lung tumor phenotypes. As shown in Fig. 1D, the subcongenic strains 3 and 4 exhibited comparable lung tumor multiplicities as congenic strain AM does and significantly lower lung tumor multiplicities than the control strain AA, defining the right boundary of Par1 candidate region at D11Mit70 and the left boundary of Par1 region at Hoxb9, respectively. The refined Par1 QTL region was 2.0 cM in genetic distance and 2.2 Mb in physical distance.
Gene content in the refined Par1 candidate region. Based on the latest Ensembl mouse genome (version 39.36), there are 49 genes located in the 2.2-Mb candidate region (Table 1
). Thirty-two genes have transcript expression in mouse lungs according to the public expression profile database (EST Profile Viewer of NCBI Unigene). Their expression in lung has been confirmed by our subsequent RT-PCR assays and indicated in Table 1. Except ribosomal protein gene Mrpl27 and several functionally unknown Riken transcript genes (e.g., A430060F13Rik and 9530033F24Rik), these lung organ-related genes are involved in a broad range of cellular biochemical or physiologic properties, such as cell growth and differentiation [e.g., Tob1, nerve growth factor receptor (Ngfr), Igf2bp1, and Phb], embryonic development (e.g., Hoxb13), protein phosphorylation (e.g. Pdk2, Ppp1r9b, and Phospho1), DNA replication and gene expression (e.g., Mycbpap and Myst2), drug resistance (e.g., 3300001p08Rik), neuron firing (e.g., Cacnalg), phototransduction (e.g., Gngt2), energy synthesis (e.g., Atp5g1), and intercellular adherence and cell migration (e.g., Itga3 and Abi3). Interestingly, several genes have shown association with cancer development in previous studies. For instance, deficiency of the Tob1 tumor suppressor gene led to development of multiple types of tumors in mice (14).
Differential transcript expression and nonsynonymous nucleotide polymorphism. We conducted expression analysis using cDNA samples prepared from age-matched A/J and M. spretus lung tissues. Among the 32 genes we investigated, majority of the genes exhibited comparable transcript expression between the two mouse strains, including the bona fide tumor suppressor Tob1 (Fig. 2A
). There were three genes (i.e., Itga3, Spop, and Myst2), however, seeming to be more abundant in A/J lungs than in M. spretus lungs (Fig. 2A). In contrast, the Phospho1 gene transcript was more abundant in M. spretus lungs than in A/J lungs (Fig. 2A). Amplification of the Phb gene, whose human homologue was shown associated with breast cancer susceptibility in previous study (15), generated a similar size of product from genomic DNA control. Searching the NCBI mouse genome database with Phb mRNA sequence (Genbank no. NM_008831) revealed multiple Phb pseudogenes whose genomic sequences are highly similar to the Phb mRNA sequence, prompting us to treat RNA samples with DNase-I enzyme. However, there was no differential expression found for the Phb gene (Fig. 2A).

View larger version (20K):
[in this window]
[in a new window]
[Download PPT slide]
|
Figure 2. Gene expression and SNP analyses. A, gene expression analysis. Individual RNA samples were prepared from age-matched A/J and M. spretus lung tissues and treated by DNase-I for semiquantitative RT-PCR analysis. The Itga3, Spop, and Myst2 genes show more transcript expression in A/J lung than in M. spretus lung. In contrast, the Phospho1 gene transcript was more abundant in M. spretus mouse lung than in the A/J lung. The Tob1 and Phb genes do not show differential expression in mouse lung tissues. The ß-actin gene was used as a control to assure equal amount of cDNA used in the reactions. DNA ladder (100 bp) was loaded before sample lanes. B, Tob1 gene nucleotide polymorphisms. Chromatograms of nonsynonymous nucleotide polymorphisms detected in the two mouse strains. Underlined, codons.
|
|
On the other hand, SNP was frequently detected for the investigated genes (data not shown), which can be attributed to the far phylogenetic distance between M. spretus and A/J mouse strains. We only focused on nonsynonymous SNPs because they are the ones likely to cause alteration of protein function. Five genes (Tob1, Xylt2, Itga3, Ngfr, and Phb) were identified with nonsynonymous SNPs (Table 2
). Despite encoding the same number of amino acids, the A/J and M. spretus Tob1 genes show three nonsynonymous SNPs. All of these nonsynonymous polymorphisms have been identified at its COOH terminus, making Tob1 alleles of A/J and M. spretus distinct from each another (Fig. 2B). At codon 246, the A/J Tob1 encode a proline (CCA) compared with a glutamine (CAG) in the M. spretus. In addition, codons 270 and 271 encode Pro-Gln (CCA-CAG) in the A/J Tob1, which is different from Gln-Pro (CAG-CCA) in the M. spretus Tob1. Furthermore, codon 301 encodes a glycine (GGT) in A/J but an aspartic acid (GAT) in M. spretus. The Xylt2 gene has two nonsynonymous SNPs identified at codons 370 and 544. At codon 370, the A/J allele encodes an arginine (CGT), whereas M. spretus allele encodes a leucine (CTT). At codon 544, the A/J allele encodes a threonine (ACC), whereas the M. spretus allele encodes a serine (AGC; Table 2). The Itga3 gene has four nonsynonymous SNPs identified (Table 2). Codon 511 (CAG/CGG: Gln/Arg), codon 537 (TTG/TCA: Leu/Ser), and codon 895 (GTC/GCC: Val/Ala) polymorphisms were identified between M. spretus and A/J (Table 2). In addition, at codon 885, a glutamine (CAG) is encoded by A/J, whereas a histidine (CAT) is encoded by M. spretus. The trinucleotide CTG polymorphism in the Ngfr gene is located in exon 1 and leads to a 6-leucine chain in the Ngfr proteins of M. spretus but a 7-leucine chain in A/J. At codon 49, the A/J Ngfr encodes a proline (CCT), whereas M. spretus Ngfr encodes a leucine (CTT). The codon 217 encodes an isoleucine (ATA) in A/J but a valine (GTA) in M. spretus (Table 2). Phb has only one nonsynonymous SNP identified at codon 66, with an isoleucine (ATC) in A/J and a leucine (CTC) in M. spretus (Table 2).
Evaluation of Tob1 tumor suppressor as Par1 candidate. In the present study, Tob1 was evaluated in particular because of its established role as a tumor suppressor and the multiple nonsynonymous SNPs identified between A/J and M. spretus. It was shown previously that exogenously expressed human TOB1 was able to suppress growth of NIH3T3 cells, and the TOB1 growth suppression pathway is negatively regulated by the ERBB2 gene product p185 (16). Researchers observed multiple types of tumors in Tob1-difficient mice, including lung tumors (14). We wanted to know if the identified nucleotide variations would cause any change in Tob1 antitumor function. We respectively cloned A/J and M. spretus Tob1 ORFs and did functional tests using in vitro cell culture systems. Normal but immortalized NIH3T3 cell line was used in colony formation and cell growth curve assays (Fig. 3A and B
). The results from these two assays suggest that the M. spretus Tob1 allele exerts stronger growth-suppressing effects than the A/J counterpart. Expression levels of the two Tob1 transfectants were confirmed to be comparable (data not shown). These results suggest that the nonsynonymous SNPs identified between the A/J and M. spretus Tob1 genes might indeed lead to some functional differences.

View larger version (24K):
[in this window]
[in a new window]
[Download PPT slide]
|
Figure 3. In vitro cell culture studies for the three mouse Tob1 isoforms. A, colony formation assay. Left, the control (pcDNA3.1) and Tob1-transfected NIH3T3 cell culture dishes stained with 0.1% crystal violet after 10 days of growth; right, the number of colonies was counted. Student's t test was used to compare the cell growth-suppressing effects of M. spretus Tob1 with A/J Tob1. **, P < 0.01. Bars, SD. B, cell growth curve assay. The cells were harvested and fixed with 4% paraformaldehyde, stained with 0.1% crystal violet for 15 min, then washed with PBS twice, and treated by 0.2% Triton X-100. The absorption value was measured by using 96-well plate reader at 595 nmol/L wavelength. Cells were counted from the 2nd day after splitting (day 1) until day 5. Consistent with the results from colony formation assay, M. spretus Tob1transfected cells have grown more slowly than the A/J Tob1transfected cells. Bars, SD. C, phosphorylation ability of the A/J and M. spretus Tob1 isoforms. Transient transfection of V5-tagged Tob1 plasmids into NIH3T3 cells did not show significant difference in phosphorylation levels of ectopically expressed Tob1 proteins.
|
|
Because antitumor function of Tob1 is largely regulated by protein phosphorylation (17), we wondered if the three mouse Tob1 isoforms were different in their abilities to get phosphorylated by cellular kinases, which might explain the differential cell growth-suppressing function observed in the cell culture studies. Indeed, one recent study showed that alteration of phosphorylation status of Tob1 is likely to contribute to early human lung tumorigenesis (18). More interestingly, in the present study, we observed positive immunohistochemical staining of phosphorylated Tob1 protein in mouse lung adenomas but not in normal mouse lung tissues (Fig. 4A and B
), suggesting that phosphorylation of Tob1 might play an important role in mouse lung adenoma development. However, transient transfection of V5-tagged Tob1 into NIH3T3 cells did not show conspicuous difference in phosphorylation levels of ectopically expressed Tob1 proteins (Fig. 3C). The result is consistent with the fact that the major phosphorylation sites of Tob1 protein are located at the NH2 terminus and remain identical for the three mouse strains (17). Therefore, the underlying mechanism for the observed functional difference between Tob1 isoforms remains to be identified.

View larger version (102K):
[in this window]
[in a new window]
[Download PPT slide]
|
Figure 4. Immunohistochemistry staining of phosphorylated Tob1 protein. Magnifications, x100 and x400. Arrows, positive staining of phosphorylated mouse Tob1. Positive immunohistochemical staining of phosphorylated Tob1 protein was readily observed in mouse lung adenomas (A) but almost undetectable in normal lung tissues (B).
|
|
 |
Discussion
|
|---|
In our previous work, we conducted positional cloning for the major mouse lung adenoma susceptibility QTL Pas1 (19). Par1 has been shown to function as a major modifier of Pas1 to reduce mouse lung adenoma susceptibility (11, 12). Conventional linkage analyses previously mapped the Par1 QTL into a rather large and undefined candidate region on mouse chromosome 11 (11, 12). To identify the underlying gene, refinement of the Par1 candidate region is a prerequisite step. In the present study, we fine mapped the Par1 QTL using a series of congenic and subcongenic mice generated from the Par1/+ strain M. spretus and Par1/ strain A/J. Analyzing these mice confirmed that the Par1 is a major mouse lung adenoma resistance QTL. Most importantly, our analyses have in the first time generated a defined small Par1 candidate region, which is
2.2 Mb in physical length and is defined by the Hoxb9 gene and D11Mit70.
There are 49 genes located in the refined Par1 candidate region and involved in a broad range of cellular activities. Thirty-two genes have shown different extents of expression in mouse lungs (Table 1). Semiquantitative RT-PCR assays revealed three genes (Itga3, Spop, and Myst2) more abundantly expressed in A/J versus M. spretus lung tissues. A decreased transcription level was observed in A/J versus M. spretus for the Phospho1 gene. Nonsynonymous SNPs were identified in five genes (Tob1, Xylt2, Itga3, Ngfr, and Phb; Table 2). The human SPOP gene encodes a 374amino acid protein with a relative molecular mass of 47,000 (20). Its amino acid sequence contains a COOH-terminal BTB/POZ domain and an NH2-terminal MATH domain (2022). It has been shown that the TRAF homology domain in SPOP protein only interacted weakly with TRAF1 and TRAF6 and did not inhibit nuclear factor-
B (NF-
B) induction as other TRAF domains do (21), and the MATH domain has been shown to bind to the X chromosome inactivation protein macroH2A1.2 (22). The Myst2 gene (also named Hbo1) is a histone acetyltransferase, which interacted with the origin recognition complex 1 protein and the minichromosome maintenance protein, suggesting a role for it in control of both DNA replication and gene expression (23, 24). Another study also showed that the Myst2 protein interacted with androgen receptor (AR) and might regulate AR-dependent genes in normal and prostate cancer cells (25). The Phospho1 gene was predicted as a putative phosphatase, but thus far, no study has been done on its function. Pdk2 presumably inactivates pyruvate dehydrogenase through phosphorylation, therefore regulating oxidative decarboxylation of pyruvate and homeostasis of carbohydrate fuels in mammals (26). The biological function of Xylt2 is not known. It is highly homologous to Xylt1, which catalyzes the transfer of UDP-xylose to serine residues within XT recognition sequences of target proteins (27). Addition of this xylose to the core protein is required for the biosynthesis of the glycosaminoglycan chains characteristic of proteoglycans. The human PHB is a 30-kDa intracellular, antiproliferative protein and has shown association with breast cancer susceptibility (15).
The human NGFR gene encodes a 75-kDa glycoprotein cell surface receptor (p75NTR), which binds with low affinity to nerve growth factor and other neurotrophins (28). The p75NTR protein shares sequence identity with other members in a receptor superfamily (such as tumor necrosis factor receptors p75TNFR, p55TNFR, and FAS) and is involved in NF-
B and apoptosis pathways. It should be noted that although NGFR is expressed in many organs outside the nervous system, recent immunohistochemistry staining studies in human lung have detected the p75NTR protein only in lung ganglionic neuron and artery smooth muscle but not in any type of pneumocytes (29). Indeed, using highly sensitive RT-PCR technology, we only detected weak transcript expression (data not shown), which is likely derived from these nonpulmonary tissues. On the other hand, both differential lung expression and nonsynonymous nucleotide polymorphisms were found for the Itga3 gene. The human ITGA3 gene encodes two isoforms of the integrin
3 subunit:
3A (1,051 amino acids) and
3B (1,066 amino acids). The mouse Itga3 gene encodes a 1,053amino acid protein, which corresponds to the human
3A isoform. Previous studies have revealed that only the
3A isoform is expressed in lung tissue (30). Functionally, the integrin
3 subunit forms a heterodimeric cell surface integrin receptor with the integrin ß1 subunit, linking the extracellular matrix to structural and functional components within the cell (31). The role of Itga3 in tumorigenesis seems to be controversial and may be dependent on tumor type and tumor stage. Some studies have revealed that reduced expression of Itga3 is associated with malignant transformation and Itga3 gene expression may promote apoptosis (3234). Clinically, it has been reported that the reduced Itga3 expression is associated with poor prognosis of patients with lung adenocarcinomas and colon cancer (35, 36). Down-regulation of the Itga3 may also contribute to the enhanced tumorigenicity of c-mycoverexpressing small cell lung cancer cells (37).
In the present study, we were particularly interested in the Tob1 gene, a tumor suppressor whose loss of function is involved in tumorigenesis (14). We found that the A/J and M. spretus encode different isoforms of the Tob1 gene. As revealed by the colony formation and growth curve studies, stronger cell growth-suppressing effects were observed for M. spretus Tob1 than for the A/J Tob1. Previous studies have shown that the TOB1 growth suppression pathway is negatively regulated by the kinase-active ERBB2 gene product p185 (16). Although the molecular mechanism by which TOB1 suppresses cell growth remains to be further clarified, phosphorylation has been suggested to be a crucial mechanism to regulate Tob1 functions. Phosphorylation of three regulatory serines (Ser152, Ser154, and Ser164) of Tob1 by extracellular signal-regulated kinase (Erk) 1 and Erk2 is required for Ras-mediated cell proliferation and transformation (17). The connection of Tob1 gene with the Ras signaling pathway is of large interest because the K-ras gene has been identified as one of the Pas1 candidates in our recent study (19). In the present study, we indeed observed a significantly elevated phosphorylation of Erk protein in A/J lung tissue compared with M. spretus.3 Because K-ras is one of critical upstream regulators of Erk activity, elevated phosphorylation of Erk may reflect a higher K-ras activity in A/J (Pas1/+) than in M. spretus (Pas1/), as we have seen for A/J and C57BL/6 (Pas1/) mice in the previous study (19). However, because we could not detect differential phosphorylation of ectopically expressed Tob1 protein, connection of Tob1 phosphorylation in A/J and M. spretus lungs with their differential Erk activities remains to be a question. Our identified amino acidchanging nucleotide polymorphisms are all located in the COOH-terminal region. It has been suggested that the COOH-terminal region of TOB1 controls protein stability (38). Furthermore, the identified nucleotide polymorphisms may also be able to affect intracellular distribution of Tob1 proteins, which was suggested to be important for the tumor-suppressing function of the gene (39). Clearly, more in vitro and in vivo studies will help to confirm Tob1 gene candidacy, although indisputable evidence for claiming Tob1 as the Par1 gene may not be available until Tob1 knockin mice are generated.
Human chromosomal region 17q21 shares a high degree of homology with the mouse Par1 region. Loss of heterozygosity (LOH) in the human 17q21 region has been observed in multiple types of human cancers, including lung cancer (Fig. 5
; refs. 4046). Most LOHs in cancers of organs other than lung were observed in the region near BRCA1 locus (38.4 Mb, NCBI human genome map, Build 36.1). However, in lung cancer, the "hottest" LOH spot seems not to be at BRCA1. It was previously found in Japanese lung cancer patients that the highest LOH frequency was found at D17S588 (45.6 Mb; ref. 46). In our recent LOH study in Caucasian lung cancer patients, we observed highest LOH frequency at D17S1869 (38%).3 Interestingly, this marker is located at 45.9-Mb chromosomal position and is proximate to the TOB1 gene (46.3 Mb), suggesting that TOB1 gene may function as a tumor suppressor in human lung cancer development.

View larger version (9K):
[in this window]
[in a new window]
[Download PPT slide]
|
Figure 5. LOH of chromosome 17q21 in human lung cancer. Percentage of LOH for each type of cancer was highest frequency observed and was derived from the reference articles. Most LOH in cancers of organs other than lung were observed in the region near BRCA1 locus (38.4 Mb, NCBI human genome map, Build 36.1). However, in lung cancer, the "hottest" LOH spot seems to be closer to Tob1 locus (46.3 Mb).
|
|
It is noteworthy that in the present study, we only studied gene expression and SNPs between A/J and M. spretus. The reason we did not test SM/J mice is because we could not see significant phenotype difference in congenic mice generated between A/J and SM/J strains.3 Therefore, it is likely that the previously identified Par1 QTL in the SMXA24 RI strain might be derived from some fixed mutations emerging in the long term of breeding process. On the other hand, although the Par1 QTL were identified in both M. spretus and SM/J at similar chromosomal positions through linkage analyses (11, 12), it may not necessarily represent the same gene. Because of the far phylogenic distance between M. spretus and other inbred mouse strains, it is perceivable that a unique Par1 gene may function in M. spretus. In future, genetic analysis on mouse populations generated from M. spretus and SM/J may provide us a conclusive answer as to whether these two strains share the same Par1 gene.
 |
Acknowledgments
|
|---|
Grant support: NIH grants CA099187, CA099147, ES012063, and ES013340.
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
 |
Footnotes
|
|---|
Note: M. Wang and Z. Zhang contributed equally to this work.
1 http://www.ensembl.org/Mus_musculus/index.html 
2 http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=unigene 
3 Unpublished data. 
Received 8/24/06.
Revised 11/ 9/06.
Accepted 1/10/07.
 |
References
|
|---|
- American Cancer Society. Cancer Statistics 2005.
- Sellers TA, Bailey-Wilson JE, Elston RC, et al. Evidence for mendelian inheritance in the pathogenesis of lung cancer. J Natl Cancer Inst 1990;82:12729.[Abstract/Free Full Text]
- Sellers TA, Potter JD, Folsom AR. Association of incident lung cancer with family history of female reproductive cancers: the Iowa Women's Health Study. Genet Epidemiol 1991;8:199208.[CrossRef][Medline]
- Sellers TA, Elston RC, Atwood LD, Rothschild H. Lung cancer histologic type and family history of cancer. Cancer (Phila) 1992;69:8691.
- Rudier HW. Genetic predisposition to the toxic effects of chemicals. In: Grandjean P, editor. Ecogenetics. New York: Chapman & Hall; 1991. p. 20516.
- Malkinson AM. The genetic basis of susceptibility to lung tumors in mice. Toxicology 1989;54:24171.[CrossRef][Medline]
- Demant P. Cancer susceptibility in the mouse: genetics, biology, and implications for human cancer. Nat Rev Genet 2003;4:72134.[Medline]
- Fijneman RJ, de Vries SS, Jansen RC, Demant P. Complex interactions of new quantitative trait loci, Sluc1, Sluc2, Sluc3, and Sluc4, that influence the susceptibility to lung cancer in the mouse. Nat Genet 1996;14:4657.[CrossRef][Medline]
- Gariboldi M, Manenti G, Canzian F, et al. A major susceptibility locus to murine lung carcinogenesis maps on chromosome 6. Nat Genet 1993;3:1326.[CrossRef][Medline]
- Herzog CR, Lubet RA, You M. Genetic alterations in mouse lung tumors: implications for cancer chemoprevention. J Cell Biochem Suppl 1997;28-29:4963.
- Manenti G, Gariboldi M, Elango R, et al. Genetic mapping of a pulmonary adenoma resistance (Par1) in mouse. Nat Genet 1996;12:4557.[CrossRef][Medline]
- Pataer A, Nishimura M, Kamoto T, Ichioka K, Sato M, Hiai H. Genetic resistance to urethan-induced pulmonary adenomas in SMXA recombinant inbred mouse strains. Cancer Res 1997;57:29048.[Abstract/Free Full Text]
- Markel P. Theoretical and empirical issues for marker-assisted breeding of congenic mouse strains. Nat Genet 1997;17:2804.[CrossRef][Medline]
- Yoshida Y, Nakamura T, Komoda M, et al. Mice lacking a transcriptional corepressor Tob are predisposed to cancer. Genes Dev 2003;7:12016.
- Sato T, Saito H, Swensen J, et al. The human prohibitin gene located on chromosome 17q21 is mutated in sporadic breast cancer. Cancer Res 1992;52:16436.[Abstract/Free Full Text]
- Matsuda S, Kawamura-Tsuzuku J, Ohsugi M, et al. Tob, a novel protein that interacts with p185erbB2, is associated with anti-proliferative activity. Oncogene 1996;12:7053.[Medline]
- Suzuki T, K-Tsuzuku J, Ajima R, Nakamura T, Yoshida Y, Yamamoto T. Phosphorylation of three regulatory serines of Tob by Erk1 and Erk2 is required for Ras-mediated cell proliferation and transformation. Genes Dev 2002;16:135670.[Abstract/Free Full Text]
- Iwanaga K, Sueoka N, Sato A, et al. Alteration of expression or phosphorylation status of tob, a novel tumor suppressor gene product, is an early event in lung cancer. Cancer Lett 2003;202:719.[CrossRef][Medline]
- Zhang Z, Futamura M, Vikis HG, et al. Positional cloning of the major quantitative trait locus underlying lung tumor susceptibility in mice. Proc Natl Acad Sci U S A 2003;100:126427.[Abstract/Free Full Text]
- Nagai Y, Kojima T, Muro Y, et al. Identification of a novel nuclear speckle-type protein, SPOP. FEBS Lett 1997;418:236.[CrossRef][Medline]
- Zapata JM, Pawlowski K, Haas E, Ware CF, Godzik A, Reed JC. A diverse family of proteins containing tumor necrosis factor receptor-associated factor domains. J Biol Chem 2001;276:2424252.[Abstract/Free Full Text]
- Takahashi I, Kameoka Y, Hashimoto K. MacroH2A1.2 binds the nuclear protein Spop. Biochim Biophys Acta 2002;1591:638.[Medline]
- Iizuka M, Stillman B. Histone acetyltransferase HBO1 interacts with the ORC1 subunit of the human initiator protein. J Biol Chem 1999;274:2302734.[Abstract/Free Full Text]
- Burke TW, Cook JG, Asano M, Nevins JR. Replication factors MCM2 and ORC1 interact with the histone acetyltransferase HBO1. J Biol Chem 2001;276:15397408.[Abstract/Free Full Text]
- Sharma M, Zarnegar M, Li X, Lim B, Sun Z. Androgen receptor interacts with a novel MYST protein, HBO1. J Biol Chem 2000;275:352008.[Abstract/Free Full Text]
- Korotchkina LG, Sidhu S, Patel MS. Characterization of testis-specific isoenzyme of human pyruvate dehydrogenase. J Biol Chem 2006;281:968896.[Abstract/Free Full Text]
- Gotting C, Kuhn J, Brinkmann T, Kleesiek K. Xylosylation of alternatively spliced isoforms of Alzheimer APP by xylosyltransferase. J Prot Chem 1998;17:295302.[CrossRef][Medline]
- Johnson D, Lanahan A, Buck CR, et al. Expression and structure of the human NGF receptor. Cell 1986;47:54554.[CrossRef][Medline]
- Ricci A, Felici L, Mariotta S, et al. Neurotrophin and neurotrophin receptor protein expression in the human lung. Am J Respir Cell Mol Biol 2004;1:129.
- Tamura RN, Cooper HM, Collo G, Quaranta V. Cell type-specific integrin variants with alternative
chain cytoplasmic domains. Proc Natl Acad Sci U S A 2001;88:101837. - Hynes RO. Integrins: versatility, modulation, and signaling in cell adhesion. Cell 1992;69:1125.[CrossRef][Medline]
- Valea FA, Haskill S, Moore DH, Fowler WC, Jr. Immunohistochemical analysis of
1-integrins in cervical cancer. Am J Obstet Gynecol 1995;173:80813.[CrossRef][Medline] - Weitzman JB, Hemler ME, Brodt P. Reduction of tumorigenicity by
3 integrin in a rhabdomyosarcoma cell line. Cell Commun Adhes 1996;4:4152. - Seewaldt VL, Mrozek K, Sigle R, et al. Suppression of p53 function in normal human mammary epithelial cells increases sensitivity to extracellular matrix-induced apoptosis. J Cell Biol 2001;155:47186.[Abstract/Free Full Text]
- Adachi M, Taki T, Huang C, et al. Reduced integrin
3 expression as a factor of poor prognosis of patients with adenocarcinoma of the lung. J Clin Oncol 1998;16:10607.[Abstract] - Hashida H, Takabayashi A, Tokuhara T, et al. Integrin
3 expression as a prognostic factor in colon cancer: association with MRP-1/CD9 and KAI1/CD82. Int J Cancer 2002;97:51825.[CrossRef][Medline] - Barr LF, Campbell SE, Bochner BS, Dang CV. Association of the decreased expression of
3ß1 integrin with the altered cell: environmental interactions and enhanced soft agar cloning ability of c-myc-overexpressing small cell lung cancer cells. Cancer Res 1998;58:553745.[Abstract/Free Full Text] - Sasajima H, Nakagawa K, Yokosawa H. Antiproliferative proteins of the BTG/Tob family are degraded by the ubiquitin-proteasome system. Eur J Biochem 2002;269:3596604.[Medline]
- Kawamura-Tsuzuku J, Suzuki T, Yoshida Y, Yamamoto T. Nuclear localization of Tob is important for regulation of its antiproliferative activity. Oncogene 2004;23:66308.[CrossRef][Medline]
- Garcia-Patino E, Gomendio B, Lleonart M, et al. Loss of heterozygosity in the region including the BRCA1 gene on 17q in colon cancer. Cancer Genet Cytogenet 1998;104:11923.[CrossRef][Medline]
- Cropp CS, Nevanlinna HA, Pyrhonen S, et al. Evidence for involvement of BRCA1 in sporadic breast carcinomas. Cancer Res 1994;54:254851.[Abstract/Free Full Text]
- Acevedo CM, Henriquez M, Emmert-Buck MR, Chuaqui RF. Loss of heterozygosity on chromosome arms 3p and 6q in microdissected adenocarcinomas of the uterine cervix and adenocarcinoma in situ. Cancer 2002;94:793802.[CrossRef][Medline]
- Sirchia SM, Sironi E, Grati FR, et al. Losses of heterozygosity in endometrial adenocarcinomas: positive correlations with histopathological parameters. Cancer Genet Cytogenet 2000;21:15662.
- Huang LW, Garrett AP, Schorge JO, et al. Distinct allelic loss patterns in papillary serous carcinoma of the peritoneum. Am J Clin Pathol 2000;114:939.[Abstract/Free Full Text]
- Shiina H, Breault JE, Basset WW, et al. Functional loss of the
-catenin gene through epigenetic and genetic pathways in human prostate cancer. Cancer Res 2005;65:21308.[Abstract/Free Full Text] - Abujiang P, Mori TJ, Takahashi T, et al. Loss of heterozygosity (LOH) at 17q and 14q in human lung cancers. Oncogene 1998;7:302933.