Global patterns of DNA polymorphism in noncoding regions in human populations

DNA sequence variation is currently a major source of data for studying human origins, evolution, and demographic history, and for detecting linkage association of complex diseases. In this dissertation, I investigated DNA variation in worldwide populations from two ∼10 kb autosomal regions on 22q11.2 (noncoding) and 1q24 (introns). A total of 75 variant sites were found among 128 human sequences in the 22q11.2 region, yielding an estimate of 0.088% for nucleotide diversity (π), and a total of 52 variant sites were found among 122 human sequences in the 1q24 region with an estimated π value of 0.057%. The data from these two regions and a 10 kb noncoding region on Xq13.3 all show a strong excess of low-frequency variants in comparison to that expected from an equilibrium population, indicating a relatively recent population expansion. The effective population sizes estimated from the three regions were 11,000, 12,700, and 8,600, respectively, which are close to the commonly used value of 10,000. In each of the two autosomal regions, the age of the most recent common ancestor (MRCA) was estimated to be older than 1 million years among all the sequences and ∼600,000 years among non-African sequences, providing first evidence from autosomal noncoding or intronic regions for a genetic history of humans much more ancient than the emergence of modern humans. The ancient genetic history of humans indicates no severe bottleneck during the evolution of humans in the last half million years; otherwise, much of the ancient genetic history would have been lost during a severe bottleneck. This study strongly suggests that both the “out of Africa” and the multiregional models are too simple for explaining the evolution of modern humans. A compilation of genome-wide data revealed that nucleotide diversity is highest in autosomal regions, intermediate in X-linked regions, and lowest in Y-linked regions. The data suggest the existence of background selection or selective sweep on Y-linked loci. In general, the nucleotide diversity in humans is low compared to that in chimpanzee and Drosophila populations. ^

Molecular evolution of primate immunodeficiency viruses and hepatitis delta virus

Primate immunodeficiency viruses, or lentiviruses (HIV-1, HIV-2, and SIV), and hepatitis delta virus (HDV) are RNA viruses characterized by rapid evolution. Infection by primate immunodeficiency viruses usually results in the development of acquired immunodeficiency syndrome (AIDS) in humans and AIDS-like illnesses in Asian macaques. Similarly, hepatitis delta virus infection causes hepatitis and liver cancer in humans. These viruses are heterogeneous within an infected patient and among individuals. Substitution rates in the virus genomes are high and vary in different lineages and among sites. Methods of phylogenetic analysis were applied to study the evolution of primate lentiviruses and the hepatitis delta virus. The following results have been obtained: (1) The substitution rate varies among sites of primate lentivirus genes according to the two parameter gamma distribution, with the shape parameter $\alpha$ being close to 1. (2) Primate immunodeficiency viruses fall into species-specific lineages. Therefore, viral transmissions across primate species are not as frequent as suggested by previous authors. (3) Primate lentiviruses have acquired or lost their pathogenicity several times in the course of evolution. (4) Evidence was provided for multiple infections of a North American patient by distinct HIV-1 strains of the B subtype. (5) Computer simulations indicate that the probability of committing an error in testing HIV transmission depends on the number of virus sequences and their length, the divergence times among sequences, and the model of nucleotide substitution. (6) For future investigations of HIV-1 transmissions, using longer virus sequences and avoiding the use of distant outgroups is recommended. (7) Hepatitis delta virus strains are usually related according to the geographic region of isolation. (8) Evolution of HDV is characterized by the rate of synonymous substitution being lower than the nonsynonymous substitution rate and the rate of evolution of the noncoding region. (9) There is a strong preference for G and C nucleotides at the third codon positions of the HDV coding region. ^

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. ^

Prognostic significance of xCT polymorphisms and expression in patients with advanced pancreatic cancer treated with chemotherapy

