Gunter Weiss

A Neutral Model of Transcriptome Evolution

Microarray technologies allow the identification of large numbers of expression differences within and between species. Although environmental and physiological stimuli are clearly responsible for changes in the expression levels of many genes, it is not known whether the majority of changes of gene expression fixed during evolution between species and between various tissues within a species are caused by Darwinian selection or by stochastic processes. We find the following: (1) expression differences between species accumulate approximately linearly with time; (2) gene expression variation among individuals within a species correlates positively with expression divergence between species; (3) rates of expression divergence between species do not differ significantly between intact genes and expressed pseudogenes; (4) expression differences between brain regions within a species have accumulated approximately linearly with time since these regions emerged during evolution. These results suggest that the majority of expression differences observed between species are selectively neutral or nearly neutral and likely to be of little or no functional significance. Therefore, the identification of gene expression differences between species fixed by selection should be based on null hypotheses assuming functional neutrality. Furthermore, it may be possible to apply a molecular clock based on expression differences to infer the evolutionary history of tissues.

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.

Reduced Y-chromosome, but not mitochondrial DNA, diversity in human populations from West New Guinea.

To investigate the paternal population history of New Guinea, 183 individuals from 11 regional populations of West New Guinea (WNG) and 131 individuals from Papua New Guinea (PNG) were analyzed at 26 binary markers and seven short-tandem-repeat loci from the nonrecombining part of the human Y chromosome and were compared with 14 populations of eastern and southeastern Asia, Polynesia, and Australia. Y-chromosomal diversity was low in WNG compared with PNG and with most other populations from Asia/Oceania; a single haplogroup (M-M4) accounts for 75% of WNG Y chromosomes, and many WNG populations have just one Y haplogroup. Four Y-chromosomal lineages (haplogroups M-M4, C-M208, C-M38, and K-M230) account for 94% of WNG Y chromosomes and 78% of all Melanesian Y chromosomes and were identified to have most likely arisen in Melanesia. Haplogroup C-M208, which in WNG is restricted to the Dani and Lani, two linguistically closely related populations from the central and western highlands of WNG, was identified as the major Polynesian Y-chromosome lineage. A network analysis of associated Y-chromosomal short-tandem-repeat haplotypes suggests two distinct population expansions involving C-M208--one in New Guinea and one in Polynesia. The observed low levels of Y-chromosome diversity in WNG contrast with high levels of mtDNA diversity reported for the same populations. This most likely reflects extreme patrilocality and/or biased male reproductive success (polygyny). Our data further provide evidence for primarily female-mediated gene flow within the highlands of New Guinea but primarily male-mediated gene flow between highland and lowland/coastal regions.

Independent histories of human y chromosomes from melanesia and Australia

To investigate the origins and relationships of Australian and Melanesian populations, 611 males from 18 populations from Australia, Melanesia, and eastern/southeastern Asia were typed for eight single-nucleotide polymorphism (SNP) loci and seven short tandem-repeat loci on the Y chromosome. A unique haplotype, DYS390.1del/RPS4Y711T, was found at a frequency of 53%-69% in Australian populations, whereas the major haplotypes found in Melanesian populations (M4G/M5T/M9G and DYS390.3del/RPS4Y711T) are absent from the Australian populations. The Y-chromosome data thus indicate independent histories for Australians and Melanesians, a finding that is in agreement with evidence from mtDNA but that contradicts some analyses of autosomal loci, which show a close relationship between Australian and Melanesian (specifically, highland Papua New Guinean) populations. Since the Australian and New Guinean landmasses were connected when first colonized by humans > or =50,000 years ago but separated some 8,000 years ago, a possible way to reconcile all the genetic data is to infer that the Y-chromosome and mtDNA results reflect the past 8,000 years of independent history for Australia and New Guinea, whereas the autosomal loci reflect the long preceding period of common origin and shared history. Two Y-chromosome haplotypes (M119C/M9G and M122C/M9G) that originated in eastern/southeastern Asia are present in coastal and island Melanesia but are rare or absent in both Australia and highland Papua New Guinea. This distribution, along with demographic analyses indicating that population expansions for both haplotypes began approximately 4,000-6,000 years ago, suggests that these haplotypes were brought to Melanesia by the Austronesian expansion. Most of the populations in this study were previously typed for mtDNA SNPs; population differentiation is greater for the Y chromosome than for mtDNA and is significantly correlated with geographic distance, a finding in agreement with results of similar analyses of European populations.

