Theoretical and experimental studies of linkage disequilibrium in human populations

Linkage disequilibrium (LD) is defined as the nonrandom association of alleles at two or more loci in a population and may be a useful tool in a diverse array of applications including disease gene mapping, elucidating the demographic history of populations, and testing hypotheses of human evolution. However, the successful application of LD-based approaches to pertinent genetic questions is hampered by a lack of understanding about the forces that mediate the genome-wide distribution of LD within and between human populations. Delineating the genomic patterns of LD is a complex task that will require interdisciplinary research that transcends traditional scientific boundaries. The research presented in this dissertation is predicated upon the need for interdisciplinary studies and both theoretical and experimental projects were pursued. In the theoretical studies, I have investigated the effect of genotyping errors and SNP identification strategies on estimates of LD. The primary importance of these two chapters is that they provide important insights and guidance for the design of future empirical LD studies. Furthermore, I analyzed the allele frequency distribution of 26,530 single nucleotide polymorphisms (SNPs) in three populations and generated the first-generation natural selection map of the human genome, which will be an important resource for explaining and understanding genomic patterns of LD. Finally, in the experimental study, I describe a novel and simple, low-cost, and high-throughput SNP genotyping method. The theoretical analyses and experimental tools developed in this dissertation will facilitate a more complete understanding of patterns of LD in human populations. ^

Chimpanzee: A predictive model for the human cytochrome P450 3A subfamily

The human cytochrome P450 3A (CYP3A) subfamily is responsible for most of the metabolism of therapeutic drugs; however, an adequate in vivo model has yet to be discovered. This study begins with an investigation of a controversial topic surrounding the human CYP3As--estrogen regulation. A novel approach to this topic was used by defining expression in the estrogen-responsive endometrium. This study shows that estrogen down-regulates CYP3A4 expression in the endometrium. On the other hand, analogous studies showed an increase in CYP3A expression as age increases in liver tissue. Following the discussion of estrogen regulation, is an investigation of the cross-species relationships among all of the CYP3As was completed. The study compares isoforms from piscines, avians, rodents, canines, ovines, bovines, and primates. Using the traditional phylogenetic analyses and employing a novel approach using exon and intron lengths, the results show that only another primate could be the best animal model for analysis of the regulation of the expression of the human CYP3As. This analysis also demonstrated that the chimpanzee seems to be the best available human model. Moreover, the study showed the presence and similarities of one additional isoform in the chimpanzee genome that is absent in humans. Based on these results, initial characterization of the chimpanzee CYP3A subfamily was begun. While the human genome contains four isoforms--CYP3A4, CYP3A5, CYP3A7, and CYP3A43--the chimpanzee genome has five, the four previously mentioned and CYP3A67. Both species express CYP3A4, CYP3A5, and CYP3A43, but humans express CYP3A7 while chimpanzees express CYP3A67. In humans, CYP3A4 is expressed at higher levels than the other isoforms, but some chimpanzee individuals express CYP3A67 at higher levels than CYP3A4. Such a difference is expected to alter significantly the total CYP3A metabolism. On the other hand, any study considering individual isoforms would still constitute a valid method of study for the human CYP3A4, CYP3A5, and CYP3A43 isoforms. ^

Functional copies of a human gene can be directly isolated by transformation-associated recombination cloning with a small 3′ end target sequence

Unique, small sequences (sequence tag sites) have been identified at the 3′ ends of most human genes that serve as landmarks in genome mapping. We investigated whether a single copy gene could be isolated directly from total human DNA by transformation-associated recombination (TAR) cloning in yeast using a short, 3′ unique target. A TAR cloning vector was constructed that, when linearized, contained a small amount (381 bp) of 3′ hypoxanthine phosphoribosyltransferase (HPRT) sequence at one end and an 189-bp Alu repeat at the other end. Transformation with this vector along with human DNA led to selective isolations of the entire HPRT gene as yeast artificial chromosomes (YACs) that extended from the 3′ end sequence to various Alu positions as much as 600 kb upstream. These YACs were retrofitted with a NeoR and a bacterial artificial chromosome (BAC) sequence to transfer the YACs to bacteria and subsequently the BACs to mouse cells by using a Neo selection. Most of the HPRT isolates were functional, demonstrating that TAR cloning retains the functional integrity of the isolated material. Thus, this modified version of TAR cloning, which we refer to as radial TAR cloning, can be used to isolate large segments of the human genome accurately and directly with only a small amount of sequence information.

Inactivation of the survival motor neuron gene, a candidate gene for human spinal muscular atrophy, leads to massive cell death in early mouse embryos

