Use of intercross outbred mice and single nucleotide polymorphisms to map skin cancer modifier loci

Car-R and Car-S outbred mouse lines, phenotypically selected for resistance and susceptibility to skin carcinogenesis respectively, show significant linkage disequilibrium (LD) at genetic markers mapping on chromosomal regions where skin cancer modifier loci (Skts3, Skts1, and Psl1 on Chrs 5, 7, and 9 respectively) have been mapped in standard crosses. Analysis of these regions for genetic linkage with skin cancer phenotypes in 245 (Car-R × Car-S)F₂ intercross mice, by using single nucleotide polymorphisms (SNPs), revealed significant linkage at a possible allelic form of the Skts1 locus, whose mapping region was shortened to a <5.5-cM interval near the Tyr locus. The Car-derived Skts1 locus was linked with papilloma multiplicity and latency by a recessive inheritance of the susceptibility allele. Putative loci on Chr 5 (Skts3) and 9 (Psl1) showed no significant linkage. These results point to the important role of the Stks1 locus in mouse skin tumorigenesis in independent crosses. The shortened Skts1 mapping region should facilitate the identification of candidate genes. Received: 23 June 2000 / Accepted: 22 November 2000     

The Impact of 3′UTR Variants on Differential Expression of Candidate Cancer Susceptibility Genes

Variants in regulatory regions are predicted to play an important role in disease susceptibility of common diseases. Polymorphisms mapping to microRNA (miRNA) binding sites have been shown to disrupt the ability of miRNAs to target genes resulting in differential mRNA and protein expression. Skin tumor susceptibility 5 (Skts5) was identified as a locus conferring susceptibility to chemically-induced skin cancer in NIH/Ola by SPRET/Outbred F1 backcrosses. To determine if polymorphisms between the strains which mapped to putative miRNA binding sites in the 3′ untranslated region (3′UTR) of genes at Skts5 influenced expression, we conducted a systematic evaluation of 3′UTRs of candidate genes across this locus. Nine genes had polymorphisms in their 3′UTRs which fit the linkage data and eight of these contained polymorphisms suspected to interfere with or introduce miRNA binding. 3′UTRs of six genes, Bcap29, Dgkb, Hbp1, Pik3cg, Twistnb, and Tspan13 differentially affected luciferase expression, but did not appear to be differentially regulated by the evaluated miRNAs predicted to bind to only one of the two isoforms. 3′UTRs from four additional genes chosen from the locus that fit less stringent criteria were evaluated. Ifrd1 and Etv1 showed differences and contained polymorphisms predicted to disrupt or create miRNA binding sites but showed no difference in regulation by the miRNAs tested. In summary, multiple 3′UTRs with putative functional variants between susceptible and resistant strains of mice influenced differential expression independent of predicted miRNA binding.     

Mapping and identifying genes for asthma and psoriasis

Susceptibility genes for complex diseases are characterized by reduced penetrance, caused by the influence of other genes, the environment or stochastic events. Recently, positional cloning efforts have yielded several candidate susceptibility genes in different complex disorders such as Crohn's disease and asthma. Within a genetic locus, however, the identification of the effector gene may pose further challenges and require functional studies. I review two examples of such challenges: the cloning of GPR154 (GPRA) and AAA1 on chromosome 7p14 at a susceptibility locus for atopy and asthma, and the study of HLA-Cw6, CCHCR1 (HCR) and CDSN on chromosome 6p21 at PSORS1, the major susceptibility locus for psoriasis. The susceptibility locus for atopy and asthma contains two genes and only one of them is protein coding. We studied its isoform-specific expression in bronchial biopsies and in a mouse model of ovalbumin-induced inflammation of bronchial epithelia. In the PSORS1 locus, strong linkage disequilibrium between genes has made it difficult to distinguish the effects of the three nearby genes. We engineered transgenic mice with either a HCR non-risk allele or the HCR*WWCC risk allele controlled by the cytokeratin-14 promoter. The results suggested that the overexpression of HCR in mouse skin was insufficient to induce a psoriasiform phenotype, but it appeared to induce allele-specific gene expression changes that were similar to those observed in psoriatic skin.     

Epistatic interactions between modifier genes confer strain-specific redundancy for Tgfb1 in developmental angiogenesis