The plasma membrane xc- cystine/glutamate transporter mediates cellular uptake of cystine in exchange for intracellular glutamate and is highly expressed by pancreatic cancer cells. The xCT gene, encoding the cystine-specific xCT protein subunit of xc-, is important in regulating intracellular glutathione (GSH) levels, critical for cancer cell protection against oxidative stress, tumor growth and resistance to chemotherapeutic agents including platinum. We examined 4 single nucleotide polymorphisms (SNPs) of the xCT gene in 269 advanced pancreatic cancer patients who received first line gemcitabine with or without cisplatin or oxaliplatin. Genotyping was performed using Taqman real-time PCR assays. A statistically significant correlation was noted between the 3' untranslated region (UTR) xCT SNP rs7674870 and overall survival (OS): Median survival time (MST) was 10.9 and 13.6 months, respectively, for the TT and TC/CC genotypes (p = 0.027). Stratified analysis showed the genotype effect was significant in patients receiving gemcitabine in combination with platinum therapy (n = 145): MST was 10.5 versus 14.1 months for the TT and TC/CC genotypes, respectively (p = 0.013). The 3' UTR xCT SNP rs7674870 may correlate with OS in pancreatic cancer patients receiving gemcitabine and platinum combination therapy. Paraffin-embedded core and surgical biopsy tumor specimens from 98 patients with metastatic pancreatic adenocarcinoma were analyzed by immunohistochemistry using an xCT specific antibody. xCT protein IHC expression scores were analyzed in relation to overall survival in 86 patients and genotype in 12 patients and no statistically significant association was found between the level of xCT IHC expression score and overall survival (p = 0.514). When xCT expression was analyzed in terms of treatment response, no statistically significant associations could be determined (p = 0.908). These data suggest that polymorphic variants of xCT may have predictive value, and that the xc- transporter may represent an important target for therapy in pancreatic cancer.

EXPRESSION OF THE HUMAN PLACENTAL LACTOGEN GENES

