Demographic history and linkage disequilibrium in human populations.

In the human genome, linkage disequilibrium (LD)--the non-random association of alleles at chromosomal loci--has been studied mainly in regions surrounding disease genes on affected chromosomes. Consequently, little information is available on the distribution of LD across anonymous genomic regions in the general population. However, demographic history is expected to influence the extent of overall LD across the genome, so a population that has been of constant size will display higher levels of LD than a population that has expanded. In support of this, the extent of LD between anonymous loci on chromosome 4 in chimpanzees (as a model of a population of constant size) has been compared to that in Finns (as a model of an expanded population; refs 8,9) and found to exhibit more LD than in the latter population. In Europe, studies of mitochondrial (mt) DNA sequences have suggested that most populations have experienced expansion, whereas the Saami in northern Fenno-Scandinavia have been of constant size (Table 1). Thus, in northern Europe, populations with radically different demographic histories live in close geographic proximity to each other. We studied the allelic associations between anonymous microsatellite loci on the X chromosome in the Saami and neighbouring populations and found dramatically higher levels of LD in the Saami than in other populations in the region. This indicates that whereas recently expanded populations, such as the Finns, are well suited to map single disease genes affected by recent mutations, populations that have been of constant size, such as the Saami, may be much better suited to map genes for complex traits that are caused by older mutations.

The genetical archaeology of the human genome.

Palaentology and archaeology are disciplines that traditionally deal with the reconstruction of human origins and history. Recently, however, molecular genetics has come to make increasing contributions to this area. In particular, several data sets indicate that variation of the human gene pool originated in Africa within the last 200,000 years. Furthermore, the study of DNA sequences allows the detection of expansions in population size. Here we briefly summarize and exemplify how DNA sequences can be used to reconstruct the history of populations.

Great ape DNA sequences reveal a reduced diversity and an expansion in humans.

The extent of DNA sequence variation of chimpanzees is several-fold greater than that of humans. It is unclear, however, if humans or chimpanzees are exceptional among primates in having low and high amounts of DNA sequence diversity, respectively. To address this, we have determined approximately 10,000 bp of noncoding DNA sequences at Xq13.3 (which has been extensively studied in both humans and chimpanzees) from 10 western lowland gorillas (Gorilla gorilla gorilla) and 1 mountain gorilla (Gorilla gorilla beringei; that is, from 2 of the 3 currently recognized gorilla subspecies), as well as 8 Bornean (Pongo pygmaeus pygmaeus) and 6 Sumatran (Pongo pygmaeus abelii) orang-utans, representing both currently recognized orang-utan subspecies. We show that humans differ from the great apes in having a low level of genetic variation and a signal of population expansion.

DNA sequence variation in a non-coding region of low recombination on the human X chromosome.

DNA sequence variation has become a major source of insight regarding the origin and history of our species as well as an important tool for the identification of allelic variants associated with disease. Comparative sequencing of DNA has to date focused mainly on mitochondrial (mt) DNA, which due to its apparent lack of recombination and high evolutionary rate lends itself well to the study of human evolution. These advantages also entail limitations. For example, the high mutation rate of mtDNA results in multiple substitutions that make phylogenetic analysis difficult and, because mtDNA is maternally inherited, it reflects only the history of females. For the history of males, the non-recombining part of the paternally inherited Y chromosome can be studied. The extent of variation on the Y chromosome is so low that variation at particular sites known to be polymorphic rather than entire sequences are typically determined. It is currently unclear how some forms of analysis (such as the coalescent) should be applied to such data. Furthermore, the lack of recombination means that selection at any locus affects all 59 Mb of DNA. To gauge the extent and pattern of point substitutional variation in non-coding parts of the human genome, we have sequenced 10 kb of non-coding DNA in a region of low recombination at Xq13.3. Analysis of this sequence in 69 individuals representing all major linguistic groups reveals the highest overall diversity in Africa, whereas deep divergences also exist in Asia. The time elapsed since the most recent common ancestor (MRCA) is 535,000+/-119,000 years. We expect this type of nuclear locus to provide more answers about the genetic origin and history of humans.

Minisatellite diversity supports a recent African origin for modern humans.

