Phylogeographic structure of the fire ant Solenopsis invicta in its native South American range: roles of natural barriers and habitat connectivity

We generated mitochondrial DNA (mtDNA) sequence data from 402 individuals of the fire ant Solenopsis invicta collected from 11 native populations and analyzed these data using a combination of demographic, phylogenetic, and phylogeographic methods to infer features of the evolutionary history of this species. Prior expectations regarding high levels of genetic structure and isolation by distance among populations were supported by the data, but we also discovered several unanticipated patterns. Our analyses revealed a major genetic break between S. invicta mtDNA haplotypes that coincides with the Mesopotamia wetlands region of South America, resulting in two higher level nested clade groupings. In addition, we identified contrasting patterns of genetic differentiation within these two major groups, which may reflect differences in connectivity of suitable habitat in different parts of the native range of S. invicta. Our study represents the first attempt to understand the phylogeographic history of S. invicta across its native range.     

Dissecting the genetic structure of Korean population using genome-wide SNP arrays

Genome-wide SNP arrays have generated unprecedented quantities of data allow the detection of human evolutionary history and dense genome-wide data also enable the identification of distance ancestry among individuals or ethnic groups. To explain wider aspects of the genetic structure of Koreans and the East Asian population, we analyzed 79 individuals from the Korean HapMap project at 555,352 common single-nucleotide polymorphism loci, and compared this data with the worldwide population groups with the 53 ethnic groups from Human Genome Diversity Panel (HGDP-CEPH). Population differentiation (F_ST), Principal Component Analyses, STRUCTURE and ADMIXTURE are examined. In general, all the individual samples studies here were classified into subset of ethnic groups according to their geographical origins. Korean HapMap individuals were grouped together with East Asian populations from HGDP panel. Recently, a sub-population structure within Korean population has been reported. Our result, however, revealed the genetic homogeneity of Korean population. The ADMIXTURE analysis showed that, overall the Korean populations derive 79 % of their genomic ancestry from southern Asia and have relatively little northern Asian ancestry (21 %). The present work, therefore, provide the evidence that the male-biased southern-to-northern migration influenced not only for the genetic make up of the Y chromosome in the Korean population but also, its autosomal composition.     

The impact of founder effects, gene flow, and European admixture on native American genetic diversity

Recent studies have concluded that the global pattern of neutral genetic diversity in humans reflects a series of founder effects and population movements associated with our recent expansion out of Africa. In contrast, regional studies tend to emphasize the significance of more complex patterns of colonization, gene flow, and secondary population movements in shaping patterns of diversity. Our objective in this study is to examine how founder effects, gene flow, and European admixture have molded patterns of neutral genetic diversity in the Americas. Our strategy is to test the fit of a serial founder effects process to the pattern of neutral autosomal genetic variation and to examine the contribution of gene flow and European admixture to departures from fit. The genetic data consist of 678 autosomal microsatellite loci assayed by Wang and colleagues in 530 individuals in 29 widely distributed Native American populations. We find that previous evidence for serial founder effects in the Americas may be driven in part by high levels of European admixture in northern North America, intermediate levels in Central America, and low levels in eastern South America. Geographically patterned admixture may also account for previously reported genetic differences between Andean and Amazonian groups. Though admixture has obscured the precise details of precontact evolutionary processes, we find that genetic diversity is still largely hierarchically structured and that gene flow between neighboring groups has had surprisingly little impact on macrogeographic patterns of genetic diversity in the Americas.     

Predicting demographic group structures based on DNA sequence data