The placenta is the site of synthesis of various peptide and steroid hormones related to pregnancy. Human placental lactogen (hPL) is the predominant peptide hormone secreted by term placenta and its synthesis is tissue-specific and coupled to placenta development. The objective of this work was to study the structure and expression of the hPL.^ Poly(A('+))RNA from human term placenta was translated in a mouse-derived cell-free system. A major band corresponding to pre-hPL and a minor band comigrating with mature hPL, represent (TURN)15% of the total radioactively labeled proteins. Analysis of the poly(A('+))RNA showed a prominent band at approximately 860 nucleotides. A corresponding band was observed in Northern blots of total RNA, hybridized with {('32)P}-labeled recombinant plasmid containing a portion of hPL cDNA. Similar analyses of nuclear RNA showed at least four additional bands at 990, 1200, 1460 and 1760 nucleotides, respectively, which are likely precursors of hPL mRNA. Poly(A('+))RNA was used to construct a cDNA library, of which approximately 5% of the clones were found to hybridize to hPL DNA sequences. Heteroduplexes constructed between a clone containing a 815 bp hPL cDNA insert and a hPL genomic DNA clone revealed four small intervening sequences which can account for the lengths observed in hnRNA molecules.^ Recombinant plasmid HCS-pBR322 containing a 550 bp insert of a cDNA transcript of human placental lactogen (hPL) mRNA was ('3)H-labeled an hybridized in situ to human chromosome preparations. These experiments allowed assignment of the hPL and growth hormone (hGH) genes, which have over 90% nucleotide homology in their coding sequences, to band q22-24 of chromosome 17. A gene copy number experiment showed that both genes are present in (TURN)3 copies per haploid genome.^ Experiments were designed to determine if all members of the hPL gene cluster, consisting of four non-allelic genes, are transcribed in term placenta. Advantage was taken of differences in restriction endonuclease sites in the coding portions of the different hPL genes, to distinguish the putative cDNAs of the transcriptionally active genes. Two genes were found to be represented in the cDNA library and their cDNA transcripts were isolated and characterized. Three independent methods showed that their corresponding mRNAs are about equally represented in the hPL mRNA population. The two cDNAs code for prehPL proteins which differ at a single amino acid position. However the secreted hPLs have identical amino acid sequences. A tetramer insertion duplication was found in a palindrome area of the 3' untranslated region of one of the hPL mRNAs. ^

GLUCAGON GENE EXPRESSION: ANALYSIS OF PREPROGLUCAGON

Glucagon is a 29 amino acid polypeptide hormone produced in the (alpha) cells of the pancreatic islets. The purpose of this research was to understand better the role of glucagon in the regulation of metabolic processes. As with other polypeptide hormones, the synthesis of glucagon is thought to involve a larger precursor, which is then enzymatically cleaved to the functional form. The specific research objectives were to obtain cloned copies of the messenger RNA (mRNA) for pancreatic glucagon, to determine their primary sequences, and from this coding information to deduce the amino acid sequence of the initial glucagon precursor. From this suggested preproglucagon sequence and prior information on possible proglucagon intermediate processing products, the overall objective of this research is to propose a possible pathway for the biosynthesis of pancreatic glucagon.^ Synthetic oligodeoxynucleotide probes of 14-nucleotides (14-mer) and 17-nucleotides (a 17-mer) complementary to codons specifying a unique sequence of mature glucagon were synthesized. The ('32)P-labeled-14-mer was hybridized with size-fractionated fetal bovine pancreatic poly(A('+))RNA bound to nitrocellulose. RNA fractions of (TURN)14S were found to hybridize specifically, resulting in an (TURN)10-fold enrichment for these sequences. These poly(A('+))RNAs were translated in a cell-free system and the products analyzed by gel electrophoresis. The translation products were found to be enriched for a protein of the putative size of mammalian preproglucagon ((TURN)21 kd). These enriched RNA fractions were used to construct a complementary DNA (cDNA) library is plasmid pBR322.^ Screening of duplicate colony filters with the ('32)P-labeled-17-mer and a ('32)P-labeled-17-mer-primed cDNA probe indicated 25 possible glucagon clones from 3100 colonies screened. Restriction mapping of 6 of these clones suggested that they represented a single mRNA species. Primary sequence analysis of one clone containing a 1200 base pair DNA insert revealed that it contained essentially a full-length copy of glucagon cDNA.^ Analaysis of the cDNA suggested that it encoded an initial translation product of 180 amino acids with an M(,r) = 21 kd. The first initiation codon (ATG, methionine) followed by the longest open reading frame of 540 nucleotides was preceded by a 5'-untranslated region of 90 nucleotides, and was followed by a longer 3'-untranslated region of 471 nucleotides, resulting in a total of 1101 nucleotides. . . . (Author's abstract exceeds stipulated maximum length. Discontinued here with permission of author.) UMI ^

GENETIC ORGANIZATION OF THE STRUCTURAL GENES FOR NITRATE REDUCTASE (TN5)

Nitrate reductase in Escherichia coli is a membrane-bound anaerobic enzyme that is repressed by oxygen and induced by nitrate. The genetic organization of the structural genes for the two larger subunits of nitrate reductase ((alpha) and (beta)) was determined by immunoprecipitation analysis of the formation of these proteins in nitrate reductase-deficient mutants resulting from transposon Tn5 mutagenesis. The results suggested that the genes encoding the (alpha) and (beta) subunits (narG and H) were arranged in an operon with transcription in the direction promoter(--->)(alpha)(--->)(beta). Segments of the chromosome containing the Tn5 inserts from several of the mutants were cloned into plasmid pBR322 and the positions of the transposons determined by restriction mapping. The Tn5 insertion sites were localized on two contiguous EcoRI fragments spanning about 6.6 kilobases of DNA. The narI gene (proposed to encode the (gamma) subunit) was positioned immediately downstream from the (beta)-gene (narH) by Southern analysis of Tn10 insertions into the narI locus. A Tn10 insertion into the narK locus, proposed to encode a nitrate-sensitive repressor of other anaerobic enzymes, was located about 1.5 kilobases upstream from the narGHI operon promoter. The narL locus, proposed to encode a nitrate-sensitive positive regulator of the narGHI operon and known to be genetically linked to the other nar genes, was demonstrated to lie outside a 19.3-kilobase region of the chromosome which encompasses the other nar genes. The physical limit of the narGHI promoter was defined by studying the effect of Tn5 insertions into a hybrid plasmid containing the functional operon. The points of origin of the coding regions for the (alpha) and (beta) genes were deduced by alignment of the chromosomal map of Tn5 insertion sites with the sizes of (alpha) and (beta) subunit fragments produced by plasmids carrying these Tn5 inserts in the nar operon. The coding region for the (alpha) subunit (143,000 daltons) begins about 250 nucleotides downstream from the deduced limit of the promoter region and includes about 4.0 kilobases of DNA; the region encoding (beta) (60,000 daltons) lies immediately downstream from the (alpha)-gene and is approximately 1.6 kilobases in length. The adjacent region encoding the (gamma) subunit (19,000 daltons) is approximately 0.5 kilobase in length. ^

Cloning, characterization and regulation of the soluble guanylyl cyclase alpha1 and beta1 genes

Cell signaling by nitric oxide (NO) through soluble guanylyl cyclase (sGC) and cGMP production regulates physiological responses such as smooth muscle relaxation, neurotransmission, and cell growth and differentiation. Although the NO receptor, sGC, has been studied extensively at the protein level, information on regulation of the sGC genes remains elusive. In order to understand the molecular mechanisms involved at the level of gene expression, cDNA and genomic fragments of the murine sGCα1 subunit gene were obtained through library screenings. Using the acquired clones, the sGCα 1 gene structure was determined following primer extension, 3 ′RACE and intron/exon boundary analyses. The basal activity of several 5′-flanking regions (putative promoter regions) for both the α1 and β1 sGC subunits were determined following their transfection into mouse N1E-115 neuroblastoma and rat RENE1Δ14 uterine epithelial cells using a luciferase reporter plasmid. Using the sGC sequences, real-time RT-PCR assays were designed to measure mRNA levels of the sGC α1 and β1 genes in rat, mouse and human. Subsequent studies found that uterine sGC mRNA and protein levels decreased rapidly in response to 17β-estradiol (estrogen) in an in vivo rat model. As early as 1 hour following treatment, mRNA levels of both sGC mRNAs decreased, and reached their lowest level of expression after 3 hours. This in vivo response was completely blocked by the pure estrogen receptor antagonist, ICI 182,780, was not seen in several other tissues examined, did not occur in response to other steroid hormones, and was due to a post-transcriptional mechanism. Additional studies ex vivo and in various cell culture models suggested that the estrogen-mediated decreased sGC mRNA expression did not require signals from other tissues, but may require cell communication or paracrine factors between different cell types within the uterus. Using chemical inhibitors and molecular targeting in other related studies, it was revealed that c-Jun-N-terminal kinase (JNK) signaling was responsible for decreased sGC mRNA expression in rat PC12 and RFL-6 cells, two models previously determined to exhibit rapid decreased sGC mRNA expression in response to different stimuli. To further investigate the post-transcriptional gene regulation, the full length sGCα1 3′-untranslated region (3′UTR) was cloned from rat uterine tissue and ligated downstream of the rabbit β-globin gene and expressed as a chimeric mRNA in the rat PC12 and RFL-6 cell models. Expression studies with the chimeric mRNA showed that the sGCα 1 3′UTR was not sufficient to mediate the post-transcriptional regulation of its mRNA by JNK or cAMP signaling in PC12 and RFL-6 cells. This study has provided numerous valuable tools for future studies involving the molecular regulation of the sGC genes. Importantly, the present results identified a novel paradigm and a previously unknown signaling pathway for sGC mRNA regulation that could potentially be exploited to treat diseases such as uterine cancers, neuronal disorders, hypertension or various inflammatory conditions. ^

Evaluation of PRESAGE® dosimeters for brachytherapy sources and the 3D dosimetry and characterization of the new AgX100 125I seed model

With continuous new improvements in brachytherapy source designs and techniques, method of 3D dosimetry for treatment dose verifications would better ensure accurate patient radiotherapy treatment. This study was aimed to first evaluate the 3D dose distributions of the low-dose rate (LDR) Amersham 6711 OncoseedTM using PRESAGE® dosimeters to establish PRESAGE® as a suitable brachytherapy dosimeter. The new AgX100 125I seed model (Theragenics Corporation) was then characterized using PRESAGE® following the TG-43 protocol. PRESAGE® dosimeters are solid, polyurethane-based, 3D dosimeters doped with radiochromic leuco dyes that produce a linear optical density response to radiation dose. For this project, the radiochromic response in PRESAGE® was captured using optical-CT scanning (632 nm) and the final 3D dose matrix was reconstructed using the MATLAB software. An Amersham 6711 seed with an air-kerma strength of approximately 9 U was used to irradiate two dosimeters to 2 Gy and 11 Gy at 1 cm to evaluate dose rates in the r=1 cm to r=5 cm region. The dosimetry parameters were compared to the values published in the updated AAPM Report No. 51 (TG-43U1). An AgX100 seed with an air-kerma strength of about 6 U was used to irradiate two dosimeters to 3.6 Gy and 12.5 Gy at 1 cm. The dosimetry parameters for the AgX100 were compared to the values measured from previous Monte-Carlo and experimental studies. In general, the measured dose rate constant, anisotropy function, and radial dose function for the Amersham 6711 showed agreements better than 5% compared to consensus values in the r=1 to r=3 cm region. The dose rates and radial dose functions measured for the AgX100 agreed with the MCNPX and TLD-measured values within 3% in the r=1 to r=3 cm region. The measured anisotropy function in PRESAGE® showed relative differences of up to 9% with the MCNPX calculated values. It was determined that post-irradiation optical density change over several days was non-linear in different dose regions, and therefore the dose values in the r=4 to r=5 cm regions had higher uncertainty due to this effect. This study demonstrated that within the radial distance of 3 cm, brachytherapy dosimetry in PRESAGE® can be accurate within 5% as long as irradiation times are within 48 hours.

Regulations of SOX9 activity during chondrogenesis

Sox9 is a transcription factor required for chondrocyte differentiation and cartilage formation. In an effort to identify SOX9 interacting protein(s), we screened a chondrocyte cDNA library with a modified yeast two-hybrid method, Son of Sevenless (SOS) recruitment system (SRS). The catalytic subunit of cyclic AMP-dependent protein kinase A (PKA-Cα) and a new long form of c-Maf transcription factor (Lc-Maf) were found to interact specifically with SOX9. We showed here that two PKA phosphorylation consensus sites of SOX9 could be phosphorylated by PKA in vitro as well as in vivo. PKA phosphorylation of SOX9 increases its DNA binding and transcriptional activities on a Col2a1 chondrocyte-specific enhancer. Mutations of these two PKA phosphorylation sites markedly decreased the activation of SOX9 by PKA. ^ To test whether parathyroid hormone-related peptide (PTHrP) signaling results in SOX9 phosphorylation, we generated a phosphospecific antibody that specifically recognizes SOX9 that is phosphorylated at serine 181 (S 181) one of the two consensus PKA phosphorylation sites. Addition of PTHrP to COS7 cells cotransfected with SOX9 and PTH/PTHrP receptor strongly increased phosphorylation of SOX9 at S181; this phosphorylation was blocked by a PKA-specific inhibitor. In similar experiments we showed that PTHrP increased the activity of a SOX9-dependent Col2a1 enhancer. This increase in activity was abolished when a SOX9 mutant was used containing serine-to-alanine substitution in the two consensus PKA phosphorylation sites of SOX9. Using our phosphospecific SOX9 antibody we showed by immunohistochemistry of mouse embryos that Sox9 phosphorylated at S181 was localized almost exclusively in the pre-hypertrophic zone of the growth plate, an area corresponding to the major site of expression of PTH/PTHrP receptor. In contrast, no phosphorylation of Sox9 at S181 was detected in growth plates of PTH/PTHrP receptor null mutant mice. Sox9, regardless of phosphorylation state, was present in all chondrocytes of both genotypes except in hypertrophic chondrocytes. Thus, Sox9 is a target of PTHrP signaling and the PTHrP-dependent phosphorylation of SOX9 enhances its transcriptional activity. ^ In order to investigate the in vivo function of Sox9 phosphorylation by PKA, we are generating a mouse model of mutant Sox9 harboring point mutations in two PKA phosphorylation sites. Preliminary results indicated that heterozygous mice containing half amount of mutant Sox9 that can not be phosphorylated by PKA have normal skeletal phenotype and homozygous mice are being generated. ^ Lc-Maf encodes an extra ten amino acids at the carboxyl terminus of c-Maf and contains a completely different 3′ untranslated region. The interaction between SOX9 and Lc-Maf was further confirmed by co-immunoprecipitation and GST-pull down assays, which mapped the interacting domains of SOX9 to HMG DNA binding domain and that of Lc-Maf to basic leusine zipper motif. In situ hybridizations showed that RNA of Lc-Maf coexpressed with those of Sox9 and Col2a1 in areas of mesenchymal condensation during the early stages of mouse embryo development. A DNA binding site of Lc-Maf was identified at the 5′ part of a 48-bp Col2a1 enhancer element near the HMG binding site of SOX9. Lc-Maf and SOX9 synergistically activated a luciferase reporter plasmid containing a Col2al enhancer and increased the transcription of endogenous Col2a1 gene. In summary, Lc-Maf is the first identified SOX9-interating protein during chondrogenesis and may be an important activator of Col2a1 gene. ^