ABERRATIONS OF A PUTATIVE TUMOR SUPPRESSOR GENE SEL1L IN PANCREATIC DUCTAL ADENOCARCINOMA

Introduction: Pancreatic cancer is the fourth leading cause of cancer-related death among males and females in the United States. Sel-1-like (SEL1L) is a putative tumor suppressor gene that is downregulated in a significant proportion of human pancreatic ductal adenocarcinoma (PDAC). It was hypothesized that SEL1L expression could be down-modulated by somatic mutation, loss of heterozygosity (LOH), CpG island hypermethylation and/or aberrantly expressed microRNAs (miRNAs). Material and methods: In 42 PDAC tumors, the SEL1L coding region was amplified using reverse transcription polymerase chain reaction (RT-PCR), and analyzed by agarose gel electrophoresis and sequenced to search for mutations. Using fluorescent fragment analysis, two intragenic microsatellites in the SEL1L gene region were examined to detect LOH in a total of 73 pairs of PDAC tumors and normal-appearing adjacent tissues. Bisulfite DNA sequencing was performed to determine the methylation status of the SEL1L promoter in 41 PDAC tumors and 6 PDAC cell lines. Using real-time quantitative PCR, the expression levels of SEL1L mRNA and 7 aberrantly upregulated miRNAs that potentially target SEL1L were assessed in 42 PDAC tumor and normal pairs. Statistical methods were applied to evaluate the correlation between SEL1L mRNA and the miRNAs. Further the interaction was determined by functional analysis using a molecular biological approach. Results: No mutations were detected in the SEL1L coding region. More than 50% of the samples displayed abnormally alternate or aberrant spliced transcripts of SEL1L. About 14.5% of the tumors displayed LOH at the CAR/CAL microsatellite locus and 10.7% at the RepIN20 microsatellite locus. However, the presence of LOH did not show significant association with SEL1L downregulation. No methylation was observed in the SEL1L promoter. Statistical analysis showed that SEL1L mRNA expression levels significantly and inversely correlated with the expression of hsa-mir-143, hsa-mir-155, and hsa-mir-223. Functional analysis indicated that hsa-mir-155 acted as a suppressor of SEL1L in PL18 and MDAPanc3 PDAC cell lines. Discussion: Evidence from these studies suggested that SEL1L was possibly downregulated by aberrantly upregulated miRNAs in PDAC. Future studies should be directed towards developing a better understanding of the mechanisms for generation of aberrant SEL1L transcripts, and further analysis of miRNAs that may downregulate SEL1L.

A novel myeloid leukemia suppressor locus at human chromosome 5q13.3

Molecular mechanisms that underlie preleukemic myelodysplasia (MDS) and acute myelogenous leukemia (AML) are poorly understood. In MDS or AML with a refractory clinical course, more than 30% of patients have acquired interstitial or complete deletions of chromosome 5. The 5q13.3 chromosomal segment is commonly lost as the result of 5q deletion. Reciprocal and unbalanced translocations of 5q13.3 can also occur as sole anomalies associated with refractory AML or MDS. This study addresses the hypothesis that a critical gene at 5q13.3 functions either as a classical tumor suppressor or as a chromosomal translocation partner and contributes to leukemogenesis. ^ Previous studies from our laboratory delineated a critical region of loss to a 2.5–3.0Mb interval at 5q13.3 between microsatellite markers D5S672 and GATA-P18104. The critical region of loss was later resolved to an interval of approximately 2Mb between the markers D5S672 and D5S2029. I, then generated a long range physical map of yeast artificial chromosomes (YACs) and developed novel sequence tagged sites (STS). To enhance the resolution of this map, bacterial artificial chromosomes (BACs) were used to construct a triply linked contig across a 1 Mb interval. These BACs were used as probes for fluorescent in situ hybridization (FISH) on an AML cell line to define the 5q13.3 critical region. A 200kb BAC, 484a9, spans the translocation breakpoint in this cell line. A novel gene, SSDP2 (single stranded DNA binding protein), is disrupted at the breakpoint because its first four exons are encoded within 140kb of BAC 484a9. This finding suggests that SSDP2 is the critical gene at 5q13.3. ^ In addition, I made an observation that deletions of chromosome 5q13 co-segregate with loss of the chromosome 17p. In some cases the deletions result from unbalanced translocations between 5q13 and 17p13. It was confirmed that the TP53 gene is deleted in patients with 17p loss, and the remaining allele harbors somatic mutation. Thus, the genetic basis for the aggressive clinical course in AML and MDS may be caused by functional cooperation between deletion or disruption of the 5q13.3 critical gene and inactivation of TP53. ^