The ability to infer relationships between groups of sequences, either by searching for their evolutionary history or by comparing their sequence similarity, can be a crucial step in hypothesis testing. Interpreting relationships of human immunodeficiency virus type 1 (HIV-1) sequences can be challenging because of their rapidly evolving genomes, but it may also lead to a better understanding of the underlying biology. Several studies have focused on the evolution of HIV-1, but there is little information to link sequence similarities and evolutionary histories of HIV-1 to the epidemiological information of the infected individual. Our goal was to correlate patterns of HIV-1 genetic diversity with epidemiological information, including risk and demographic factors. These correlations were then used to predict epidemiological information through analyzing short stretches of HIV-1 sequence. Using standard phylogenetic and phenetic techniques on 100 HIV-1 subtype B sequences, we were able to show some correlation between the viral sequences and the geographic area of infection and the risk of men who engage in sex with men. To help identify more subtle relationships between the viral sequences, the method of multidimensional scaling (MDS) was performed. That method identified statistically significant correlations between the viral sequences and the risk factors of men who engage in sex with men and individuals who engage in sex with injection drug users or use injection drugs themselves. Using tree construction, MDS, and newly developed likelihood assignment methods on the original 100 samples we sequenced, and also on a set of blinded samples, we were able to predict demographic/risk group membership at a rate statistically better than by chance alone. Such methods may make it possible to identify viral variants belonging to specific demographic groups by examining only a small portion of the HIV-1 genome. Such predictions of demographic epidemiology based on sequence information may become valuable in assigning different treatment regimens to infected individuals.     

Is Cupressus sempervirens native in Italy? An answer from genetic and palaeobotanical data

This study represents the first large-scale analysis using nuclear molecular markers to assess genetic diversity and structure of Cupressus sempervirens L.. Genetic and fossil data were combined to infer the possible role of human activity and evolutionary history in shaping the diversity of cypress populations. We analysed 30 populations with six polymorphic nuclear microsatellite markers. Dramatic reductions in heterozygosity and allelic richness were observed from east to west across the species range. Structure analysis assigned individuals to two main groups separating central Mediterranean and eastern populations. The two main groups could be further divided into five subgroups which showed the following geographical distributions: Turkey with the Greek islands Rhodes and Samos, Greece (Crete), Southern Italy, Northern Italy, Tunisia with Central Italy. This pattern of genetic structure is also supported by samova and Barrier analyses. Palaeobotanical data indicated that Cupressus was present in Italy in the Pliocene, Pleistocene and Holocene. Furthermore, our molecular survey showed that Italian cypress populations experienced bottlenecks that resulted in reduced genetic diversity and allelic richness and greater genetic differentiation. Recent colonization or introduction may also have influenced levels of diversity detected in the Italian populations, as most individuals found in this range today have multilocus genotypes that are also present in the eastern range of the species. The data reveal a new interpretation of the history of cypress distribution characterized by ancient eastern populations (Turkey and Greek islands) and a mosaic of recently introduced trees and remnants of ancient, depauperate populations in the central Mediterranean range.     

Evolutionary relationships among parvoviruses: virus-host coevolution among autonomous primate parvoviruses and links between adeno-associated and avian parvoviruses

The current classification of parvoviruses is based on virus host range and helper virus dependence, while little data on evolutionary relationships among viruses are available. We identified and analyzed 472 sequences of parvoviruses, among which there were (virtually) full-length genomes of all 41 viruses currently recognized as individual species within the family Parvoviridae. Our phylogenetic analysis of full-length genomes as well as open reading frames distinguished three evolutionary groups of parvoviruses from vertebrates: (i) the human helper-dependent adeno-associated virus (AAV) serotypes 1 to 6 and the autonomous avian parvoviruses; (ii) the bovine, chipmunk, and autonomous primate parvoviruses, including human viruses B19 and V9; and (iii) the parvoviruses from rodents (except for chipmunks), carnivores, and pigs. Each of these three evolutionary groups could be further subdivided, reflecting both virus-host coevolution and multiple cross-species transmissions in the evolutionary history of parvoviruses. No parvoviruses from invertebrates clustered with vertebrate parvoviruses. Our analysis provided evidence for negative selection among parvoviruses, the independent evolution of their genes, and recombination among parvoviruses from rodents. The topology of the phylogenetic tree of autonomous human and simian parvoviruses matched exactly the topology of the primate family tree, as based on the analysis of primate mitochondrial DNA. Viruses belonging to the AAV group were not evolutionarily linked to other primate parvoviruses but were linked to the parvoviruses of birds. The two lineages of human parvoviruses may have resulted from independent ancient zoonotic infections. Our results provide an argument for reclassification of Parvovirinae based on evolutionary relationships among viruses.     

Phylogenetic relationship analysis of Iranians and other world populations using allele frequencies at 12 polymorphic markers