Proximal spinal muscular atrophy is an autosomal recessive human disease of spinal motor neurons leading to muscular weakness with onset predominantly in infancy and childhood. With an estimated heterozygote frequency of 1/40 it is the most common monogenic disorder lethal to infants; milder forms represent the second most common pediatric neuromuscular disorder. Two candidate genes—survival motor neuron (SMN) and neuronal apoptosis inhibitory protein have been identified on chromosome 5q13 by positional cloning. However, the functional impact of these genes and the mechanism leading to a degeneration of motor neurons remain to be defined. To analyze the role of the SMN gene product in vivo we generated SMN-deficient mice. In contrast to the human genome, which contains two copies, the mouse genome contains only one SMN gene. Mice with homozygous SMN disruption display massive cell death during early embryonic development, indicating that the SMN gene product is necessary for cellular survival and function.

Complete sequencing of the Fugu WAGR region from WT1 to PAX6: Dramatic compaction and conservation of synteny with human chromosome 11p13

The pufferfish Fugu rubripes has a genome ≈7.5 times smaller than that of mammals but with a similar number of genes. Although conserved synteny has been demonstrated between pufferfish and mammals across some regions of the genome, there is some controversy as to what extent Fugu will be a useful model for the human genome, e.g., [Gilley, J., Armes, N. & Fried, M. (1997) Nature (London) 385, 305–306]. We report extensive conservation of synteny between a 1.5-Mb region of human chromosome 11 and <100 kb of the Fugu genome in three overlapping cosmids. Our findings support the idea that the majority of DNA in the region of human chromosome 11p13 is intergenic. Comparative analysis of three unrelated genes with quite different roles, WT1, RCN1, and PAX6, has revealed differences in their structural evolution. Whereas the human WT1 gene can generate 16 protein isoforms via a combination of alternative splicing, RNA editing, and alternative start site usage, our data predict that Fugu WT1 is capable of generating only two isoforms. This raises the question of the extent to which the evolution of WT1 isoforms is related to the evolution of the mammalian genitourinary system. In addition, this region of the Fugu genome shows a much greater overall compaction than usual but with significant noncoding homology observed at the PAX6 locus, implying that comparative genomics has identified regulatory elements associated with this gene.

An ancient family of human endogenous retroviruses encodes a functional homolog of the HIV-1 Rev protein

The human endogenous retrovirus K (HERV-K) family of endogenous retroviruses consists of ≈50 proviral copies per haploid human genome. Herein, the HERV-Ks are shown to encode a sequence-specific nuclear RNA export factor, termed K-Rev, that is functionally analogous to the HIV-1 Rev protein. Like HIV-1 Rev, K-Rev binds to both the Crm1 nuclear export factor and to a cis-acting viral RNA target to activate nuclear export of unspliced RNAs. Surprisingly, this HERV-K RNA sequence, which is encoded within the HERV-K long terminal repeat, is also recognized by HIV-1 Rev. These data provide surprising evidence for an evolutionary link between HIV-1 and a group of endogenous retroviruses that first entered the human genome ≈30 million years ago.

Whole genome amplification using a degenerate oligonucleotide primer allows hundreds of genotypes to be performed on less than one nanogram of genomic DNA

Genetic analysis of limiting quantities of genomic DNA play an important role in DNA forensics, paleoarcheology, genetic disease diagnosis, genetic linkage analysis, and genetic diversity studies. We have tested the ability of degenerate oligonucleotide primed polymerase chain reaction (DOP-PCR) to amplify picogram quantities of human genomic DNA for the purpose of increasing the amount of template for genotyping with microsatellite repeat markers. DNA was uniformly amplified at a large number of typable loci throughout the human genome with starting template DNAs from as little as 15 pg to as much as 400 ng. A much greater-fold enrichment was seen for the smaller genomic DOP-PCRs. All markers tested were amplified from starting genomic DNAs in the range of 0.6–40 ng with amplifications of 200- to 600-fold. The DOP-PCR-amplified genomic DNA was an excellent and reliable template for genotyping with microsatellites, which give distinct bands with no increase in stutter artifact on di-, tri-, and tetranucleotide repeats. There appears to be equal amplification of genomic DNA from 55 of 55 tested discrete microsatellites implying near complete coverage of the human genome. Thus, DOP-PCR appears to allow unbiased, hundreds-fold whole genome amplification of human genomic DNA for genotypic analysis.

Human trophoblast and choriocarcinoma expression of the growth factor pleiotrophin attributable to germ-line insertion of an endogenous retrovirus

