Populationsgenetik der ersten Bauern Mitteleuropas

Die vorliegende Dissertation ist eine molekulargenetische Studie an humanem neolithischem Skelettmaterial. Im zentralen Blickpunkt stand die Bestimmung der Variabilität der mitochondrialen Haplogruppen einer frühneolithischen Stichprobe aus drei unterschiedlichen Kulturkreisen, welche die Linearbandkeramik (LBK und AVK), die Körös-Kultur und eine Sammelkategorie osteuropäischer spätmeso- und frühneolithischer Kulturen umfasste. Im Vergleich dieser Gruppen untereinander sowie mit Rezentdaten moderner Populationen aus vergleichbaren Gebieten Mittel- und Osteuropas sowie dem Nahen Osten sollten bestehende Modelle und Hypothesen zur Neolithisierung Mitteleuropas geprüft werden. Insgesamt konnte für 43 neolithische Individuen aus 16 Fundorten der reproduzierbare Nachweis endogener DNA erbracht werden. Eine eindeutige Haplogruppenbestimmung konnte durch die Sequenzierung vier überlappender Fragmente der mitochondrialen Hypervariablen Region I sowie durch RFLP-Analyse zusätzlicher charakteristischer Nukleotidpositionen für alle 43 Individuen durchgeführt werden. Die neolithischen Individuen der Linearbandkeramik sowie der Körös-Kultur zeigten eine hohe Diversität an bekannten europäischen Haplogruppen, wohingegen die kleinere Stichprobe aus dem Gebiet Osteuropas eine auffällige Homogenität aufwies. Neben Frequenzunterschieden zur modernen mitteleuropäischen Bevölkerung war innerhalb der LBK/AVK-Stichprobe eine hohe Frequenz der Haplogruppe N1a festzustellen, welche nicht in den beiden anderen neolithischen Stichproben zu finden war und auch in der heutigen Rezentbevölkerung Eurasiens und Nordafrikas nur mit einer durchschnittlichen Frequenz von 0,2% vertreten ist. Innerhalb der Individuen der Körös-Kultur fanden sich zwei Haplotypen, welche heute nicht in Europa bekannt sind, dagegen jedoch in Süd- bzw. Nordostasien gehäuft vorkommen. Die Ergebnisse der aDNA-Analysen bestätigten im Wesentlichen das komplexe Bild der Neolithischen Transition in Mitteleuropa und konnten die, für diesen Raum postulierte, Hypothese der leap frog colonization weitestgehend unterstützen. Auch für den geographischen Vergleichsraum des nördlichen Osteuropa konnten Parallelen zur etablierten Sichtweise der archäologischen Forschung zu diesem Gebiet und den vorliegenden Ergebnissen der aDNA-Analysen aufgezeigt werden. Die zeitlich verzögerte Annahme der neolithischen Lebensweise im waldreichen nördlichen Osteuropa spiegelt sich in der reduzierten Diversität an mtDNA-Haplogruppen wider. Die vorliegende Dissertation konnte nicht nur durch die Ergebnisse der Haplogruppen-Bestimmung, sondern vor allem durch die umfangreichen und elaborierten Reproduktions- und Authentifizierungprotokolle deutlich machen, dass der Nachweis von humaner alter DNA trotz der allgegenwärtigen, methodenimmanenten Kontaminationsgefahr unter streng kontrollierten Bedingungen möglich ist. Gleichermaßen konnte veranschaulicht werden, dass die aDNA-Analyse wertvolle Hinweise auf das genetische status quo einer Population liefern kann, welche nicht bzw. nur in sehr eingeschränkten Maße von rezenten DNA-Daten abgeleitet werden können. Als sekundäres Ergebnis erlaubte der bislang größte vorliegende Datensatz von ~2500 Klonsequenzen zudem einen detaillierten Einblick in Häufigkeiten und Verteilungsmuster von post mortem Sequenzveränderungen. Es konnten für den mitochondrialen Bereich der Nukleotidpositionen 15997-16409 so genannte hot bzw. cold spots definiert werden, welche für die Auswertung und Interpretation von zukünftigen Sequenzierungen der Hypervariablen Region I des mt-Genoms von entscheidender Bedeutung sein werden.

Evidence for selection on coloration in a Panamanian poison frog: a coalescent-based approach