The estimation of genetic distance between populations could improve our viewpoint about human migration and its genetic origin. In this study, we used allele frequency data of 12 polymorphic markers on 250 individuals (500 alleles) from the Iranian population to estimate genetic distance between the Iranians and other world populations. The phylogenetic trees for three different sets of allele frequency data were constructed. Our results revealed the genetic similarity between the Iranians and European populations. The lowest genetic distance was observed between the Iranians and some populations reside in Russia. Furthermore, the high genetic distance was observed between the Iranians and East Asian populations. The data suggested that the Iranians might have relatively close evolutionary history with Europeans, but historically independent from East Asian populations. The evaluation of genetic distance between Indians populations and Iranians was also performed. The Indian groups showed low genetic distance with others, but high genetic distance with the Iranians. This study could provide a new insight into the evolutionary history of the Iranian population.     

Estimating divergence times from DNA sequences

The patterns of genetic variation within and among individuals and populations can be used to make inferences about the evolutionary forces that generated those patterns. Numerous population genetic approaches have been developed in order to infer evolutionary history. Here, we present the “Two-Two (TT)” and the “Two-Two-outgroup (TTo)” methods; two closely related approaches for estimating divergence time based in coalescent theory. They rely on sequence data from two haploid genomes (or a single diploid individual) from each of two populations. Under a simple population-divergence model, we derive the probabilities of the possible sample configurations. These probabilities form a set of equations that can be solved to obtain estimates of the model parameters, including population split times, directly from the sequence data. This transparent and computationally efficient approach to infer population divergence time makes it possible to estimate time scaled in generations (assuming a mutation rate), and not as a compound parameter of genetic drift. Using simulations under a range of demographic scenarios, we show that the method is relatively robust to migration and that the TTo method can alleviate biases that can appear from drastic ancestral population size changes. We illustrate the utility of the approaches with some examples, including estimating split times for pairs of human populations as well as providing further evidence for the complex relationship among Neandertals and Denisovans and their ancestors.     

Mitochondrial phylogeny of African wood mice, genus Hylomyscus (Rodentia, Muridae): implications for their taxonomy and biogeography

This paper investigates the usefulness of two mitochondrial genes (16S rRNA and cytochrome b) to solve taxonomical difficulties within the genus Hylomyscus and to infer its evolutionary history. Both genes proved to be suitable molecular markers for diagnosis of Hylomyscus species. Nevertheless the resolving powers of these two genes differ, and with both markers (either analyzed singly or in combination), some nodes remain unresolved. This is probably related to the fact that the species emerged during a rapid diversification event that occurred 2-6 Myr ago (4-5 Myr ago for most divergence events). Our molecular data support the recognition of an "aeta" group, while the "alleni" and "parvus" groups are not fully supported. Based on tree topology and genetic divergence, two taxa generally recognized as subspecies should be elevated at the species level (H. simus and H. cf kaimosae). H. stella populations exhibit ancient haplotype segregation that may represent currently unrecognized allopatric species. The existence of cryptic species within H. parvus is questioned. Finally, three potentially new species may occur in West Central Africa. The Congo and Oubangui Rivers, as well as the Volta and Niger Rivers and/or the Dahomey gap could have formed effective barriers to Hylomyscus species dispersal, favoring their speciation in allopatry. The pronounced shifts in African climate during the late Pliocene and Miocene, which resulted in major changes in the distribution and composition of the vegetation, could have promoted speciation within the genus (refuge theory). Future reports should focus on the geographic distribution of Hylomyscus species in order to get a better understanding of the evolutionary history of the genus.     

Genetic variation of recent Alu insertions in human populations

The Alu family of intersperesed repeats is comprised of ovr 500,000 members which may be divided into discrete subfamilies based upon mutations held in common between members. Distinct subfamilies of Alu sequences have amplified within the human genome in recent evolutionary history. Several individual Alu family members have amplified so recently in human evolution that they are variable as to presence and absence at specific loci within different human populations. Here, we report on the distribution of six polymorphic Alu insetions in a survey of 563 individuals from 14 human population groups across several continents. Our results indicate that these polymorphic Alu insertions probably have an African origin and that there is a much smaller amount of genetic variation between European populations than that found between other populations groups. Present address: Department of Pathology, Stanley S. Scott Cancer Center, Louisiana State University Medical Center, 1901 Perdido St., New Orleans, LA 70112 Correspondence to: M.A. Batzer     