Retroviral elements are found in abundance throughout the human genome but only rarely have alterations of endogenous genes by retroviral insertions been described. Herein we report that a human endogenous retrovirus (HERV) type C is inserted in the human growth factor gene pleiotrophin (PTN) between the 5′ untranslated and the coding region. This insert in the human genome expands the region relative to the murine gene. Studies with promoter-reporter constructs show that the HERV insert in the human PTN gene generates an additional promoter with trophoblast-specific activity. Due to this promoter function, fusion transcripts between HERV and the open reading frame of PTN (HERV-PTN) were detected in all normal human trophoblast cell cultures as early as 9 weeks after gestation (n = 7) and in all term placenta tissues (n = 5) but not in other normal adult tissues. Furthermore, only trophoblast-derived choriocarcinoma cell lines expressed HERV-PTN mRNA whereas tumor cell lines derived from the embryoblast (teratocarcinoma) or from other lineages failed to do so. We investigated the significance of HERV-PTN mRNA in a choriocarcinoma model by targeting this transcript with ribozymes and found that the depletion of HERV-PTN mRNA prevents human choriocarcinoma growth, invasion, and angiogenesis in mice. This suggests that the tissue-specific expression of PTN due to the HERV insertion in the human genome supports the highly aggressive growth of human choriocarcinoma and possibly of the human trophoblast.

Yeast and human genes that affect the Escherichia coli SOS response

The sequencing of the human genome has led to the identification of many genes whose functions remain to be determined. Because of conservation of genetic function, microbial systems have often been used for identification and characterization of human genes. We have investigated the use of the Escherichia coli SOS induction assay as a screen for yeast and human genes that might play a role in DNA metabolism and/or in genome stability. The SOS system has previously been used to analyze bacterial and viral genes that directly modify DNA. An initial screen of meiotically expressed yeast genes revealed several genes associated with chromosome metabolism (e.g., RAD51 and HHT1 as well as others). The SOS induction assay was then extended to the isolation of human genes. Several known human genes involved in DNA metabolism, such as the Ku70 end-binding protein and DNA ligase IV, were identified, as well as a large number of previously unknown genes. Thus, the SOS assay can be used to identify and characterize human genes, many of which may participate in chromosome metabolism.

The Mouse Genome Database (MGD): integration nexus for the laboratory mouse