Tgfbm1 (chromosome 5, P = 8 x 10(-5)) and Tgfbm3 (chromosome 12, P = 6 x 10(-11)) were identified as loci that modify developmental angiogenesis of Tgfb1 -/- mice. Congenic mice validated these loci and demonstrated epistatic interaction between them. The novel locus, Tgfbm3, encompasses approximately 22 genes, colocalizes with both tumor susceptibility and atherosclerosis susceptibility loci, and is enriched in genes regulating cell growth and morphogenesis. The use of gene knockout and/or transgenic mice that predispose to a complex trait, such as vascular development/angiogenesis, facilitates the identification of modifiers by simplifying genetic analysis. Identification of genes that modify response to lack of transforming growth factor beta1 (TGFbeta1) will enhance the understanding of TGFbeta1 action in vivo and may help predict which patients would respond well to anti-TGFbeta therapy. Identification of angiogenesis-modifying genes may provide new targets for angiogenesis therapies and analysis of polymorphisms therein may contribute to assessment of risk for diseases involving angiogenesis.     

Comparative genomic sequencing reveals a strikingly similar architecture of a conserved syntenic region on human chromosome 11p15.3 (including gene ST5) and mouse chromosome 7

Comparative genomics is a superior way to identify phylogenetically conserved features like genes or regions involved in gene regulation. The comparison of extended orthologous chromosomal regions should also reveal other characteristic traits essential for chromosome or gene function. In the present study we have sequenced and compared a region of conserved synteny from human chromosome 11p15.3 and mouse chromosome 7. In human, this region is known to contain several genes involved in the development of various disorders like Beckwith-Wiedemann overgrowth syndrome and other tumor diseases. Furthermore, in the neighboring chromosome region 11p15.5 extensive imprinting of genes has been reported which might extend to region 11p15.3. The analysis of approximately 730 kb in human and 620 kb in mouse led to the identification of eleven genes. All putative genes found in the mouse DNA were also present in the same order and orientation in the human chromosome. However, in the human DNA one putative gene of unknown function could be identified which is not present in the orthologous position of the mouse chromosome. The sequence similarity between human and mouse is higher in transcribed and exon regions than in non-transcribed segments. Dot plot analysis, however, reveals a surprisingly well-conserved sequence similarity over the entire analyzed region. In particular, the positions of CpG islands, short regions of very high GC content in the 5' region of putative genes, are similar in human and mouse. With respect to base composition, two distinct segments of significantly different GC content exist as well in human as in the mouse. With a GC content of 45% the one segment would correspond to "isochore H1" and the other segment (39% GC in human, 40% GC in mouse) to "isochore L1/L2". The gene density (one gene per 66 kb) is slightly higher than the average calculated for the complete human genome (one gene per 90 kb). The comparison of the number and distribution of repetitive elements shows that the proportion of human DNA made up by interspersed repeats (43.8%) is significantly higher than in the corresponding mouse DNA (30.1%). This partly explains why the human DNA is longer between the landmark genes used to define the orthologous positions in human and mouse.     

Genetic mapping of Prm-1, Igl-1, Smst, Mtv-6, Sod-1, and Ets-2 and localization of the Down syndrome region on mouse chromosome 16

Molecular probes were used as markers in the backcross (Czech II X BALB/cPt) X Czech II to determine the positions of six genes on mouse chromosome 16 (MMU 16). The order of the genes mapped is (centromere), protamine-1 (Prm-1), immunoglobulin lambda 1 light chain (Igl-1), preprosomatostatin (Smst), an endogenous mouse mammary tumor virus locus (Mtv-6), and two more distal sequences, superoxide dismutase, cytoplasmic form (Sod-1), and the proto-oncogene sequence Ets-2. The largest recombination frequency between any two adjacent markers is 24 cM, and thus the position of any marker on MMU 16 that is polymorphic between these two strains can be readily determined in this backcross. A region of MMU 16 which corresponds to the Down syndrome region of human chromosome 21 is located near the distal end of the chromosome.     

Genetic identification of multiple loci that control breast cancer susceptibility in the rat.