Anthropological Genetics: Inferring the History of Our Species Through the Analysis of DNA

The genetic material, deoxyribonucleic acid (DNA), contains information about the evolutionary history of life. Both the relationships amongst organisms and the times of their divergence can be inferred from DNA sequences. Anthropological geneticists use DNA sequences to infer the evolutionary history of humans and their primate relatives. We review the basic methodology used to infer these relationships. We then review the anthropological genetic evidence for modern human origins. We conclude that modern humans evolved recently in Africa and then left to colonize the rest of the world within the last 50,000 years, largely replacing the other human groups that they encountered. Modern humans likely exchanged genes with Neanderthals prior to or early during their expansion out of Africa.     

Impacts of gene flow and logging history on the local genetic structure of a scattered tree species, Sorbus torminalis L. Crantz

Sorbus torminalis L. Crantz is a colonizing tree species usually found at low density in managed European forests. Using six microsatellite markers, we investigated spatial and temporal patterns of genetic structure within a 472-ha population of 185 individuals to infer processes shaping the distribution of genetic diversity. Only eight young stems were found to be the result of vegetative reproduction. Despite high levels of gene flow (standard deviation of gene dispersal = 360 m), marked patterns of isolation by distance were detected, associated with an aggregated distribution of individuals in approximately 100-m patches. This spatial structure of both genes and individuals is likely to result from patterns of seedling recruitment combined with low tree density. Our results suggest that landscape factors and logging cycles markedly shape the distribution of favourable sites for seedling establishment, which are then colonized by sibling cohorts as a result of joint seed transportation by frugivores. These combined genetic and demographic processes result in similar genetic structure both within and among logging units. However, conversion to high forest may enhance genetic structuring.     

Variable patterns of introgression in two sculpin hybrid zones suggest that genomic isolation differs among populations

Theory predicts that reproductive isolation may be due to intrinsic genetic incompatibilities or extrinsic ecological factors. Therefore, an understanding of the genetic basis of isolation may require analyses of evolutionary processes in situ to include environmental factors. Here we study genetic isolation between populations of sculpins ( Cottus ) at 168 microsatellites. Genomic clines were fit using 480 individuals sampled across independent natural hybrid zones that have formed between one invading species and two separate populations of a resident species. Our analysis tests for deviations from neutral patterns of introgression at individual loci based on expectations given genome-wide admixture. Roughly 51% of the loci analysed displayed significant deviations. An overall deficit of interspecific heterozygotes in 26% and 21% of the loci suggests that widespread underdominance drives genomic isolation. At the same time, selection promotes introgression of almost 30% of the markers, which implies that hybridization may increase the fitness of admixed individuals. Cases of overdominance or epistatic interactions were relatively rare. Despite the similarity of the two hybrid zones in their overall genomic composition, patterns observed at individual loci show little correlation between zones and many fit different genotypic models of fitness. At this point, it remains difficult to determine whether these results are due to differences in external selection pressures or cryptic genetic differentiation of distinct parental populations. In the future, data from mapped genetic markers and on variation of ecological factors will provide additional insights into the contribution of these factors to variation in the evolutionary consequences of hybridization.     

A molecular phylogeny of apple snails (Gastropoda, Caenogastropoda, Ampullariidae) with an emphasis on African species

Ampullariids are widespread in Africa, Asia, South- and Central America, and the Caribbean Islands. Basal phylogenetic relationships of the African genera Afropomus and Saulea have been inferred based on anatomical evidence. Until recently the Viviparidae was regarded as the sister-group of Ampullariidae, but recent molecular data infer a sister-group relationship with Campanilidae. We have used members of both families as outgroups in the present investigation on ampullariid phylogeny. We have used data from portions of five molecular loci, that is, the nuclear genes 18S rRNA, 28S rRNA and H3, and the mitochondrial genes 16S rRNA and COI. Our data most often infer a basal position of Afropomus . The West African species Saulea is inferred as the basal member of a clade including the South American Marisa and Pomacea . We hypothesize that evolutionary lineages leading to Saulea and the American genera were isolated from each other by vicariance events (Gondwanaland break-up 130–110 Mya). Our individual gene analyses inferred two major clades of the African Lanistes . However, in some analyses they were not inferred as sister-groups making Lanistes paraphyletic. The African and Asian genus Pila is most often inferred to be monophyletic (except for the generally unresolved 28S). Our analyses most often inferred a sister-group relationship between Lanistes and Pila . The very low genetic diversity of the endemic radiation of Lanistes in Lake Malawi suggests that the morphological divergence has happened much faster than the molecular divergence as is also evidenced from the cichlid radiations.     