In a study of human diversity at a highly variable locus, we have mapped the internal structures of tandem-repetitive alleles from different populations at the minisatellite MS205 (D16S309). The results give an unusually detailed view of the different allelic structures represented on modern human chromosomes, and of the ancestral relationships between them. There was a clear difference in allelic diversity between African and non-African populations. A restricted set of allele families was found in non-African populations, and formed a subset of the much greater diversity seen on African chromosomes. The data strongly support a recent African origin for modern human diversity at this locus.

Targeted resequencing of a genomic region influencing tameness and aggression reveals multiple signals of positive selection.

The identification of the causative genetic variants in quantitative trait loci (QTL) influencing phenotypic traits is challenging, especially in crosses between outbred strains. We have previously identified several QTL influencing tameness and aggression in a cross between two lines of wild-derived, outbred rats (Rattus norvegicus) selected for their behavior towards humans. Here, we use targeted sequence capture and massively parallel sequencing of all genes in the strongest QTL in the founder animals of the cross. We identify many novel sequence variants, several of which are potentially functionally relevant. The QTL contains several regions where either the tame or the aggressive founders contain no sequence variation, and two regions where alternative haplotypes are fixed between the founders. A re-analysis of the QTL signal showed that the causative site is likely to be fixed among the tame founder animals, but that several causative alleles may segregate among the aggressive founder animals. Using a formal test for the detection of positive selection, we find 10 putative positively selected regions, some of which are close to genes known to influence behavior. Together, these results show that the QTL is probably not caused by a single selected site, but may instead represent the joint effects of several sites that were targets of polygenic selection.

DNA sequence and comparative analysis of chimpanzee chromosome 22.

Human-chimpanzee comparative genome research is essential for narrowing down genetic changes involved in the acquisition of unique human features, such as highly developed cognitive functions, bipedalism or the use of complex language. Here, we report the high-quality DNA sequence of 33.3 megabases of chimpanzee chromosome 22. By comparing the whole sequence with the human counterpart, chromosome 21, we found that 1.44% of the chromosome consists of single-base substitutions in addition to nearly 68,000 insertions or deletions. These differences are sufficient to generate changes in most of the proteins. Indeed, 83% of the 231 coding sequences, including functionally important genes, show differences at the amino acid sequence level. Furthermore, we demonstrate different expansion of particular subfamilies of retrotransposons between the lineages, suggesting different impacts of retrotranspositions on human and chimpanzee evolution. The genomic changes after speciation and their biological consequences seem more complex than originally hypothesized.

Human evolution.

The origin, history, and singularity of our species has fascinated storytellers, philosophers and scientists throughout, and doubtless before, recorded history. Anthropology, the modern-era discipline that deals with these issues, is a notoriously contentious field, perhaps because the topic at hand - the nature of our own species - is one that is difficult or impossible to approach in an unbiased way. Recently, molecular genetics has increasingly contributed to this field. Here, I briefly discuss three areas where I believe molecular studies are likely to be of decisive importance in the future. These concern the questions of where and when our species originated, what the genetic background for characters that differ between us and apes is, and how the phenotypic traits that vary among human groups have evolved.

Molecular analyses of oral polio vaccine samples.

It has been suggested that the human immunodeficiency virus (HIV), and thus the acquired immunodeficiency syndrome (AIDS) it causes, was inadvertently introduced to humans by the use of an oral polio vaccine (OPV) during a vaccination campaign launched by the Wistar Institute, Philadelphia, PA, USA, in the Belgian Congo in 1958 and 1959. The "OPV/AIDS hypothesis" suggests that the OPV used in this campaign was produced in chimpanzee kidney epithelial cell cultures rather than in monkey kidney cell cultures, as stated by H. Koprowski and co-workers, who produced the OPV. If chimpanzee cells were indeed used, this would lend support to the OPV/AIDS hypothesis, since chimpanzees harbor a simian immunodeficiency virus, widely accepted to be the origin of HIV-1. We analyzed several early OPV pools and found no evidence for the presence of chimpanzee DNA; by contrast, monkey DNA is present.

Extensive nuclear DNA sequence diversity among chimpanzees.