FUNC: A package for detecting significant associations between gene sets and ontological annotations

BackgroundGenome-wide expression, sequence and association studies typically yield large sets of gene candidates, which must then be further analysed and interpreted. Information about these genes is increasingly being captured and organized in ontologies, such as the Gene Ontology. Relationships between the gene sets identified by experimental methods and biological knowledge can be made explicit and used in the interpretation of results. However, it is often difficult to assess the statistical significance of such analyses since many inter-dependent categories are tested simultaneously.ResultsWe developed the program package FUNC that includes and expands on currently available methods to identify significant associations between gene sets and ontological annotations. Implemented are several tests in particular well suited for genome wide sequence comparisons, estimates of the family-wise error rate, the false discovery rate, a sensitive estimator of the global significance of the results and an algorithm to reduce the complexity of the results.ConclusionFUNC is a versatile and useful tool for the analysis of genome-wide data. It is freely available under the GPL license and also accessible via a web service.

Sudden replacement of cave bear mitochondrial DNA in the late Pleistocene.

In the absence of interaction with genetically distinct populations, changes in frequencies of distinct mitochondrial DNA (mtDNA) sequences (haplotypes) within a population are caused by mutations and genetic drift, both of which are comparatively slow processes. By contrast, interactions between populations may lead to rapid changes in haplotype frequencies [1–3], for instance through extinction and recolonialisation. So far, the only case of a direct replacement of mtDNA sequences within a continuous population has been described in historical mouse populations [2].

Recent Origin and Cultural Reversion of a Hunter–Gatherer Group

Contemporary hunter–gatherer groups are often thought to serve as models of an ancient lifestyle that was typical of human populations prior to the development of agriculture. Patterns of genetic variation in hunter–gatherer groups such as the !Kung and African Pygmies are consistent with this view, as they exhibit low genetic diversity coupled with high frequencies of divergent mtDNA types not found in surrounding agricultural groups, suggesting long-term isolation and small population sizes. We report here genetic evidence concerning the origins of the Mlabri, an enigmatic hunter–gatherer group from northern Thailand. The Mlabri have no mtDNA diversity, and the genetic diversity at Y-chromosome and autosomal loci are also extraordinarily reduced in the Mlabri. Genetic, linguistic, and cultural data all suggest that the Mlabri were recently founded, 500–800 y ago, from a very small number of individuals. Moreover, the Mlabri appear to have originated from an agricultural group and then adopted a hunting–gathering subsistence mode. This example of cultural reversion from agriculture to a hunting–gathering lifestyle indicates that contemporary hunter–gatherer groups do not necessarily reflect a pre-agricultural lifestyle.

Modeling the polymerase chain reaction.

We introduce a mathematical model to treat the polymerase chain reaction (PCR), where we regard the accumulation of new molecules during a PCR cycle as a randomly bifurcating tree. This model enables us to compute an approximate formula for the distribution of the number of replications that have occurred between a pair of molecules, which depends on the efficiency lambda of the reaction, the number N0 of template molecules at the beginning of the PCR and the number c of PCR cycles. The reliability of the approximation is tested by computer simulations. Finally, to model the effect of the intrinsic error rate of the polymerase, we superimpose a substitution process on the tree. The resulting closed formula for the distribution of pairwise differences of sequences as a function of error rate mu and efficiency lambda can be used to estimate the error rate, if lambda is known.