Combining genetic structure and demographic analyses to estimate persistence in endangered Key deer (<Emphasis Type="Italic">Odocoileus virginianus clavium</Emphasis>)

Recent improvements in genetic analyses have paved the way in using molecular data to answer questions regarding evolutionary history, genetic structure, and demography. Key deer are a federally endangered subspecies assumed to be genetically unique, homogeneous, and have a female-biased population of approximately 900 deer. We used 985 bp of the mitochondrial cytochrome b gene and 12 microsatellite loci to test two hypotheses: (1) that Key deer are isolated and have reduced diversity compared to mainland deer and (2) that isolation of the Florida Keys has led to a small population size and a high risk of extinction. Our results indicate that Key deer are indeed genetically isolated from mainland white-tailed deer and that there is a lack of genetic substructure between islands. While Key deer exhibit reduced levels of genetic diversity compared to their mainland counterparts, they contain enough diversity to uniquely identify individual deer. Based on genetic identification, we estimated a census size of around 1000 individuals with a heavily skewed female-biased adult sex ratio. Furthermore, our genetic and contemporary demographic data were used to generate a species persistence model of the Key deer. Sensitivity tests within the population viability analysis brought to light the importance of fetal sex ratio and female survival as the primary factors at risk of driving the subspecies to extinction. This study serves as a prime example of how persistence models can be used to evaluate population viability in natural populations of endangered organisms.     

Testing local-scale panmixia provides insights into the cryptic ecology, evolution, and epidemiology of metazoan animal parasites

When every individual has an equal chance of mating with other individuals, the population is classified as panmictic. Amongst metazoan parasites of animals, local-scale panmixia can be disrupted due to not only non-random mating, but also non-random transmission among individual hosts of a single host population or non-random transmission among sympatric host species. Population genetics theory and analyses can be used to test the null hypothesis of panmixia and thus, allow one to draw inferences about parasite population dynamics that are difficult to observe directly. We provide an outline that addresses 3 tiered questions when testing parasite panmixia on local scales: is there greater than 1 parasite population/species, is there genetic subdivision amongst infrapopulations within a host population, and is there asexual reproduction or a non-random mating system? In this review, we highlight the evolutionary significance of non-panmixia on local scales and the genetic patterns that have been used to identify the different factors that may cause or explain deviations from panmixia on a local scale. We also discuss how tests of local-scale panmixia can provide a means to infer parasite population dynamics and epidemiology of medically relevant parasites.     

Genome-wide Insights into the Patterns and Determinants of Fine-Scale Population Structure in Humans

Studying genomic patterns of human population structure provides important insights into human evolutionary history and the relationship among populations, and it has significant practical implications for disease-gene mapping. Here we describe a principal component (PC)-based approach to studying intracontinental population structure in humans, identify the underlying markers mediating the observed patterns of fine-scale population structure, and infer the predominating evolutionary forces shaping local population structure. We applied this methodology to a data set of 650K SNPs genotyped in 944 unrelated individuals from 52 populations and demonstrate that, although typical PC analyses focus on the top axes of variation, substantial information about population structure is contained in lower-ranked PCs. We identified 18 significant PCs, some of which distinguish individual populations. In addition to visually representing sample clusters in PC biplots, we estimated the set of all SNPs significantly correlated with each of the most informative axes of variation. These polymorphisms, unlike ancestry-informative markers (AIMs), constitute a much larger set of loci that drive genomic signatures of population structure. The genome-wide distribution of these significantly correlated markers can largely be accounted for by the stochastic effects of genetic drift, although significant clustering does occur in genomic regions that have been previously implicated as targets of recent adaptive evolution.     

