Bayesian inference of ancient human demography from individual genome sequences

Whole-genome sequences provide a rich source of information about human evolution. Here we describe an effort to estimate key evolutionary parameters based on the whole-genome sequences of six individuals from diverse human populations. We used a Bayesian, coalescent-based approach to obtain information about ancestral population sizes, divergence times and migration rates from inferred genealogies at many neutrally evolving loci across the genome. We introduce new methods for accommodating gene flow between populations and integrating over possible phasings of diploid genotypes. We also describe a custom pipeline for genotype inference to mitigate biases from heterogeneous sequencing technologies and coverage levels. Our analysis indicates that the San population of southern Africa diverged from other human populations approximately 108-157 thousand years ago, that Eurasians diverged from an ancestral African population 38-64 thousand years ago, and that the effective population size of the ancestors of all modern humans was ～9,000.     

Targeted breakage of a human chromosome mediated by cloned human telomeric DNA.

Novel approaches to the structural and functional analysis of mammalian chromosomes would be possible if the gross structure of the chromosomes in living cells could be engineered. Controlled modifications can be engineered by conventional targeting techniques based on homologous recombination. Large but uncontrolled modifications can be made by the integration of cloned human telomeric DNA. We describe here the combined use of gene targeting and telomere-mediated chromosome breakage to generate a defined truncation of a human chromosome. Telomeric DNA was targeted to the 6-16 gene on the short arm of chromosome 1 in a human cell line. Molecular and cytogenetic analyses showed that, of eight targeted clones that were isolated, one clone had the predicted truncation of chromosome 1.     

Y chromosome sequence variation and the history of human populations

Binary polymorphisms associated with the non-recombining region of the human Y chromosome (NRY) preserve the paternal genetic legacy of our species that has persisted to the present, permitting inference of human evolution, population affinity and demographic history. We used denaturing high-performance liquid chromatography (DHPLC; ref. 2) to identify 160 of the 166 bi-allelic and 1 tri-allelic site that formed a parsimonious genealogy of 116 haplotypes, several of which display distinct population affinities based on the analysis of 1062 globally representative individuals. A minority of contemporary East Africans and Khoisan represent the descendants of the most ancestral patrilineages of anatomically modern humans that left Africa between 35,000 and 89,000 years ago.     

Automated DNA sequencing and analysis of 106 kilobases from human chromosome 19q13.3.

A total of 116,118 basepairs (bp) derived from three cosmids spanning the ERCC1 locus of human chromosome 19q13.3 have been sequenced with automated fluorescence-based sequencers and analysed by polymerase chain reaction amplification and computer methods. The assembled sequence forms two contigs totalling 105,831 bp, which contain a human fosB proto-oncogene, a gene encoding a protein phosphatase, two genes of unknown function and the previously-characterized ERCC1 DNA repair gene. This light band region has a high average density of 1.4 Alu repeats per kilobase. Human chromosome light bands could therefore contain up to 75,000 genes and 1.5 million Alu repeats.     

Isolation of chromosome 21-specific yeast artificial chromosomes from a total human genome library.

A new approach for the isolation of chromosome-specific subsets from a human genomic yeast artificial chromosome (YAC) library is described. It is based on the hybridization with an Alu polymerase chain reaction (PCR) probe. We screened a 1.5 genome equivalent YAC library of megabase insert size with Alu PCR products amplified from hybrid cell lines containing human chromosome 21, and identified a subset of 63 clones representative of this chromosome. The majority of clones were assigned to chromosome 21 by the presence of specific STSs and in situ hybridization. Twenty-nine of 36 STSs that we tested were detected in the subset, and a contig spanning 20 centimorgans in the genetic map and containing 8 STSs in 4 YACs was identified. The proposed approach can greatly speed efforts to construct physical maps of the human genome.     

Polymorphism for a 1.6-Mb deletion of the human Y chromosome persists through balance between recurrent mutation and haploid selection

Many human Y-chromosomal deletions are thought to severely impair reproductive fitness, which precludes their transmission to the next generation and thus ensures their rarity in the population. Here we report a 1.6-Mb deletion that persists over generations and is sufficiently common to be considered a polymorphism. We hypothesized that this deletion might affect spermatogenesis because it removes almost half of the Y chromosome's AZFc region, a gene-rich segment that is critical for sperm production. An association study established that this deletion, called gr/gr, is a significant risk factor for spermatogenic failure. The gr/gr deletion has far lower penetrance with respect to spermatogenic failure than previously characterized Y-chromosomal deletions; it is often transmitted from father to son. By studying the distribution of gr/gr-deleted chromosomes across the branches of the Y chromosome's genealogical tree, we determined that this deletion arose independently at least 14 times in human history. We suggest that the existence of this deletion as a polymorphism reflects a balance between haploid selection, which culls gr/gr-deleted Y chromosomes from the population, and homologous recombination, which continues to generate new gr/gr deletions.     

Mutational and functional analyses reveal that ST7 is a highly conserved tumor-suppressor gene on human chromosome 7q31