Physical and functional mapping of a novel tumor suppressor locus for prostate cancer within chromosome 10p12.31-q11

Prostate cancer remains the second leading cause of male cancer deaths in the United States, yet the molecular mechanisms underlying this disease remain largely unknown. Cytogenetic and molecular analyses of prostate tumors suggest a consistent association with the loss of chromosome 10. Previously, we have defined a novel tumor suppressor locus PAC-1 within chromosome 10pter-q11. Introduction of the short arm of chromosome 10 into a prostatic adenocarcinoma cell line PC-3H resulted in dramatic tumor suppression and restoration of a programmed cell death pathway. Using a combined approach of comparative genomic hybridization and microsatellite analysis of PC-3H, I have identified a region of hemizygosity within 10p12-p15. This region has been shown to be involved in frequent loss of heterozygosity in gliomas and melanoma. To functionally dissect the region within chromosome 10p containing PAC-1, we developed a strategy of serial microcell fusion, a technique that allows the transfer of defined fragments of chromosome 10p into PC-3H. Serial microcell fusion was used to transfer defined 10p fragments into a mouse A9 fibrosarcoma cell line. Once characterized by FISH and microsatellite analyses, the 10p fragments were subsequently transferred into PC-3H to generate a panel of microcell hybrid clones containing overlapping deletions of chromosome 10p. In vivo and microsatellite analyses of these PC hybrids identified a small chromosome 10p fragment (an estimated 31 Mb in size inclusive of the centromere) that when transferred into the PC-3H background, resulted in significant tumor suppression and limited a region of functional tumor suppressor activity to chromosome 10p12.31-q11. This region coincides with a region of LOH demonstrated in prostate cancer. These studies demonstrate the utility of this approach as a powerful tool to limit regions of functional tumor suppressor activity. Furthermore, these data used in conjunction with data generated by the Human Genome Project lent a focused approach to identify candidate tumor suppressor genes involved in prostate cancer. ^

Genome -wide and locus -specific instability

The discovery of expanded simple repeated sequences causing or associated with human disease has lead to a new area of research involved in the elucidation of how the expanded repeat causes disease and how the repeat becomes unstable. ^ To study the genetic basis of the (CTG)n repeat instability in the DMPK gene in myotonic dystrophy (DM1) patients, somatic cell hybrids were constructed between the lymphocytes of DM1 patients and a variety of Chinese hamster ovary (CHO) cell DNA repair gene deficient mutants. By using small pool PCR (SP-PCR), the instability of the (CTG)n can be quantitated for both the frequency and sizes of length change mutations. ^ Additional SP-PCR analysis on 2/11 subclones generated from this original hybrid showed a marked increase in large repeat deletions, ∼50%. A bimodal distribution of repeats was seen around the progenitor allele and at a large deleted product (within the normal range) with no intermediate products present. ^ To determine if the repair capacity of the CHO cell led to a mutator phenotype in the hamster and hybrid clones, SP-PCR was also done on 3 hamster microsatellites in a variety of hamster cell backgrounds. No variant alleles were seen in over 2500 genome equivalents screened. ^ Human-hamster hybrids have long been shown to be chromosomally unstable, yet information about the stability of repeated sequences was not known. To test if repeat instability was associated with either intact or non-intact human chromosomes, more than 300 microsatellite repeats on 13 human chromosomes (intact and non-intact) were analyzed in eight hybrid cells. No variants were seen between the hybrid and patient alleles in the hybrids. ^ To identify whether DM1 patients have a previously undetected level of genome wide instability or if the instability is truly locus specific, SP-PCR was done on 6 human microsatellites within the patient used to make the hybrid cells. No variants were seen in over 1000 genomes screened. ^ These studies show that the somatic cell hybrid approach is a genetically stable system that allows for the determination of factors that could lead to changes in microsatellite instability. It also shows that there is something inherent about the DM1 expanded (CTG)n repeat that it is solely targeted by, as of yet, and unknown mechanism that causes the repeat to be unstable. (Abstract shortened by UMI.)^