Aim The strawberry poison frog, Oophaga pumilio, has undergone a remarkable radiation of colour morphs in the Bocas del Toro archipelago in Panama. This species shows extreme variation in colour and pattern between populations that have been geographically isolated for < 10,000 years. While previous research has suggested the involvement of divergent selection, to date no quantitative test has examined this hypothesis. Location Bocas del Toro archipelago, Panama. Methods We use a combination of population genetics, phylogeography and phenotypic analyses to test for divergent selection in coloration in O. pumilio. Tissue samples of 88 individuals from 15 distinct populations were collected. Using these data, we developed a gene tree using the mitochondrial DNA (mtDNA) d-loop region. Using parameters derived from our mtDNA phylogeny, we predicted the coalescence of a hypothetical nuclear gene underlying coloration. We collected spectral reflectance and body size measurements on 94 individuals from four of the populations and performed a quantitative analysis of phenotypic divergence. Results The mtDNA d-loop tree revealed considerable polyphyly across populations. Coalescent reconstructions of gene trees within population trees revealed incomplete genotypic sorting among populations. The quantitative analysis of phenotypic divergence revealed complete lineage sorting by colour, but not by body size: populations showed non-overlapping variation in spectral reflectance measures of body coloration, while variation in body size did not separate populations. Simulations of the coalescent using parameter values derived from our empirical analyses demonstrated that the level of sorting among populations seen in colour cannot reasonably be attributed to drift. Main conclusions These results imply that divergence in colour, but not body size, is occurring at a faster rate than expected under neutral processes. Our study provides the first quantitative support for the claim that strong diversifying selection underlies colour variation in the strawberry poison frog.

ABCtoolbox: a versatile toolkit for approximate Bayesian computations

Background The estimation of demographic parameters from genetic data often requires the computation of likelihoods. However, the likelihood function is computationally intractable for many realistic evolutionary models, and the use of Bayesian inference has therefore been limited to very simple models. The situation changed recently with the advent of Approximate Bayesian Computation (ABC) algorithms allowing one to obtain parameter posterior distributions based on simulations not requiring likelihood computations. Results Here we present ABCtoolbox, a series of open source programs to perform Approximate Bayesian Computations (ABC). It implements various ABC algorithms including rejection sampling, MCMC without likelihood, a Particle-based sampler and ABC-GLM. ABCtoolbox is bundled with, but not limited to, a program that allows parameter inference in a population genetics context and the simultaneous use of different types of markers with different ploidy levels. In addition, ABCtoolbox can also interact with most simulation and summary statistics computation programs. The usability of the ABCtoolbox is demonstrated by inferring the evolutionary history of two evolutionary lineages of Microtus arvalis. Using nuclear microsatellites and mitochondrial sequence data in the same estimation procedure enabled us to infer sex-specific population sizes and migration rates and to find that males show smaller population sizes but much higher levels of migration than females. Conclusion ABCtoolbox allows a user to perform all the necessary steps of a full ABC analysis, from parameter sampling from prior distributions, data simulations, computation of summary statistics, estimation of posterior distributions, model choice, validation of the estimation procedure, and visualization of the results.

East Asian ethnolinguistic phylogeography

A polyphyletic understanding of Asian linguistic diversity was first propagated in 1823. Since 1901, various scholars have proposed larger linguistic phyla uniting two or more recognised Asian language families. The most recent proposal in this tradition, Starosta’s 2001 East Asian phylum, comprising the Trans-Himalayan, Hmong-Mien, Austroasiatic, Austronesian and Kradai language families, is reassessed in light of linguistic and non-linguistic evidence. Ethnolinguistically informed inferences based on Asian Y chromosomal phylogeography lead to a reconstruction of various episodes of ethnolinguistic prehistory which lie beyond the linguistic event horizon, i.e. at a time depth empirically inaccessible to historical linguistics. The Father Tongue correlation in population genetics, the evidence for refugia during the Last Glacial Maximum and the hypothesis of language families having arisen as the result of demographic bottlenecks in prehistory are shown to be crucial to an understanding of the ethnogenesis of East Asian linguistic phyla. The prehistory of several neighbouring Asian language families is discussed, and the Centripetal Migration model is opposed to the Farming Language Dispersal theory.

Molecular epidemiology of vancomycin-resistant Enterococcus faecium: a prospective, multicenter study in South American hospitals.

Enterococcus faecium has emerged as an important nosocomial pathogen worldwide, and this trend has been associated with the dissemination of a genetic lineage designated clonal cluster 17 (CC17). Enterococcal isolates were collected prospectively (2006 to 2008) from 32 hospitals in Colombia, Ecuador, Perú, and Venezuela and subjected to antimicrobial susceptibility testing. Genotyping was performed with all vancomycin-resistant E. faecium (VREfm) isolates by pulsed-field gel electrophoresis (PFGE) and multilocus sequence typing. All VREfm isolates were evaluated for the presence of 16 putative virulence genes (14 fms genes, the esp gene of E. faecium [espEfm], and the hyl gene of E. faecium [hylEfm]) and plasmids carrying the fms20-fms21 (pilA), hylEfm, and vanA genes. Of 723 enterococcal isolates recovered, E. faecalis was the most common (78%). Vancomycin resistance was detected in 6% of the isolates (74% of which were E. faecium). Eleven distinct PFGE types were found among the VREfm isolates, with most belonging to sequence types 412 and 18. The ebpAEfm-ebpBEfm-ebpCEfm (pilB) and fms11-fms19-fms16 clusters were detected in all VREfm isolates from the region, whereas espEfm and hylEfm were detected in 69% and 23% of the isolates, respectively. The fms20-fms21 (pilA) cluster, which encodes a putative pilus-like protein, was found on plasmids from almost all VREfm isolates and was sometimes found to coexist with hylEfm and the vanA gene cluster. The population genetics of VREfm in South America appear to resemble those of such strains in the United States in the early years of the CC17 epidemic. The overwhelming presence of plasmids encoding putative virulence factors and vanA genes suggests that E. faecium from the CC17 genogroup may disseminate in the region in the coming years.