Although data on nucleotide sequence variation in the human nuclear genome have begun to accumulate, little is known about genomic diversity in chimpanzees (Pan troglodytes) and bonobos (Pan paniscus). A 10,154-base pair sequence on the chimpanzee X chromosome is reported, representing all major subspecies and bonobos. Comparison to humans shows the diversity of the chimpanzee sequences to be almost four times as high and the age of the most recent common ancestor three times as great as the corresponding values of humans. Phylogenetic analyses show the sequences from the different chimpanzee subspecies to be intermixed and the distance between some chimpanzee sequences to be greater than the distance between them and the bonobo sequences.

Amino acid racemization and the preservation of ancient DNA.

The extent of racemization of aspartic acid, alanine, and leucine provides criteria for assessing whether ancient tissue samples contain endogenous DNA. In samples in which the D/L ratio of aspartic acid exceeds 0.08, ancient DNA sequences could not be retrieved. Paleontological finds from which DNA sequences purportedly millions of years old have been reported show extensive racemization, and the amino acids present are mainly contaminates. An exception is the amino acids in some insects preserved in amber.

Molecular coproscopy: dung and diet of the extinct ground sloth Nothrotheriops shastensis.

DNA from excrements can be amplified by means of the polymerase chain reaction. However, this has not been possible with ancient feces. Cross-links between reducing sugars and amino groups were shown to exist in a Pleistocene coprolite from Gypsum Cave, Nevada. A chemical agent, N-phenacylthiazolium bromide, that cleaves such cross-links made it possible to amplify DNA sequences. Analyses of these DNA sequences showed that the coprolite is derived from an extinct sloth, presumably the Shasta ground sloth Nothrotheriops shastensis. Plant DNA sequences from seven groups of plants were identified in the coprolite. The plant assemblage that formed part of the sloth's diet exists today at elevations about 800 meters higher than the cave.

Evidence for import of a lysyl-tRNA into marsupial mitochondria.

The mitochondrial tRNA gene for lysine was analyzed in 11 different marsupial mammals. Whereas its location is conserved when compared with other vertebrate mitochondrial genomes, its primary sequence and inferred secondary structure are highly unusual and variable. For example, eight species lack the expected anticodon. Because the corresponding transcripts are not altered by any RNA-editing mechanism, the lysyl-tRNA gene seems to represent a mitochondrial pseudogene. Purification of marsupial mitochondria and in vitro aminoacylation of isolated tRNAs with lysine, followed by analysis of aminoacylated tRNAs, show that a nuclear-encoded tRNA(Lys) is associated with marsupial mitochondria. We conclude that a functional tRNA(Lys) encoded in the nuclear genome is imported into mitochondria in marsupials. Thus, tRNA import is not restricted to plant, yeast, and protozoan mitochondria but also occurs also in mammals.

DNA sequences from multiple amplifications reveal artifacts induced by cytosine deamination in ancient DNA.

We show that DNA molecules amplified by PCR from DNA extracted from animal bones and teeth that vary in age between 25 000 and over 50 000 years carry C-->T and G-->A substitutions. These substitutions can reach high proportions among the molecules amplified and are due to the occurrence of modified deoxycytidine residues in the template DNA. If the template DNA is treated with uracil N-glycosylase, these substitutions are dramatically reduced. They are thus likely to result from deamination of deoxycytidine residues. In addition, 'jumping PCR', i.e. the occurrence of template switching during PCR, may contribute to these substitutions. When DNA sequences are amplified from ancient DNA extracts where few template molecules initiate the PCR, precautions such as DNA sequence determination of multiple clones derived from more than one independent amplification are necessary in order to reduce the risk of determination of incorrect DNA sequences. When such precautionary measures are taken, errors induced by damage to the DNA template are unlikely to be more frequent than approximately 0.1% even under the unlikely scenario where each amplification starts from a single template molecule.

Adenovirus proteins and MHC expression.