CHARACTERIZING A NOVEL GENETIC LOCUS ASSOCIATED WITH FAMILIAL CO-OCCURRENCE OF THORACIC AORTIC ANEURYSMS AND INTRACRANIAL ANEURYSMS

The Mendelian inheritance of genetic mutations can lead to adult-onset cardiovascular disease. Several genetic loci have been mapped for the familial form of Thoracic Aortic Aneurysms (TAA), and many causal mutations have been identified for this disease. Intracranial Aneurysms (ICA) also show linkage heterogeneity, but no mutations have been identified causing familial ICA alone. Here, we characterized a large family (TAA288) with an autosomal dominant pattern of inherited aneurysms. It is intriguing that female patients predominantly present with ICA and male patients predominantly with TAA in this family. To identify a causal mutation in this family, a genome-wide linkage analysis was previously performed on nine members of this family using the 50k GenChips Hind array from Affymetrix. This analysis eventually identified a single disease-segregating locus, on chromosome 5p15. We build upon this previous analysis in this study, hypothesizing that a genetic mutation inherited in this locus leads to the sex-specific phenotype of TAA and ICA in this family First we refined the boundaries of the 5p15 disease linked locus down to the genomic coordinates 5p15: 3,424,465- 6,312,925 (GRCh37/hg19 Assembly). This locus was named the TAA288 critical interval. Next, we sequenced candidate genes within the TAA288 critical interval. The selection of genes was simplified by the relatively small number of well-characterized genetic elements within the region. Seeking novel or rare disease-segregating variants, we initially observed a single point alteration in the metalloproteinase gene ADAMTS16 fulfilling this criteria. This variant was later classified as a low-frequency population polymorphism (rs72647757), but we continued to explore the potential role of the ADAMTS16 as the cause of disease in TAA288. We observed that fibroblasts cultured from TAA288 patients consistently upregulated the expression of this gene more strongly compared to matched control fibroblasts when treated with the cytokine TGF-β1, though there was some variation in the exact nature of this expression. We also observed evidence that this protein is expressed at elevated levels in aortic aneurysm tissue from patients with mutations in the gene TGFBR2 and Marfan syndrome, shown by immunohistochemical detection of this protein.

Bicarbonate enhances expression of the endocarditis and biofilm associated pilus locus, ebpR-ebpABC, in Enterococcus faecalis.