HGDP and HapMap Analysis by Ancestry Mapper Reveals Local and Global Population Relationships

Knowledge of human origins, migrations, and expansions is greatly enhanced by the availability of large datasets of genetic information from different populations and by the development of bioinformatic tools used to analyze the data. We present Ancestry Mapper, which we believe improves on existing methods, for the assignment of genetic ancestry to an individual and to study the relationships between local and global populations. The principle function of the method, named Ancestry Mapper, is to give each individual analyzed a genetic identifier, made up of just 51 genetic coordinates, that corresponds to its relationship to the HGDP reference population. As a consequence, the Ancestry Mapper Id (AMid) has intrinsic biological meaning and provides a tool to measure similarity between world populations. We applied Ancestry Mapper to a dataset comprised of the HGDP and HapMap data. The results show distinctions at the continental level, while simultaneously giving details at the population level. We clustered AMids of HGDP/HapMap and observe a recapitulation of human migrations: for a small number of clusters, individuals are grouped according to continental origins; for a larger number of clusters, regional and population distinctions are evident. Calculating distances between AMids allows us to infer ancestry. The number of coordinates is expandable, increasing the power of Ancestry Mapper. An R package called Ancestry Mapper is available to apply this method to any high density genomic data set.     

Integrating Pool-seq uncertainties into demographic inference

Next-generation sequencing of pooled samples (Pool-seq) is a popular method to assess genome-wide diversity patterns in natural and experimental populations. However, Pool-seq is associated with specific sources of noise, such as unequal individual contributions. Consequently, using Pool-seq for the reconstruction of evolutionary history has remained underexplored. Here we describe a novel Approximate Bayesian Computation (ABC) method to infer demographic history, explicitly modelling Pool-seq sources of error. By jointly modelling Pool-seq data, demographic history and the effects of selection due to barrier loci, we obtain estimates of demographic history parameters accounting for technical errors associated with Pool-seq. Our ABC approach is computationally efficient as it relies on simulating subsets of loci (rather than the whole-genome) and on using relative summary statistics and relative model parameters. Our simulation study results indicate Pool-seq data allows distinction between general scenarios of ecotype formation (single versus parallel origin) and to infer relevant demographic parameters (e.g. effective sizes and split times). We exemplify the application of our method to Pool-seq data from the rocky-shore gastropod Littorina saxatilis, sampled on a narrow geographical scale at two Swedish locations where two ecotypes (Wave and Crab) are found. Our model choice and parameter estimates show that ecotypes formed before colonization of the two locations (i.e. single origin) and are maintained despite gene flow. These results indicate that demographic modelling and inference can be successful based on pool-sequencing using ABC, contributing to the development of suitable null models that allow for a better understanding of the genetic basis of divergent adaptation.     

Forest refugia and riverine barriers promote diversification in the West African pygmy shrew (Crocidura obscurior complex,Soricomorpha)

The Crocidura obscurior or West African pygmy shrew complex is endemic to West African forests from south‐eastern Guinea, eastern Liberia, southern Côte d'Ivoire and south‐western Ghana. We explore the genetic and morphometric diversity of 239 individuals of the C. obscurior complex from 17 localities across its geographical range. Using genetic data from three mitochondrial (16S, cytochrome b and COI) and four nuclear markers (BRCA1, STAT5A, HDAC2 and RIOK3) and skull geometric morphometrics, we show that this complex is composed of two cryptic and sympatric species, C. obscurior and C. eburnea. We then test several hypotheses to infer their evolutionary history. The observed phylogeographical pattern based on cytochrome b and COI sequences fits the forest refuge theory: during arid phases of the Plio‐Pleistocene, around 3.5, 2.1, 1 and 0.5 Mya, a small number of populations survived in isolated forest patches and diverged allopatrically. During wetter climatic periods, forests expanded, leading to secondary contacts between previously isolated populations. Our results also suggest the possible contribution of episodes of isolation in subrefuges. Historical variation of the West African hydrographic network could also have contributed to the observed patterns of genetic differentiation. Rivers such as the Volta and Sassandra may act as past and/or current barriers to gene flow. Although these two species have sympatric distributions, their phylogeographical histories are somewhat dissimilar due to small differences in their dispersal abilities and ecological requirements.