Dominance and the maintenance of polymorphism in multiallelic migration-selection models with two demes

The maintenance of genetic variation in a spatially heterogeneous environment has been one of the main research themes in theoretical population genetics. Despite considerable progress in understanding the consequences of spatially structured environments on genetic variation, many problems remain unsolved. One of them concerns the relationship between the number of demes, the degree of dominance, and the maximum number of alleles that can be maintained by selection in a subdivided population. In this work, we study the potential of maintaining genetic variation in a two-deme model with deme-independent degree of intermediate dominance, which includes absence of G x E interaction as a special case. We present a thorough numerical analysis of a two-deme three-allele model, which allows us to identify dominance and selection patterns that harbor the potential for stable triallelic equilibria. The information gained by this approach is then used to construct an example in which existence and asymptotic stability of a fully polymorphic equilibrium can be proved analytically. Noteworthy, in this example the parameter range in which three alleles can coexist is maximized for intermediate migration rates. Our results can be interpreted in a specialist-generalist context and (among others) show when two specialists can coexist with a generalist in two demes if the degree of dominance is deme independent and intermediate. The dominance relation between the generalist allele and the specialist alleles play a decisive role. We also discuss linear selection on a quantitative trait and show that G x E interaction is not necessary for the maintenance of more than two alleles in two demes.

Distance from Sub-Saharan Africa Predicts Mutational Load in Diverse Human Genomes

The Out-of-Africa (OOA) dispersal ∼50,000 y ago is characterized by a series of founder events as modern humans expanded into multiple continents. Population genetics theory predicts an increase of mutational load in populations undergoing serial founder effects during range expansions. To test this hypothesis, we have sequenced full genomes and high-coverage exomes from seven geographically divergent human populations from Namibia, Congo, Algeria, Pakistan, Cambodia, Siberia, and Mexico. We find that individual genomes vary modestly in the overall number of predicted deleterious alleles. We show via spatially explicit simulations that the observed distribution of deleterious allele frequencies is consistent with the OOA dispersal, particularly under a model where deleterious mutations are recessive. We conclude that there is a strong signal of purifying selection at conserved genomic positions within Africa, but that many predicted deleterious mutations have evolved as if they were neutral during the expansion out of Africa. Under a model where selection is inversely related to dominance, we show that OOA populations are likely to have a higher mutation load due to increased allele frequencies of nearly neutral variants that are recessive or partially recessive.

Development and characterization of novel microsatellite markers for Arion slug species

Seventeen polymorphic microsatellite markers were isolated and characterized in Arion vulgaris/lusitanicus, which belongs to the worst European slug pests with serious economic and ecological impact. These markers were tested on 23 individuals collected in a population in Switzerland. Numbers of alleles ranged from 2 to 14 per locus, observed and expected heterozygosities ranged from 0.174 to 0.87, and from 0.162 to 0.903, respectively. These loci were also successfully amplified and were polymorphic in the closely related species A. rufus and A. ater. These loci represent the first highly polymorphic nuclear markers described for A. vulgaris and pave the way for population genetics and molecular ecology research of the important Arion pest slugs.

A coalescent analysis for modeling the mutation process in colorectal cancer