BACKGROUND: We previously identified ebpR, encoding a potential member of the AtxA/Mga transcriptional regulator family, and showed that it is important for transcriptional activation of the Enterococcus faecalis endocarditis and biofilm associated pilus operon, ebpABC. Although ebpR is not absolutely essential for ebpABC expression (100-fold reduction), its deletion led to phenotypes similar to those of an ebpABC mutant such as absence of pili at the cell surface and, consequently, reduced biofilm formation. A non-piliated ebpABC mutant has been shown to be attenuated in a rat model of endocarditis and in a murine urinary tract infection model, indicating an important participation of the ebpR-ebpABC locus in virulence. However, there is no report relating to the environmental conditions that affect expression of the ebpR-ebpABC locus. RESULTS: In this study, we examined the effect of CO2/HCO3(-), pH, and the Fsr system on the ebpR-ebpABC locus expression. The presence of 5% CO2/0.1 M HCO3(-) increased ebpR-ebpABC expression, while the Fsr system was confirmed to be a weak repressor of this locus. The mechanism by which the Fsr system repressed the ebpR-ebpABC locus expression appears independent of the effects of CO2(-) bicarbonate. Furthermore, by using an ebpA::lacZ fusion as a reporter, we showed that addition of 0.1 M sodium bicarbonate to TSBG (buffered at pH 7.5), but not the presence of 5% CO2, induced ebpA expression in TSBG broth. In addition, using microarray analysis, we found 73 genes affected by the presence of sodium bicarbonate (abs(fold) > 2, P < 0.05), the majority of which belong to the PTS system and ABC transporter families. Finally, pilus production correlated with ebpA mRNA levels under the conditions tested. CONCLUSIONS: This study reports that the ebp locus expression is enhanced by the presence of bicarbonate with a consequential increase in the number of cells producing pili. Although the molecular basis of the bicarbonate effect remains unclear, the pathway is independent of the Fsr system. In conclusion, E. faecalis joins the growing family of pathogens that regulates virulence gene expression in response to bicarbonate and/or CO2.

Importance of the ebp (endocarditis- and biofilm-associated pilus) locus in the pathogenesis of Enterococcus faecalis ascending urinary tract infection.

BACKGROUND: We recently demonstrated that the ubiquitous Enterococcus faecalis ebp (endocarditis- and biofilm-associated pilus) operon is important for biofilm formation and experimental endocarditis. Here, we assess its role in murine urinary tract infection (UTI) by use of wild-type E. faecalis OG1RF and its nonpiliated, ebpA allelic replacement mutant (TX5475). METHODS: OG1RF and TX5475 were administered transurethrally either at an ~1 : 1 ratio (competition assay) or individually (monoinfection). Kidney pairs and urinary bladders were cultured 48 h after infection. These strains were also tested in a peritonitis model. RESULTS: No differences were observed in the peritonitis model. In mixed UTIs, OG1RF significantly outnumbered TX5475 in kidneys (P=.0033) and bladders (P< or =.0001). More OG1RF colony-forming units were also recovered from the kidneys of monoinfected mice at the 4 inocula tested (P=.015 to P=.049), and 50% infective doses of OG1RF for kidneys and bladder (9.1x10(1) and 3.5x10(3) cfu, respectively) were 2-3 log(10) lower than those of TX5475. Increased tropism for the kidney relative to the bladder was observed for both OG1RF and TX5475. CONCLUSION: The ebp locus, part of the core genome of E. faecalis, contributes to infection in an ascending UTI model and is the first such enterococcal locus shown to be important in this site.

Origins and evolution of VNTR loci: The apolipoprotein B 3$\sp\prime$ VNTR

I studied the apolipoprotein (apo) B 3$\sp\prime$ variable number tandem repeat (VNTR) and did computer simulations of the stepwise mutation model to address four questions: (1) How did the apo B VNTR originate? (2) What is the mutational mechanism of repeat number change at the apo B VNTR? (3) To what extent are population and molecular level events responsible for the determination of the contemporary apo B allele frequency distribution? (4) Can VNTR allele frequency distributions be explained by a simple and conservative mutation-drift model? I used three general approaches to address these questions: (1) I characterized the apo B VNTR region in non-human primate species; (2) I constructed haplotypes of polymorphic markers flanking the apo B VNTR in a sample of individuals from Lorrain, France and studied the associations between the flanking-marker haplotypes and apo B VNTR size; (3) I did computer simulations of the one-step stepwise mutation model and compared the results to real data in terms of four allele frequency distribution characteristics.^ The results of this work have allowed me to conclude that the apo B VNTR originated after an initial duplication of a sequence which is still present as a single copy sequence in New World monkey species. I conclude that this locus did not originate by the transposition of an array of repeats from somewhere else in the genome. It is unlikely that recombination is the primary mutational mechanism. Furthermore, the clustered nature of these associations implicates a stepwise mutational mechanism. From the high frequencies of certain haplotype-allele size combinations, it is evident that population level events have also been important in the determination of the apo B VNTR allele frequency distribution. Results from computer simulations of the one-step stepwise mutation model have allowed me to conclude that bimodal and multimodal allele frequency distributions are not unexpected at loci evolving via stepwise mutation mechanisms. Short tandem repeat loci fit the stepwise mutation model best, followed by microsatellite loci. I therefore conclude that there are differences in the mutational mechanisms of VNTR loci as classed by repeat unit size. (Abstract shortened by UMI.) ^