Loss of heterozygosity (LOH) of markers on human chromosome 7q31 is frequently encountered in a variety of human neoplasias, indicating the presence of a tumor-suppressor gene (TSG). By a combination of microcell-fusion and deletion-mapping studies, we previously established that this TSG resides within a critical region flanked by the genetic markers D7S522 and D7S677. Using a positional cloning strategy and aided by the availability of near-complete sequence of this genomic interval, we have identified a TSG within 7q31, named ST7 (for suppression of tumorigenicity 7; this same gene was recently reported in another context and called RAY1). ST7 is ubiquitously expressed in human tissues. Analysis of a series of cell lines derived from breast tumors and primary colon carcinomas revealed the presence of mutations in ST7. Introduction of the ST7 cDNA into the prostate-cancer-derived cell line PC3 had no effect on the in vitro proliferation of the cells, but abrogated their in vivo tumorigenicity. Our data indicate that ST7 is a TSG within chromosome 7q31 and may have an important role in the development of some types of human cancer.     

Dosage compensation of the active X chromosome in mammals

Monosomy of the X chromosome owing to divergence between the sex chromosomes leads to dosage compensation mechanisms to restore balanced expression between the X and the autosomes. In Drosophila melanogaster, upregulation of the male X leads to dosage compensation. It has been hypothesized that mammals likewise upregulate their active X chromosome. Together with X inactivation, this mechanism would maintain balanced expression between the X chromosome and autosomes and between the sexes. Here, we show that doubling of the global expression level of the X chromosome leads to dosage compensation in somatic tissues from several mammalian species. X-linked genes are highly expressed in brain tissues, consistent with a role in cognitive functions. Furthermore, the X chromosome is expressed but not upregulated in spermatids and secondary oocytes, preserving balanced expression of the genome in these haploid cells. Upon fertilization, upregulation of the active X must occur to achieve the observed dosage compensation in early embryos.     

Dynamics of the p53-Mdm2 feedback loop in individual cells

The tumor suppressor p53, one of the most intensely investigated proteins, is usually studied by experiments that are averaged over cell populations, potentially masking the dynamic behavior in individual cells. We present a system for following, in individual living cells, the dynamics of p53 and its negative regulator Mdm2 (refs. 1,4-7): this system uses functional p53-CFP and Mdm2-YFP fusion proteins and time-lapse fluorescence microscopy. We found that p53 was expressed in a series of discrete pulses after DNA damage. Genetically identical cells had different numbers of pulses: zero, one, two or more. The mean height and duration of each pulse were fixed and did not depend on the amount of DNA damage. The mean number of pulses, however, increased with DNA damage. This approach can be used to study other signaling systems and suggests that the p53-Mdm2 feedback loop generates a 'digital' clock that releases well-timed quanta of p53 until damage is repaired or the cell dies.     

Stochastic yet biased expression of multiple Dscam splice variants by individual cells

The Drosophila melanogaster gene Dscam is essential for axon guidance and has 38,016 possible alternative splice forms. This diversity can potentially be used to distinguish cells. We analyzed the Dscam mRNA isoforms expressed by different cell types and individual cells. The choice of splice variants expressed is regulated both spatially and temporally. Different subtypes of photoreceptors express broad yet distinctive spectra of Dscam isoforms. Single-cell RT-PCR documented that individual cells express several different Dscam isoforms and allowed an estimation of the diversity that is present. For example, we estimate that each R3/R4 photoreceptor cell expresses 14-50 distinct mRNAs chosen from the spectrum of thousands of splice variants distinctive of its cell type. Thus, the Dscam repertoire of each cell is different from those of its neighbors, providing a potential mechanism for generating unique cell identity in the nervous system and elsewhere.     

Genetic basis of individual differences in the response to small-molecule drugs in yeast

Individual response to small-molecule drugs is variable; a drug that provides a cure for some may confer no therapeutic benefit or trigger an adverse reaction in others. To begin to understand such differences systematically, we treated 104 genotyped segregants from a cross between two yeast strains with a collection of 100 diverse small molecules. We used linkage analysis to identify 124 distinct linkages between genetic markers and response to 83 compounds. The linked markers clustered at eight genomic locations, or quantitative-trait locus 'hotspots', that contain one or more polymorphisms that affect response to multiple small molecules. We also experimentally verified that a deficiency in leucine biosynthesis caused by a deletion of LEU2 underlies sensitivity to niguldipine, which is structurally related to therapeutic calcium channel blockers, and that a natural coding-region polymorphism in the inorganic phosphate transporter PHO84 underlies sensitivity to two polychlorinated phenols that uncouple oxidative phosphorylation. Our results provide a step toward a systematic understanding of small-molecule drug action in genetically distinct individuals.     

Whole-genome sequencing and comprehensive variant analysis of a Japanese individual using massively parallel sequencing

We report the analysis of a Japanese male using high-throughput sequencing to × 40 coverage. More than 99% of the sequence reads were mapped to the reference human genome. Using a Bayesian decision method, we identified 3,132,608 single nucleotide variations (SNVs). Comparison with six previously reported genomes revealed an excess of singleton nonsense and nonsynonymous SNVs, as well as singleton SNVs in conserved non-coding regions. We also identified 5,319 deletions smaller than 10 kb with high accuracy, in addition to copy number variations and rearrangements. De novo assembly of the unmapped sequence reads generated around 3 Mb of novel sequence, which showed high similarity to non-reference human genomes and the human herpesvirus 4 genome. Our analysis suggests that considerable variation remains undiscovered in the human genome and that whole-genome sequencing is an invaluable tool for obtaining a complete understanding of human genetic variation.