We have used a rat model of induced mammary carcinomas in an effort to identify breast cancer susceptibility genes. Using genetic crosses between the carcinoma-resistant Copenhagen (COP) and carcinoma-sensitive Wistar-Furth rats, we have confirmed the identification of the Mcs1 locus that modulates tumor number. We have now also identified two additional loci, Mcs2 and Mcs3. These three loci map to chromosomes 2, 7, and 1, respectively, and interact additively to suppress mammary carcinoma development in the COP strain. They are responsible for a major portion of the tumor-resistant phenotype of the COP rat. No loss of heterozygosity was observed surrounding the three loci. A fourth COP locus, Mcs4, has also been identified on chromosome 8 and acts in contrast to increase the number of carcinomas. These results show that mammary carcinoma susceptibility in the COP rat is a polygenic trait. Interestingly, a polymorphism in the human genomic region homologous to the rat Mcs4 region is associated with an increased breast cancer risk in African-American women. The isolation of the Mcs genes may help elucidate novel mechanisms of carcinogenesis, provide information important for human breast cancer risk estimation, and also provide unique drug discovery targets for breast cancer prevention.     

Genetic mapping of a cloned sequence responsible for susceptibility to ecotropic murine leukemia viruses.

A mouse cDNA that confers susceptibility to ecotropic murine leukemia viruses following transfection into human EJ cells has been cloned and sequenced. We show that this sequence is likely to be Rec-1, the chromosome 5 locus originally defined by studies with somatic cell hybrids as responsible for virus susceptibility, and provide a specific chromosomal map position for this locus by analysis of an interspecies backcross. This locus maps in the distal region of chromosome 5 and is thus not within the cluster of retrovirus-related genes near the centromere.     

Construction of a 1.2-Mb BAC/PAC contig of the porcine gene RYR1 region on SSC 6q1.2 and comparative analysis with HSA 19q13.13

We screened a porcine bacterial artificial chromosome (BAC) and a P1 derived artificial chromosome (PAC) library to construct a sequence-ready approximately 1.2-Mb BAC/PAC contig of the ryanodine receptor-1 gene (RYR1) region on porcine chromosome (SSC) 6q1.2. This genomic segment is of special interest because it harbors the locus for stress susceptibility in pigs and a putative quantitative trait locus for muscle growth. Detailed physical mapping of this gene-rich region allowed us to assign to this contig 17 porcine genes orthologous to known human chromosome 19 genes. Apart from the relatively well-characterized porcine gene RYR1, the other 16 genes represent novel chromosomal assignments and 14 genes have been cloned for the first time in pig. Comparative analysis of the porcine BAC/PAC contig with the human chromosome (HSA) 19q13.13 map revealed a completely conserved gene order of this segment between pig and human. A detailed porcine-human-mouse comparative map of this region was constructed.     

An integrated radiation hybrid map of bovine chromosome 18 that refines a critical region associated with multiple ocular defects in cattle

Congenital multiple ocular defects (MOD) of Japanese black cattle is a hereditary ocular disorder with an autosomal recessive mode of inheritance showing developmental defects of the lens, retina and iris, persistent embryonic eye vascularization and microphthalmia. The MOD locus has been mapped by linkage analysis to a 6.6-cM interval on the proximal end of bovine chromosome 18, which corresponds to human chromosome 16q and mouse chromosome 8. To refine the MOD region in cattle, we constructed an integrated radiation hybrid (RH) map of the proximal region of bovine chromosome 18, which consisted of 17 genes and 10 microsatellite markers, using the SUNbRH7000 panel. Strong conservation of gene order was found among the corresponding chromosomal regions in cattle, human and mouse. The MOD-critical region was fine mapped to a 59.5-cR region that corresponds to a 6.3-Mb segment of human chromosome 16 and a 4.8-Mb segment of mouse chromosome 8. Several positional candidate genes, including FOXC2 and USP10, were identified in this region.     

Overexpression of Hoxc13 in differentiating keratinocytes results in downregulation of a novel hair keratin gene cluster and alopecia

Studying the roles of Hox genes in normal and pathological development of skin and hair requires identification of downstream target genes in genetically defined animal models. We show that transgenic mice overexpressing Hoxc13 in differentiating keratinocytes of hair follicles develop alopecia, accompanied by a progressive pathological skin condition that resembles ichthyosis. Large-scale analysis of differential gene expression in postnatal skin of these mice identified 16 previously unknown and 13 known genes as presumptive Hoxc13 targets. The majority of these targets are downregulated and belong to a subgroup of genes that encode hair-specific keratin-associated proteins (KAPs). Genomic mapping using a mouse hamster radiation hybrid panel showed these genes to reside in a novel KAP gene cluster on mouse chromosome 16 in a region of conserved linkage with human chromosome 21q22.11. Furthermore, data obtained by Hoxc13/lacZ reporter gene analysis in mice that overexpress Hoxc13 suggest negative autoregulatory feedback control of Hoxc13 expression levels, thus providing an entry point for elucidating currently unknown mechanisms that are required for regulating quantitative levels of Hox gene expression. Combined, these results provide a framework for understanding molecular mechanisms of Hoxc13 function in hair growth and development.     