A molecular study of the rhodopsin locus in patients with retinitis pigmentosa

DNA for this study was collected from a sample of 133 retinitis pigmentosa (RP) patients and the rhodopsin locus molecularly analyzed by linkage and for disease specific mutations. The cohort of patients consisted of 85 individuals diagnosed with autosomal dominant RP (adRP), and 48 patients representing other forms of retinitis pigmentosa or retinal dystrophy related disease. In three large families with adRP rhodopsin was excluded from linkage to the disease locus. A search for subtle mutations in the rhodopsin coding region using single strand conformational polymorphisms (SSCP) and sequencing detected a total of 14 unique sequence variants in 24 unrelated patients. These variants included one splicing variant, 5168 -1G-A, one deletion variant of 17 base pairs causing a frame shift at codon 332, and 12 misense variants: Pro23His, Leu46Arg, Gly106Trp, Arg135Pro, Pro171Glu, Pro180Ala, Glu181Lys, Asp190Asn, His211Arg, Ser270Arg, Leu328Pro and Pro347Thr. All but three of the missense variants change amino acids that are evolutionarily conserved. The Pro23His mutation was found in 10 unrelated individuals with family histories of adRP and not in any normal controls (over 80 chromosomes tested). The Pro180Ala mutation was present in a patient with simplex RP and probably represents a new mutation. Three normal polymorphic nucleotide substitutions, A-269-G, T-3982-C, and G-5145-A, were also identified. We conclude, based on this study, that 25% of adRP cases are attributable to rhodopsin mutations.^ Clinical data, including ERG results and visual field testing, was available for patients with eleven different mutations. The eleven patients were all diagnosed with RP, however the severity of the disease varied with five patients mildly affected and diagnosed with type II adRP and 5 patients severely affected and diagnosed with type I adRP. The patient with simplex RP was mildly affected. The location of the mutations within the rhodopsin protein was randomly associated with the severity of the disease in those patients evaluated. However, four mutations, Pro23His, Leu46Arg, Pro347Thr, and 5168 -1G-A, are particularly interesting. The Pro23His mutation appears to have radiated from a recent common ancestor of the affected patients as all of them share a common haplotype at the rhodopsin locus. The Leu46Arg mutation causes an unusually severe form of RP. Hydropathy analysis of the mutated sequence revealed a marked change in the hydrophobicity of this first transmembrane spanning region. Codon 347 has been the target of multiple mutations with at least six documented changes at the position, significantly more than expected by a random distribution of mutations. Finally the splice-site variant is extremely variable in its expression in the family studied. Similar mutations have been reported in other cases of adRP and postulated to be involved in autosomal recessive RP (arRP). Mechanisms to account for the variable expression of rhodopsin mutations in relation to RP heterogeneity are discussed. (Abstract shortened by UMI.) ^

Mitotic mapping of a novel tumor suppressor locus on human chromosome 3q important in osteosarcoma tumorigenesis

Tumor-specific loss of constitutional heterozygosity by deletion, mitotic recombination or nondisjunction is a common mechanism for tumor suppressor allele inactivation. When loss of heterozygosity is the result of mitotic recombination, or a segmental deletion event, only a portion of the chromosome is lost. This information can be used to map the location of new tumor suppressor genes. In osteosarcoma, the highest frequencies of loss of heterozygosity have been reported for chromosomes 3q, 13q, 17p. On chromosomes 13q and 17p, allelic losses are associated with loss of function at the retinoblastoma susceptibility locus (RB1) and the p53 locus, respectively. Chromosome 3q is also of particular interest because the high percent of loss of heterozygosity (62%-75%) suggests the presence of another tumor suppressor important for osteosarcoma tumorigenesis. To localize this putative tumor suppressor gene, we used polymorphic markers on chromosome 3q to find the smallest common region of allele loss. This putative tumor suppressor was localized to a 700 kb region on chromosome 3q26.2 between the polymorphic loci D3S1282 and D3S1246. ^