Colorectal cancer is the forth most common diagnosed cancer in the United States. Every year about a hundred forty-seven thousand people will be diagnosed with colorectal cancer and fifty-six thousand people lose their lives due to this disease. Most of the hereditary nonpolyposis colorectal cancer (HNPCC) and 12% of the sporadic colorectal cancer show microsatellite instability. Colorectal cancer is a multistep progressive disease. It starts from a mutation in a normal colorectal cell and grows into a clone of cells that further accumulates mutations and finally develops into a malignant tumor. In terms of molecular evolution, the process of colorectal tumor progression represents the acquisition of sequential mutations. ^ Clinical studies use biomarkers such as microsatellite or single nucleotide polymorphisms (SNPs) to study mutation frequencies in colorectal cancer. Microsatellite data obtained from single genome equivalent PCR or small pool PCR can be used to infer tumor progression. Since tumor progression is similar to population evolution, we used an approach known as coalescent, which is well established in population genetics, to analyze this type of data. Coalescent theory has been known to infer the sample's evolutionary path through the analysis of microsatellite data. ^ The simulation results indicate that the constant population size pattern and the rapid tumor growth pattern have different genetic polymorphic patterns. The simulation results were compared with experimental data collected from HNPCC patients. The preliminary result shows the mutation rate in 6 HNPCC patients range from 0.001 to 0.01. The patients' polymorphic patterns are similar to the constant population size pattern which implies the tumor progression is through multilineage persistence instead of clonal sequential evolution. The results should be further verified using a larger dataset. ^

A computerized statistical framework for coalescent analysis

Coalescent theory represents the most significant progress in theoretical population genetics in the past three decades. The coalescent theory states that all genes or alleles in a given population are ultimately inherited from a single ancestor shared by all members of the population, known as the most recent common ancestor. It is now widely recognized as a cornerstone for rigorous statistical analyses of molecular data from population [1]. The scientists have developed a large number of coalescent models and methods[2,3,4,5,6], which are not only applied in coalescent analysis and process, but also in today’s population genetics and genome studies, even public health. The thesis aims at completing a statistical framework based on computers for coalescent analysis. This framework provides a large number of coalescent models and statistic methods to assist students and researchers in coalescent analysis, whose results are presented in various formats as texts, graphics and printed pages. In particular, it also supports to create new coalescent models and statistical methods. ^

Description of two mass species of Pacific copepods in terms of two polymorphic loci

Population genetics of two species of mass copepods Undinula darwini and Calanus australis, with different range types, is investigated. Both species exhibit considerable genetic diversity, especially C. australis (observed heterozygoticity = 0.36), which inhabits a variable biotope in the zone of the Peru current. Samples of both species exhibited highly significant genetic heterogeneity as well as heterozygote deficiency compared with the situation expected from the Hardy-Weinberg law. Contribution of distance isolation to genetic differentiation of populations is estimated. Gene drift is discussed as a source of heterogeneity in populations of planktic copepods. Possible aspects of population genetic research on marine plank-tic crustaceans are discussed.

Cloning of inversion breakpoints in the Anopheles gambiae complex traces a transposable element at the inversion junction

Anopheles arabiensis, one of the two most potent malaria vectors of the gambiae complex, is characterized by the presence of chromosomal paracentric inversions. Elucidation of the nature and the dynamics of these inversions is of paramount importance for the understanding of the population genetics and evolutionary biology of this mosquito and of the impact on malaria epidemiology. We report here the cloning of the breakpoints of the naturally occurring polymorphic inversion 2Rd′ of A. arabiensis. A cDNA clone that cytologically mapped on the proximal breakpoint was the starting material for the isolation of a cosmid clone that spanned the breakpoint. Analysis of the surrounding sequences demonstrated that adjacent to the distal breakpoint lies a repetitive element that exhibits distinct distribution in different A. arabiensis strains. Sequencing analysis of that area revealed elements characteristic of transposable element terminal repeats. We called this presumed transposable element Odysseus. The presence of Odysseus at the junction of the naturally occuring inversion 2Rd′ suggests that the inversion may be the result of the transposable element’s activity. Characteristics of Odysseus’ terminal region as well as its cytological distribution in different strains may indicate a relatively recent activity of Odysseus.

Resolving taxonomic uncertainty in vulnerable elasmobranchs : are the Madeira skate (Raja maderensis) and the thornback ray (Raja clavata) distinct species?

We are grateful to all those who helped with sample collection. This includes the MEDITS survey programme (IEO Mallorca) for Mediterranean samples. Portugal mainland samples were collected under the EU Data Collection Framework (DCF, PNAB). Azores specimens from the Department of Oceanography and Fisheries (DOP) of the University of the Azores (UAc) were collected under the project DEMERSAIS “Monitorização das espécies demersais dos Açores” financed by the Azorean government, and the project DEECON “Unravelling population connectivity for sustainable fisheries in the Deep Sea” project approved by the European Science Foundation (ESF) under the EUROCORES programme (proposal No 06-EuroDEEP-FP-008 & SFRH-EuroDEEP/0002/2007). This study was funded by the UK Natural Environment Research Council (NERC) Oceans 2025 Strategic Research Programme Theme 6 (Science for Sustainable Marine Resources). REB was supported by the Fisheries Society of the British Isles (FSBI), BSP was funded by the Fundação para a Ciência e a Tecnologia, SFRH/BPD/72351/2010 and JL was supported by The Alasdair Downes Marine Conservation Fund.