Adenoviruses are able to specifically down-regulate the cell surface expression of MHC class I antigens. Most viral serotypes achieve these ends by synthesizing a protein that binds to class I antigens in the endoplasmic reticulum (ER) and impedes the transport of these molecules to the cell surface. However, viruses belonging to the highly oncogenic subgenus A do not affect the class I antigen expression during acute infection. Instead, they are distinct from other adenoviruses in that they specifically down-regulate the level of mRNAs, encoding MHC class I antigens, in virally transformed cells. The virus-induced reduction of class I antigen expression drastically diminishes the ability of CTLs to recognize cells infected or transformed by adenovirus. A number of issues concerning these viral mechanisms for class I antigen modulation need to be addressed. The molecular mechanism by which the E1A gene product of subgenus A viruses diminishes class I mRNA levels has not been elucidated. Also, the details of the interaction between the E19 protein and class I molecules should be studied, preferably by X-ray crystallography of the complexes. This would clarify the role of the antigen-binding site as well as other portions of the class I molecule in the binding to the E19 protein. Of general importance for our understanding of the sorting and intracellular transport of proteins is the exact delimitation of the signal for ER localization, which is present in the COOH-terminus of the E19 protein. The putative interaction of this peptide sequence with components of the ER membrane should also be studied. Finally, the study of the pathophysiological role of the MHC class I down-regulation will undoubtedly yield new insights into how the immune system combats virally infected and transformed cells.

Polyadenylation creates the discriminator nucleotide of chicken mitochondrial tRNA(Tyr).

In the chicken mitochondrial genome, the gene for tRNA(Tyr) overlaps by one nucleotide with the downstream tRNA(Cys) gene, which is located on the same strand. The overlapping nucleotide, a guanosine residue, thus encodes both the discriminator base of the tRNA(Tyr) and the 5'base of the tRNA(Cys). When cDNA clones of circularized forms of the tRNA(Tyr) are analyzed, the discriminator nucleotide is an adenosine residue rather than the genomically encoded guanosine. Thus, the tRNA(Tyr) is subjected to an RNA editing activity similar to that shown to exist in the mitochondria of two other animal species. Interestingly, some cDNA clones have several adenosine residues at their 3'-ends instead of the expected CCA-sequence. Furthermore, a review of sequence data from animal mitochondrial genomes suggests that only tRNAs whose discriminator bases are adenosines tend to have genes that overlap with downstream genes. Thus, polyadenylation seems to be a major component of the RNA editing machinery that affects overlapping genes in animal mitochondria.

The mitochondrial genome of the hemichordate Balanoglossus carnosus and the evolution of deuterostome mitochondria.

The complete nucleotide sequence of the mitochondrial genome of the hemichordate Balanoglossus carnosus (acorn worm) was determined. The arrangement of the genes encoding 13 protein, 22 tRNA, and 2 rRNA genes is essentially the same as in vertebrates, indicating that the vertebrate and hemichordate mitochondrial gene arrangement is close to that of their common ancestor, and, thus, that it has been conserved for more than 600 million years, whereas that of echinoderms has been rearranged extensively. The genetic code of hemichordate mitochondria is similar to that of echinoderms in that ATA encodes isoleucine and AGA serine, whereas the codons AAA and AGG, whose amino acid assignments also differ between echinoderms and vertebrates, are absent from the B. carnosus mitochondrial genome. There are three noncoding regions of length 277, 41, and 32 bp: the larger one is likely to be equivalent to the control region of other deuterostomes, while the two others may contain transcriptional promoters for genes encoded on the minor coding strand. Phylogenetic trees estimated from the inferred protein sequences indicate that hemichordates are a sister group of echinoderms.

The marsupial mitochondrial genome and the evolution of placental mammals.

The entire nucleotide sequence of the mitochondrial genome of the American opossum, Didelphis virginiana, was determined. Two major features distinguish this genome from those of other mammals. First, five tRNA genes around the origin of light strand replication are rearranged. Second, the anticodon of tRNA(Asp) is posttranscriptionally changed by an RNA editing process such that its coding capacity is altered. When the complete protein-coding region of the mitochondrial genome is used as an outgroup for placental mammals it can be shown that rodents represent an earlier branch among placental mammals than primates and artiodactyls and that artiodactyls share a common ancestor with carnivores. The overall rates of evolution of most of the mitochondrial genome of placentals are clock-like. Furthermore, the data indicate that the lineages leading to the mouse and rat may have diverged from each other as much as 35 million years ago.