The Mouse Genome Database (MGD) is the community database resource for the laboratory mouse, a key model organism for interpreting the human genome and for understanding human biology and disease (http://www.informatics.jax.org). MGD provides standard nomenclature and consensus map positions for mouse genes and genetic markers; it provides a curated set of mammalian homology records, user-defined chromosomal maps, experimental data sets and the definitive mouse ‘gene to sequence’ reference set for the research community. The integration and standardization of these data sets facilitates the transition between mouse DNA sequence, gene and phenotype annotations. A recent focus on allele and phenotype representations enhances the ability of MGD to organize and present data for exploring the relationship between genotype and phenotype. This link between the genome and the biology of the mouse is especially important as phenotype information grows from large mutagenesis projects and genotype information grows from large-scale sequencing projects.

Identification of proacrosin binding protein sp32 precursor as a human cancer/testis antigen

Serological expression cloning of antigens eliciting a humoral immune response to a syngeneic mouse sarcoma identified pem (mouse placenta and embryonic expression gene) as a new member of the cancer/testis family. To identify the human homologue of pem, mouse pem sequences and pem-related expressed sequence tags from human testis were used as PCR primers for amplification using human testis cDNA. However, rather than pem, another gene, designated OY-TES-1, was isolated and found to be the human homologue of proacrosin binding protein sp32 precursor originally identified in mouse, guinea pig, and pig. OY-TES-1 maps to chromosome 12p12-p13 and contains 10 exons. Southern blot analysis suggests the presence of two OY-TES-1-related genes in the human genome. In normal tissues, OY-TES-1 mRNA was expressed only in testis, whereas in malignant tissues, a variable proportion of a wide array of cancers, including bladder, breast, lung, liver, and colon cancers, expressed OY-TES-1. Serological survey of 362 cancer patients with a range of different cancers showed antibody to OY-TES-1 in 25 patients. No OY-TES-1 sera reactivity was found in 20 normal individuals. These findings indicate that OY-TES-1 is an additional member of the cancer/testis family of antigens and that OY-TES-1 is immunogenic in humans.

Importance of anchor genomes for any plant genome project

Progress in agricultural and environmental technologies is hampered by a slower rate of gene discovery in plants than animals. The vast pool of genes in plants, however, will be an important resource for insertion of genes, via biotechnological procedures, into an array of plants, generating unique germ plasms not achievable by conventional breeding. It just became clear that genomes of grasses have evolved in a manner analogous to Lego blocks. Large chromosome segments have been reshuffled and stuffer pieces added between genes. Although some genomes have become very large, the genome with the fewest stuffer pieces, the rice genome, is the Rosetta Stone of all the bigger grass genomes. This means that sequencing the rice genome as anchor genome of the grasses will provide instantaneous access to the same genes in the same relative physical position in other grasses (e.g., corn and wheat), without the need to sequence each of these genomes independently. (i) The sequencing of the entire genome of rice as anchor genome for the grasses will accelerate plant gene discovery in many important crops (e.g., corn, wheat, and rice) by several orders of magnitudes and reduce research and development costs for government and industry at a faster pace. (ii) Costs for sequencing entire genomes have come down significantly. Because of its size, rice is only 12% of the human or the corn genome, and technology improvements by the human genome project are completely transferable, translating in another 50% reduction of the costs. (iii) The physical mapping of the rice genome by a group of Japanese researchers provides a jump start for sequencing the genome and forming an international consortium. Otherwise, other countries would do it alone and own proprietary positions.

All kinesin superfamily protein, KIF, genes in mouse and human

Intracellular transport is essential for morphogenesis and functioning of the cell. The kinesin superfamily proteins (KIFs) have been shown to transport membranous organelles and protein complexes in a microtubule- and ATP-dependent manner. More than 30 KIFs have been reported in mice. However, the nomenclature of KIFs has not been clearly established, resulting in various designations and redundant names for a single KIF. Here, we report the identification and classification of all KIFs in mouse and human genome transcripts. Previously unidentified murine KIFs were found by a PCR-based search. The identification of all KIFs was confirmed by a database search of the total human genome. As a result, there are a total of 45 KIFs. The nomenclature of all KIFs is presented. To understand the function of KIFs in intracellular transport in a single tissue, we focused on the brain. The expression of 38 KIFs was detected in brain tissue by Northern blotting or PCR using cDNA. The brain, mainly composed of highly differentiated and polarized cells such as neurons and glia, requires a highly complex intracellular transport system as indicated by the increased number of KIFs for their sophisticated functions. It is becoming increasingly clear that the cell uses a number of KIFs and tightly controls the direction, destination, and velocity of transportation of various important functional molecules, including mRNA. This report will set the foundation of KIF and intracellular transport research.

Genomic analysis of orthologous mouse and human olfactory receptor loci

Olfactory receptor (OR) genes represent ≈1% of genomic coding sequence in mammals, and these genes are clustered on multiple chromosomes in both the mouse and human genomes. We have taken a comparative genomics approach to identify features that may be involved in the dynamic evolution of this gene family and in the transcriptional control that results in a single OR gene expressed per olfactory neuron. We sequenced ≈350 kb of the murine P2 OR cluster and used synteny, gene linkage, and phylogenetic analysis to identify and sequence ≈111 kb of an orthologous cluster in the human genome. In total, 18 mouse and 8 human OR genes were identified, including 7 orthologs that appear to be functional in both species. Noncoding homology is evident between orthologs and generally is confined within the transcriptional unit. We find no evidence for common regulatory features shared among paralogs, and promoter regions generally do not contain strong promoter motifs. We discuss these observations, as well as OR clustering, in the context of evolutionary expansion and transcriptional regulation of OR repertoires.

Complete structure of the human alpha-albumin gene, a new member of the serum albumin multigene family.

The nucleotide sequence of the human alpha-albumin gene, including 887 bp of the 5'-flanking region and 1311 bp of the 3-flanking region (24,454 in total), was determined from three overlapping lambda phage clones. The sequence spans 22,256 bp from the cap site to the polyadenylylation site, revealing a gene structure of 15 exons separated by 14 introns. The methionine initiation codon ATG is within exon 1; the termination codon TGA is within exon 14. Exon 15 is entirely untranslated and contains the polyadenylylation signal AATAAA. The deduced polypeptide chain is composed of a 21-amino-acid leader peptide, followed by 578 amino acids of the mature protein. There are seven repetitive DNA elements (Alu and Kpn) in the introns and 3-flanking region. The sizes of the 15 alpha-albumin exons match closely those of the albumin, alpha-fetoprotein, and vitamin D-binding protein genes. The exons are symmetrically placed within the three domains of the individual proteins, and they share a characteristic codon splitting pattern that is conserved among members of the gene family. The results provide strong evidence that alpha-albumin belongs to, and most likely completes with, the serum albumin gene family. Based on structural similarity, alpha-albumin appears to be most closely related to alpha-fetoprotein. The complete structure of this family of four tandemly linked genes provides a well-characterized approximately 200 kb locus in the 4q subcentromeric region of the human genome.