Characterization and subchromosomal location of a novel tumor suppressor gene on human chromosome 7

Alterations in oncogenes and tumor suppressor genes (TSGs) are considered to be critical steps in oncogenesis. Consistent deletions and loss of heterozygosity (LOH) of polymorphic markers in a determinate chromosomal fragment are known to be indicative of a closely mapping TSG. Deletion of the long arm of chromosome 7 (hchr 7) is a frequent trait in many kinds of human primary tumors. LOH was studied with an extensive set of markers on chromosome 7q in several types of human neoplasias (primary breast, prostate, colon, ovarian and head and neck carcinomas) to determine the location of a putative TSG. The extent of LOH varied depending the type of tumor studied but all the LOH curves we obtained had a peak at (C-A)$\sb{\rm n}$ microsatellite repeat D7S522 at 7q31.1 and showed a Gaussian distribution. The high incidence of LOH in all tumor types studied suggests that a TSG relevant to the development of epithelial cancers is present on the 7q31.1. To investigate whether the putative TSG is conserved in the syntenic mouse locus, we studied LOH of 30 markers along mouse chromosome 6 (mchr 6) in chemically induced squamous cell carcinomas (SCCs). Tumors were obtained from SENCAR and C57BL/6 x DBA/2 F1 females by a two-stage carcinogenesis protocol. The high incidence of LOH in the tumor types studied suggests that a TSG relevant to the development of epithelial cancers is present on mchr 6 A1. Since this segment is syntenic with the hchr 7q31, these data indicate that the putative TSG is conserved in both species. Functional evidence for the existence of a TSG in hchr 7 was obtained by microcell fusion transfer of a single hchr 7 into a murine SCC-derived cell line. Five out of seven hybrids had two to three-fold longer latency periods for in vivo tumorigenicity assays than parental cells. One of the unrepressed hybrids had a deletion in the introduced chromosome 7 involving q31.1-q31.3, confirming the LOH data. ^

Linkage and association of candidate obesity genes in Mexican-Americans from Starr County, Texas

Obesity and related chronic diseases represent a tremendous public health burden among Mexican Americans, a young and rapidly-expanding population. This study investigated the impact of variation within eight candidate obesity genes, which include leptin (LEP), leptin receptor (LEPR), neuropeptide Y (NPY), NPYY1 receptor (NPYY1), glucagon-like peptide-1 (GLP-1), GLP-1 receptor (GLP1R), beta-3 adrenergic receptor (β3AR), and uncoupling protein (UCP1), on variation in human obesity status and/or quantitative traits related to obesity in Mexican Americans from Starr County, Texas. The Trp64Arg polymorphism within β3AR was typed in 820 random individuals and 240 pedigrees (N = 2,044). The Arg allele frequency was significantly greater in obese versus non-obese individuals (0.20 versus 0. 15, respectively). In addition, within the random sample, the Arg allele was associated with significantly greater body weight (p = 0.031) and body mass index (BMI, p = 0.008) than the Trp allele. In the family sample, the Trp64Arg locus was also linked to percent fat (p = 0.045) but not to body weight or BMI. No linkage between obesity, diabetes, hypertension, or gallbladder disease and the Trp64Arg mutation was observed in families using affected sib pair linkage analysis or the transmission disequilibrium test. Microsatellite markers proximate to the remaining seven genes were typed in 302 individuals from 59 families. Sib pair linkage analysis provided evidence for linkage between obesity and NPY within affected sibling pairs (p = 0.042; n = 170 pairs). NPY was also linked to weight (p = 0.020), abdominal circumference (p = 0.031), hip circumference (p = 0.012), DBP (p ≤ 0.005), and a composite measure of body mass/fat (p ≤ 0.048) in all sibling pairs (n = 545 pairs). Additionally, LEP was linked to waist/hip ratio (p ≤ 0.009), total cholesterol (p ≤ 0.030), and HDL cholesterol (p ≤ 0.026), and LEPR was linked to fasting blood glucose (p ≤ 0.018) and DBP (p ≤ 0.003). Subsequent to the linkage analyses, the NPY gene was sequenced and eight variant sites identified. Two variant sites (-880I/D and 69I/D) were typed in a random sample of 914 individuals. The 880I/D variant was significantly associated with waist/hip ratio (p = 0.035) in the entire sample (N = 914) and with BMI (p = 0. 031), abdominal circumference (p = 0.044), and waist/hip ratio (p = 0.041) in a non-obese subsample (BW < 30 kg/m2, n = 594). The 69I/D variant was a rare mutation observed in only one pedigree and was not associated with obesity or body size/mass within this pedigree. Results of this study indicate that variation at or near β3AR, LEP, LEPR, and NPY may exert effects which increase obesity susceptibility and influence obesity-related measures in this population. ^