Physical mapping of the pink-eyed dilution complex in mouse chromosome 7 shows that Atp10c is the only transcript between Gabrb3 and Ube3a.

Phenotypic analyses of a set of homozygous-lethal deletion mutants at the pink-eyed dilution (p) locus has resulted in the identification of p-linked obesity locus 1 (plo 1), distal to the p locus, as a locus involved in the modulation of body fat and/or affecting lipid metabolism in these mice. The plo 1 region maps to mouse chromosome 7 (MMU 7) between two genes, Gabrb3 and Ube3a, which have been used as anchor points to generate an integrated deletion and physical map of plo 1 that encompasses about 1.2-1.3 Mb. A deletion/physical map was constructed and the genomic DNA between the two loci was sequenced to identify genes mapping to this region. Data show that Atp10c, a novel type IV ATPase a putative phospholipid transporter, is the only coding unit in this region of the chromosome.     

High-resolution linkage mapping for susceptibility genes in human polygenic disease: Insulin-dependent diabetes mellitus and chromosome 11q

Insulin-dependent diabetes mellitus (IDDM) has a complex pattern of genetic inheritance. In addition to genes mapping to the major histocompatibility complex (MHC), several lines of evidence point to the existence of other genetic susceptibility factors. Recent studies of the nonobese diabetic mouse (NOD) model of IDDM have suggested the presence, on mouse chromosome 9, of a susceptibility gene linked to the locus encoding the T-cell antigen, Thy-1. A region on human chromosome 11q is syntenic to this region on mouse chromosome 9. We have used a set of polymorphic DNA markers from chromosome 11q to investigate this region for linkage to a susceptibility gene in 81 multiplex diabetic pedigrees. The data were investigated by maximization of lod scores over genetic models and by multiple-locus affected-sib-pair analysis. We were able to exclude the presence of a susceptibility gene (location scores less than -2) throughout greater than 90% of the chromosome 11q homology region, under the assumption that the susceptibility factor would cause greater than 50% of affected sib pairs to share two alleles identical by descent. Theoretical estimates of the power to map susceptibility genes with a high-resolution map of linked markers in a candidate region were made, using HLA as a model locus. This result illustrates the feasibility that IDDM linkage studies using mapped sets of polymorphic DNA markers have, both for other areas of the genome in IDDM and for other polygenic diseases. The analytic approaches introduced here will be useful for affected-sib-pair studies of other complex phenotypes.     

Molecular genetic analysis of the Idd4 locus implicates the IFN response in type 1 diabetes susceptibility in nonobese diabetic mice

High-resolution mapping and identification of the genes responsible for type 1 diabetes (T1D) has proved difficult because of the multigenic etiology and low penetrance of the disease phenotype in linkage studies. Mouse congenic strains have been useful in refining Idd susceptibility loci in the NOD mouse model and providing a framework for identification of genes underlying complex autoimmune syndromes. Previously, we used NOD and a nonobese diabetes-resistant strain to map the susceptibility to T1D to the Idd4 locus on chromosome 11. Here, we report high-resolution mapping of this locus to 1.4 megabases. The NOD Idd4 locus was fully sequenced, permitting a detailed comparison with C57BL/6 and DBA/2J strains, the progenitors of T1D resistance alleles found in the nonobese diabetes-resistant strain. Gene expression arrays and quantitative real-time PCR were used to prioritize Idd4 candidate genes by comparing macrophages/dendritic cells from congenic strains where allelic variation was confined to the Idd4 interval. The differentially expressed genes either were mapped to Idd4 or were components of the IFN response pathway regulated in trans by Idd4. Reflecting central roles of Idd4 genes in Ag presentation, arachidonic acid metabolism and inflammation, phagocytosis, and lymphocyte trafficking, our combined analyses identified Alox15, Alox12e, Psmb6, Pld2, and Cxcl16 as excellent candidate genes for the effects of the Idd4 locus.     