Population and molecular evolution studies on the Melanocortin 1 receptor locus implicate its role in human pigmentation variation

Human pigmentation is a complex trait with the observed variation caused by the varied production of eumelanin (brown/black melanins) and phaeomelanin (red/yellow melanins) by the melanocytes. The melanocortin 1 receptor (MC1R), a G protein-coupled receptor expressed in the melanocytes, is a regulator eu- and phaeomelanin synthesis, and MC1R mutations causing skin and coat color changes are known in many mammals. To understand the role of MC1R in human pigmentation variation, I have sequenced the MC1R gene in 121 individuals sampled from world populations. In addition, I have sequenced the MC1R gene in common and pygmy chimpanzees, gorilla, orangutan, and baboon to study the evolution of MC1R and to infer the ancestral human MC1R sequence. The ancestral MC1R sequence is observed in all 25 African individuals studied, but at lower frequencies in the other populations examined, especially in East and Southeast Asians. The Arg163Gln variant is absent in the Africans studied, almost absent in Europeans, and at a low frequency in Indians, but is at an exceptionally high frequency (70%) in East and Southeast Asians. To further evaluate the role of MC1R variants in human pigmentation variation, I have combined these molecular evolution and population studies with functional assays on MC1R variants and primate MC1Rs. ^

Detecting the signature of natural selection with microsatellites

Natural selection is one of the major factors in the evolution of all organisms. Detecting the signature of natural selection has been a central theme in evolutionary genetics. With the availability of microsatellite data, it is of interest to study how natural selection can be detected with microsatellites. ^ The overall aim of this research is to detect signatures of natural selection with data on genetic variation at microsatellite loci. The null hypothesis to be tested is the neutral mutation theory of molecular evolution, which states that different alleles at a locus have equivalent effects on fitness. Currently used tests of this hypothesis based on data on genetic polymorphism in natural populations presume that mutations at the loci follow the infinite allele/site models (IAM, ISM), in the sense that at each site at most only one mutation event is recorded, and each mutation leads to an allele not seen before in the population. Microsatellite loci, which are abundant in the genome, do not obey these mutation models, since the new alleles at such loci can be created either by contraction or expansion of tandem repeat sizes of core motifs. Since the current genome map is mainly composed of microsatellite loci and this class of loci is still most commonly studied in the context of human genome diversity, this research explores how the current test procedures for testing the neutral mutation hypothesis should be modified to take into account a generalized model of forward-backward stepwise mutations. In addition, recent literature also suggested that past demographic history of populations, presence of population substructure, and varying rates of mutations across loci all have confounding effects for detecting signatures of natural selection. ^ The effects of the stepwise mutation model and other confounding factors on detecting signature of natural selection are the main results of the research. ^