Establishment of the mouse chromosome 7 region with homology to the myotonic dystrophy region of human chromosome 19q

A number of genetic markers, including ATP1A3, TGFB, CKMM, and PRKCG, define the genetic region on human chromosome 19 containing the myotonic dystrophy locus. These and a number of other DNA probes have been mapped to mouse chromosome 7 utilizing a mouse Mus domesticus/Mus spretus interspecific backcross segregating for the genetic markers pink-eye dilution (p) and chinchilla (cch). The establishment of a highly syntenic group conserved between mouse chromosome 7 and human chromosome 19q indicates the likely position of the homologous gene locus to the human myotonic dystrophy gene on proximal mouse chromosome 7. In addition, we have mapped the muscle ryanodine receptor gene (Ryr) to mouse chromosome 7 and demonstrated its close linkage to the Atpa-2, Tgfb-1, and Ckmm cluster of genes. In humans, the malignant hyperthermia susceptibility locus (MHS) also maps close to this gene cluster. The comparative mapping data support Ryr as a candidate gene for MHS.     

A family of type I keratin genes and the homeobox-2 gene complex are closely linked to the rex locus on mouse chromosome 11

Type I and type II keratins are major constituents of intermediate filaments that play a fundamental role in the cytoskeletal network. By using both somatic cell hybrids and conventional and interspecific linkage crosses, several genes encoding type I keratins, including the epidermal keratin K10, were shown to be closely linked to the homeobox-2 complex and the rex locus on mouse chromosome 11. The absence of crossovers between type I keratin-encoding genes and rex (N = 239), a locus affecting hair development, raises the possibility that mutations at rex and neighboring loci affecting skin and hair development involve type I keratin genes.     

09/15: Comparative genomics of a conserved chromosomal region associated with a complex human phenotype

Three genes that encode related immunoglobulin superfamily molecules have recently been mapped to human chromosome 15 in the region q22.3-q23 and to the syntenic region on mouse chromosome 9. These genes presumably derived from gene duplications, and they are highly similar to Deleted in Colorectal Cancer (DCC), which functions as an axon guidance molecule during development of the nervous system. To find out whether additional genes of this class were present in a chromosomal cluster, we produced a comparative physical map within the region of synteny between mouse chromosome 9 and human chromosome 15. This interval overlaps the critical region for the fourth genetic locus for Bardet-Biedl syndrome (BBS4) in humans. Bardet-Biedl syndrome (OMIM 600374) is characterized by poly/syn/brachydactyly, retinal degeneration, hypogonadism, mental retardation, obesity, diabetes, and kidney abnormalities. A detailed map of this locus will help to identify candidate genes for this disorder.     

Mouse DESC1 is located within a cluster of seven DESC1-like genes and encodes a type II transmembrane serine protease that forms serpin inhibitory complexes

We report the identification and functional analysis of a type II transmembrane serine protease encoded by the mouse differentially expressed in squamous cell carcinoma (DESC) 1 gene, and the definition of a cluster of seven homologous DESC1-like genes within a 0.5-Mb region of mouse chromosome 5E1. This locus is syntenic to a region of human chromosome 4q13.3 containing the human orthologues of four of the mouse DESC1-like genes. Bioinformatic analysis indicated that all seven DESC1-like genes encode functional proteases. Direct cDNA cloning showed that mouse DESC1 encodes a multidomain serine protease with an N-terminal signal anchor, a SEA (sea urchin sperm protein, enterokinase, and agrin) domain, and a C-terminal serine protease domain. The mouse DESC1 mRNA was present in epidermal, oral, and male reproductive tissues and directed the translation of a membrane-associated 60-kDa N-glycosylated protein with type II topology. Mouse DESC1 was synthesized in insect cells as a zymogen that could be activated by exposure to trypsin. The purified activated DESC1 hydrolyzed synthetic peptide substrates, showing a preference for Arg in the P1 position. DESC1 proteolytic activity was abolished by generic inhibitors of serine proteases but not by other classes of protease inhibitors. Most interestingly, DESC1 formed stable inhibitory complexes with both plasminogen activator inhibitor-1 and protein C inhibitor that are expressed in the same tissues with DESC1, suggesting that type II transmembrane serine proteases may be novel targets for serpin inhibition. Together, these data show that mouse DESC1 encodes a functional cell surface serine protease that may have important functions in the epidermis, oral, and reproductive epithelium.