Search

Creole sheep represent a strategic genetic resource for populations living in marginal areas under financial restrictions on the American continent. Six Colombian sheep breeds (two wool (BCL-Boyacá and NCL-Nariño, 12 and 14 samples) and four hair (OPCE-Ethiopian, 54 samples; OPCS-Sudan, 74 samples; OPCP-Pelibeuy, 59 samples; OPCW-Wayúu, 24 samples) were genotyped using the Illumina Ovine SNP50 BeadChip. Data was also included from international 44 breeds from International Sheep Genomics Consortium (ISGC) and from data published in previous a previous work on the Caribbean and African breeds. Although geographically separated, wool (NCL, BCL) and hair types (OPCE, OPCS, OPCW) presented little genetic differentiation (FST 0.05) at a global level but several groups of animals separated suggesting local clustering due to geographical isolation. The OPCP underwent a recent crossing with Mexican Pelibuey, explaining its differentiation. Findings in this work such as the proximity to West African Djallonké (WAD) and Barbados Black Belly (BBB), suggest different introductions of African type animals from the Caribbean region on a pre-existing genetic basis formed by animals deriving from the first importations coming from Europe in colonial times. As expected, Colombian wool breeds showed, in particular in Admixture software results, a greater genomic component in common with European breeds and in particular with Iberian ones (Churra). This study provides a basis for future research into the genetic diversity within and between the Colombian sheep breeds analysed, and scientific data for policy decisions on Farm Animal Genetic Resources (FAnGR).

Demography determines the strength of genetic drift, which generally reduces genetic variation and the efficacy of selection. Here, we disentangled the importance of demographic processes at a local scale (census size and mating system) and at a species-range scale (old split between population clusters, recolonization after the last glaciation cycle, and admixture) in determining within-population genomic diversity and genomic signatures of positive selection. Analyses were based on re-sequence data from 52 populations of North American Arabidopsis lyrata collected across its entire distribution. The mating system and range dynamics since the last glaciation cycle explained around 60% of the variation in genomic diversity among populations and 52% of the variation in the signature of positive selection. Diversity was lowest in selfing compared with outcrossing populations and in areas further away from glacial refugia. In parallel, reduced positive selection was found in selfing populations and in populations with a longer route of postglacial range expansion. The signature of positive selection was also reduced in populations without admixture. We conclude that recent range expansion can have a profound influence on diversity in coding and non-coding DNA, similar in magnitude to the shift toward selfing. Distribution limits may in fact be caused by reduced effective population size and compromised positive selection in recently colonized parts of the range., (© 2020 The Author(s).)

Hetaerina americana Fabricius, 1798 has a long and irresolute taxonomic history. Several synonyms have been suggested (H. californica Hagen in Selys-Longchamps, 1859, H. basalis Hagen in Selys-Longchamps, 1859, H. texana Walsh, 1863, H. scelerata Walsh, 1863, etc.), related to the variation in the size of the wing spots as well as to the morphology of the male cerci. However, Calvert (1901) suggested that H. americana represents one variable species. Nevertheless, Vega-Sánchez et al. (2019) through a genetic and morphological analysis presented evidence to propose that H. americana represents a species complex. In the present work, we describe a new species that belongs to this complex: H. calverti sp. nov. The morphological characteristics by which males and females of H. calverti differ from H. americana are highlighted. The most important character for the differentiation of males is the shape of the cerci and the size of the individuals (when the two species are in sympatry). In females, the main differences are in the shape of the intersternites and the medio-dorsal carina of the last segment of the abdomen. Some generalities about the biology of the species are presented, including geographical distribution patterns and genetic divergence data. [urn:lsid:zoobank.org:act:F5C329AE-7A00-4979-8A0D-A13D869E54B1].

Background: Emerging influenza A/H3N2 clades have been associated with M1 gene mutations which affect the performance of commercial PCR assays., Objectives and Study Design: The evolution and prevalence of problematic M1 mutations, and their associated viral clades, were investigated. All European and USA isolates from the GISAID database with both HA and M1 sequences available, collected during the respiratory seasons from the Fall of 2007 through January of 2018, were analyzed., Results: Five M1 target region patterns, designated A-E, were observed in more than 10% of the isolates during a season, with patterns that appeared sequentially, each having one additional mutation. The C153T mutation was universal. Pattern A, which only had the single mutation, predominated between 2007/08 and 2009/10. Dual- and triple-mutation patterns (B and C) emerged in 2010/11 and 2011/12 respectively, and pattern C predominated for one season (2012/13). In 2012/13, the problematic quadruple-mutation containing C163T first appeared in 3C.2 viruses. Seasons 2013/14 and 14/15 were associated with significant viral diversity with five clades and four M1 patterns co-circulating, with different rates in Europe and the USA. Since 2014, clade 3C.2a with M1 pattern D has emerged as the predominant type. During 2016/17 season, a new quintuplet mutation pattern (E) emerged in cluster 3C.2a1 isolates., Conclusions: M1 target region mutations have been prevalent for more than ten years, with the number of mutations continually increasing. Often population inferences of M1 mutations can be made based on viral clade. However, gene segment reassortment can affect predictive abilities., (Copyright © 2018 Elsevier B.V. All rights reserved.)

The USA300 North American epidemic (USA300-NAE) clone of methicillin-resistant Staphylococcus aureus has caused a wave of severe skin and soft tissue infections in the United States since it emerged in the early 2000s, but its geographic origin is obscure. Here we use the population genomic signatures expected from the serial founder effects of a geographic range expansion to infer the origin of USA300-NAE and identify polymorphisms associated with its spread. Genome sequences from 357 isolates from 22 U.S. states and territories and seven other countries are compared. We observe two significant signatures of range expansion, including decreases in genetic diversity and increases in derived allele frequency with geographic distance from the Pennsylvania region. These signatures account for approximately half of the core nucleotide variation of this clone, occur genome wide, and are robust to heterogeneity in temporal sampling of isolates, human population density, and recombination detection methods. The potential for positive selection of a gyrA fluoroquinolone resistance allele and several intergenic regions, along with a 2.4 times higher recombination rate in a resistant subclade, is noted. These results are the first to show a pattern of genetic variation that is consistent with a range expansion of an epidemic bacterial clone, and they highlight a rarely considered but potentially common mechanism by which genetic drift may profoundly influence bacterial genetic variation. IMPORTANCE The process of geographic spread of an origin population by a series of smaller populations can result in distinctive patterns of genetic variation. We detect these patterns for the first time with an epidemic bacterial clone and use them to uncover the clone's geographic origin and variants associated with its spread. We study the USA300 clone of methicillin-resistant Staphylococcus aureus , which was first noticed in the early 2000s and subsequently became the leading cause of skin and soft tissue infections in the United States. The eastern United States is the most likely origin of epidemic USA300. Relatively few variants, which include an antibiotic resistance mutation, have persisted during this clone's spread. Our study suggests that an early chapter in the genetic history of this epidemic bacterial clone was greatly influenced by random subsampling of isolates during the clone's geographic spread., (Copyright © 2018 Challagundla et al.)

Biological invasions rank among the main global threats for biodiversity. The Eastern mosquitofish (Gambusia holbrooki) is considered one of the 100 world worst invasive species due to its high adaptation capability to new environments. Using the restriction-site-associated DNA tags (RADtags), introduced European locations were compared against native US mosquitofish populations to analyse genomic changes that occurred during invasive process of European locations. After filtering, 7724 RADtags containing only one SNP were retained for population studies. Comparative genomics indicated that 186 of these RADtags matched sequences in the transcriptome of Xyphophorus maculatus, the most closely related genome available. Genomic analyses showed that invasive populations show high reductions in diversity. Further, analyses of population structuring based on these data are concordant with previous analyses based on microsatellites. It is concluded that during the invasion process genetic drift was the main evolutionary force affecting patterns of diversity and population structure. While recognizing that positive selection could be masked by the strong drift during founder events, adaptive processes were evidenced in a reduced number of RADtags (<2%), with only one of these in a putative coding region. Surprisingly, balancing selection was detected in several coding RADtags, suggesting that the preservation of polymorphism in specific genes could be more important than the average population diversity for the population maintenance at any location, particularly for the survival of introduced populations., (© 2016 John Wiley & Sons Ltd.)

Sclerotinia sclerotiorum is an important pathogen of numerous crops in the North Central region of the United States. The objective of this study was to examine the genetic diversity of 145 isolates of the pathogen from multiple hosts in the region. Mycelial compatibility groups (MCG) and microsatellite haplotypes were determined and analyzed for standard estimates of population genetic diversity and the importance of host and distance for genetic variation was examined. MCG tests indicated there were 49 different MCGs in the population and 52 unique microsatellite haplotypes were identified. There was an association between MCG and haplotype such that isolates belonging to the same MCG either shared identical haplotypes or differed at no more than 2 of the 12 polymorphic loci. For the majority of isolates, there was a one-to-one correspondence between MCG and haplotype. Eleven MCGs shared haplotypes. A single haplotype was found to be prevalent throughout the region. The majority of genetic variation in the isolate collection was found within rather than among host crops, suggesting little genetic divergence of S. sclerotiorum among hosts. There was only weak evidence of isolation by distance. Pairwise population comparisons among isolates from canola, dry bean, soybean and sunflower suggested that gene flow between host-populations is more common for some crops than others. Analysis of linkage disequilibrium in the isolates from the four major crops indicated primarily clonal reproduction, but also evidence of genetic recombination for isolates from canola and sunflower. Accordingly, genetic diversity was highest for populations from canola and sunflower. Distribution of microsatellite haplotypes across the study region strongly suggest that specific haplotypes of S. sclerotiorum are often found on multiple crops, movement of individual haplotypes among crops is common and host identity is not a barrier to gene flow for S. sclerotiorum in the north central United States.

Studies of species with continental distributions continue to identify intraspecific lineages despite continuous habitat. Lineages may form due to isolation by distance, adaptation, divergence across barriers, or genetic drift following range expansion. We investigated lineage diversification and admixture within American black bears (Ursus americanus) across their range using 22 k single nucleotide polymorphisms and mitochondrial DNA sequences. We identified three subcontinental nuclear clusters which we further divided into nine geographic regions: Alaskan (Alaska-East), eastern (Central Interior Highlands, Great Lakes, Northeast, Southeast), and western (Alaska-West, West, Pacific Coast, Southwest). We estimated that the western cluster diverged 67 ka, before eastern and Alaskan divergence 31 ka; these divergence dates contrasted with those from the mitochondrial genome where clades A and B diverged 1.07 Ma, and clades A-east and A-west diverged 169 ka. We combined estimates of divergence timing with hindcast species distribution models to infer glacial refugia for the species in Beringia, Pacific Northwest, Southwest, and Southeast. Our results show a complex arrangement of admixture due to expansion out of multiple refugia. The delineation of the genomic population clusters was inconsistent with the ranges for 16 previously described subspecies. Ranges for U. a. pugnax and U. a. cinnamomum were concordant with admixed clusters, calling into question how to order taxa below the species level. Additionally, our finding that U. a. floridanus has not diverged from U. a. americanus also suggests that morphology and genetics should be reanalyzed to assess taxonomic designations relevant to the conservation management of the species., (© The Author 2015. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution. All rights reserved.For permissions please email: journals.permissions@oup.com.)

Animal models are crucial components in the search for better understanding of the biological basis of psychiatric disorders and for the development of novel drugs. Research, in general, and research with animal models, in particular, relies on the consistency of effects of investigated drugs or manipulations across experiments. In that context, it had been noted that behavioral responses to lithium in ICR (CD-1) mice from Harlan Israel have changed across the last years. To examine this change, the present study compared the effect of lithium treatment in ICR mice from Harlan Israel with the ICR mice from Harlan USA. The mice were treated with chronic oral lithium. Their lithium serum levels were measured and their behavior in the forced swim test (FST) was evaluated. The mice were also treated with [(3)H]-inositol ICV and lithium injection and their frontal cortex [(3)H]-phosphoinositols accumulation was measured. Results show that lithium serum levels in Israeli mice were significantly lower compared with the USA mice, that lithium had no behavioral effect in the Israeli mice but significantly reduced FST immobility time of the USA mice, and that phosphoinositols accumulation was much more strongly affected by lithium in the USA mice compared with the Israeli mice. These results suggest that the Israeli Harlan colony of ICR mice changed significantly from the original ICR colony in Harlan USA and that the differences might be related to absorption or secretion of lithium., (Copyright © 2014 Elsevier Inc. All rights reserved.)

Dramatic changes in the North American landscape over the last 12 000 years have shaped the genomes of the small mammals, such as the white-footed mouse (Peromyscus leucopus), which currently inhabit the region. However, very recent interactions of populations with each other and the environment are expected to leave the most pronounced signature on rapidly evolving nuclear microsatellite loci. We analyzed landscape characteristics and microsatellite markers of P. leucopus populations along a transect from southern Ohio to northern Michigan, in order to evaluate hypotheses about the spatial distribution of genetic heterogeneity. Genetic diversity increased to the north and was best approximated by a single-variable model based on habitat availability within a 0.5-km radius of trapping sites. Interpopulation differentiation measured by clustering analysis was highly variable and not significantly related to latitude or habitat availability. Interpopulation differentiation measured as FST values and chord distance was correlated with the proportion of habitat intervening, but was best explained by agricultural distance and by latitude. The observed gradients in diversity and interpopulation differentiation were consistent with recent habitat availability being the major constraint on effective population size in this system, and contradicted the predictions of both the postglacial expansion and core-periphery hypotheses.

Bottlenecks, founder events, and genetic drift often result in decreased genetic diversity and increased population differentiation. These events may follow abundance declines due to natural or anthropogenic perturbations, where translocations may be an effective conservation strategy to increase population size. American black bears (Ursus americanus) were nearly extirpated from the Central Interior Highlands, USA by 1920. In an effort to restore bears, 254 individuals were translocated from Minnesota, USA, and Manitoba, Canada, into the Ouachita and Ozark Mountains from 1958 to 1968. Using 15 microsatellites and mitochondrial haplotypes, we observed contemporary genetic diversity and differentiation between the source and supplemented populations. We inferred four genetic clusters: Source, Ouachitas, Ozarks, and a cluster in Missouri where no individuals were translocated. Coalescent models using approximate Bayesian computation identified an admixture model as having the highest posterior probability (0.942) over models where the translocation was unsuccessful or acted as a founder event. Nuclear genetic diversity was highest in the source (AR = 9.11) and significantly lower in the translocated populations (AR = 7.07-7.34; P = 0.004). The Missouri cluster had the lowest genetic diversity (AR = 5.48) and served as a natural experiment showing the utility of translocations to increase genetic diversity following demographic bottlenecks. Differentiation was greater between the two admixed populations than either compared to the source, suggesting that genetic drift acted strongly over the eight generations since the translocation. The Ouachitas and Missouri were previously hypothesized to be remnant lineages. We observed a pretranslocation remnant signature in Missouri but not in the Ouachitas., (Published 2014. This article is a U.S. Government work and is in the public domain in the USA.)

Genetic drift (GD) randomly impacts small breeds and imported populations. Therefore, it can impact policies that affect conservation of animal genetic resources. This paper evaluates GD for a population of Meishan pigs imported into the United States and explores the ramifications of GD on access and benefit sharing of genetic resources under the Nagoya Protocol (NP) of the United Nations' Convention on Biological Diversity. The NP was motivated by concerns about fair and equitable benefit sharing of genetic resources across life forms. In this experiment, 35 microsatellite markers were used to quantify the level of GD that occurred between Meishan pigs (Meishan-China; n = 22) imported into the United States in the late 1980s and contemporary Meishan (Meishan-US; n = 42), which have been randomly bred since importation. The Meishan-US consisted of 2 subpopulations (Meishan-MARC and Meishan-ISU). Five other breeds were also included in the analysis to serve as reference populations: Fengjing and Minzhu, which were imported with Meishan-China, and Duroc, Berkshire, and Yorkshire from the United States. Mean shift in allele frequency was 0.11 (SE = 0.019) due to GD for Meishan-US vs. Meishan-China with some loci having changed allele frequencies by greater than 0.20. Principle coordinate analysis confirmed divergence among the Meishan populations. Model-based clustering tended to place the U.S. and Chinese breeds into 2 distinct clusters, likely due to differences in allele frequencies between U.S. and Chinese breeds. Contemporary Meishan-US has become differentiated from the original imported animals due to GD. Attributing future performance of Meishan-US to Meishan-China, as set forth by NP, is problematic due to GD. As an imported breed becomes established there will be an increasing number of breeders who may have different selection goals and private treaty contracts will govern the exchange of stock between them. Therefore, considering biological phenomena and social dynamics simultaneously draws into question whether a rigorous access and benefit sharing protocol as envisioned in the NP will be operational.

Two distinct antigenic clusters were previously identified among the H3N2 swine influenza A viruses (IAVs) and were designated H3N2SIV-alpha and H3N2SIV-beta (Feng et al., 2013. Journal of Virology 87 (13), 7655-7667). A consistent mutation was observed at the position 189 of hemagglutinin (R189K) between H3N2SIV-alpha and H3N2SIV-beta fair isolates. To evaluate the contribution of R189K mutation to the antigenic drift from H3N2SIV-alpha to H3N2SIV-beta, four reassortant viruses with 189R or 189K were generated. The antigenic cartography demonstrated that the R189K mutation in the hemagglutinin of H3N2 IAV contributed to the antigenic drift, separating these viruses into H3N2SIV-alpha to H3N2SIV-beta. This R189K mutation was also found to contribute to the cross-reaction with several ferret sera raised against historical human IAVs with hemagglutinin carrying 189K. This study suggests that the R189K mutation plays a vital role in the antigenicity of swine and human H3N2 IAVs and identification of this antigenic determinant will help us rapidly identify antigenic variants in influenza surveillance., (© 2013 Elsevier Inc. All rights reserved.)

The potential for widespread, mobile species to exhibit genetic structure without clear geographic barriers is a topic of growing interest. Yet the patterns and mechanisms of structure--particularly over broad spatial scales--remain largely unexplored for these species. Bobcats occur across North America and possess many characteristics expected to promote gene flow. To test whether historical, topographic or ecological factors have influenced genetic differentiation in this species, we analysed 1 kb mtDNA sequence and 15 microsatellite loci from over 1700 samples collected across its range. The primary signature in both marker types involved a longitudinal cline with a sharp transition, or suture zone, occurring along the Great Plains. Thus, the data distinguished bobcats in the eastern USA from those in the western half, with no obvious physical barrier to gene flow. Demographic analyses supported a scenario of expansion from separate Pleistocene refugia, with the Great Plains representing a zone of secondary contact. Substructure within the two main lineages likely reflected founder effects, ecological factors, anthropogenic/topographic effects or a combination of these forces. Two prominent topographic features, the Mississippi River and Rocky Mountains, were not supported as significant genetic barriers. Ecological regions and environmental correlates explained a small but significant proportion of genetic variation. Overall, results implicate historical processes as the primary cause of broad-scale genetic differentiation, but contemporary forces seem to also play a role in promoting and maintaining structure. Despite the bobcat's mobility and broad niche, large-scale landscape changes have contributed to significant and complex patterns of genetic structure., (© 2012 Blackwell Publishing Ltd.)

Switchgrass (Panicum virgatum), a central and Eastern USA native, is highly valued as a component in tallgrass prairie and savanna restoration and conservation projects and a potential bioenergy feedstock. The purpose of this study was to identify regional diversity, gene pools, and centers-of-diversity of switchgrass to gain an understanding of its post-glacial evolution and to identify both the geographic range and potential overlap between functional gene pools. We sampled a total of 384 genotypes from 49 accessions that included the three main taxonomic groups of switchgrass (lowland 4x, upland 4x, and upland 8x) along with one accession possessing an intermediate phenotype. We identified primary centers of diversity for switchgrass in the eastern and western Gulf Coast regions. Migration, drift, and selection have led to adaptive radiation in switchgrass, creating regional gene pools within each of the main taxa. We estimate that both upland-lowland divergence and 4x-to-8x polyploidization within switchgrass began approximately 1.5-1 M ybp and that subsequent ice age cycles have resulted in gene flow between ecotype lineages and between ploidy levels. Gene flow has resulted in "hot spots" of genetic diversity in the southeastern USA and along the Atlantic Seaboard.

A rice mini-core collection consisting of 217 accessions has been developed to represent the USDA core and whole collections that include 1,794 and 18,709 accessions, respectively. To improve the efficiency of mining valuable genes and broadening the genetic diversity in breeding, genetic structure and diversity were analyzed using both genotypic (128 molecular markers) and phenotypic (14 numerical traits) data. This mini-core had 13.5 alleles per locus, which is the most among the reported germplasm collections of rice. Similarly, polymorphic information content (PIC) value was 0.71 in the mini-core which is the highest with one exception. The high genetic diversity in the mini-core suggests there is a good possibility of mining genes of interest and selecting parents which will improve food production and quality. A model-based clustering analysis resulted in lowland rice including three groups, aus (39 accessions), indica (71) and their admixtures (5), upland rice including temperate japonica (32), tropical japonica (40), aromatic (6) and their admixtures (12) and wild rice (12) including glaberrima and four other species of Oryza. Group differentiation was analyzed using both genotypic distance Fst from 128 molecular markers and phenotypic (Mahalanobis) distance D(2) from 14 traits. Both dendrograms built by Fst and D(2) reached similar-differentiative relationship among these genetic groups, and the correlation coefficient showed high value 0.85 between Fst matrix and D(2) matrix. The information of genetic and phenotypic differentiation could be helpful for the association mapping of genes of interest. Analysis of genotypic and phenotypic diversity based on genetic structure would facilitate parent selection for broadening genetic base of modern rice cultivars via breeding effort.

Although vaccines pose the best means of preventing influenza infection, strain selection and optimal implementation remain difficult due to antigenic drift and a lack of understanding global spread. Detecting viral movement by sequence analysis is complicated by skewed geographic and seasonal distributions in viral isolates. We propose a probabilistic method that accounts for sampling bias through spatiotemporal clustering and modeling regional and seasonal transmission as a binomial process. Analysis of H3N2 not only confirmed East-Southeast Asia as a source of new seasonal variants, but also increased the resolution of observed transmission to a country level. H1N1 data revealed similar viral spread from the tropics. Network analysis suggested China and Hong Kong as the origins of new seasonal H3N2 strains and the United States as a region where increased vaccination would maximally disrupt global spread of the virus. These techniques provide a promising methodology for the analysis of any seasonal virus, as well as for the continued surveillance of influenza.

Recently introduced, exotic plant pathogens may exhibit low genetic diversity and be limited to clonal reproduction. However, rapidly mutating molecular markers such as microsatellites can reveal genetic variation within these populations and be used to model putative migration patterns. Phytophthora ramorum is the exotic pathogen, discovered in the late 1990s, that is responsible for sudden oak death in California forests and ramorum blight of common ornamentals. The nursery trade has moved this pathogen from source populations on the West Coast to locations across the United States, thus risking introduction to other native forests. We examined the genetic diversity of P. ramorum in United States nurseries by microsatellite genotyping 279 isolates collected from 19 states between 2004 and 2007. Of the three known P. ramorum clonal lineages, the most common and genetically diverse lineage in the sample was NA1. Two eastward migration pathways were revealed in the clustering of NA1 isolates into two groups, one containing isolates from Connecticut, Oregon, and Washington and the other isolates from California and the remaining states. This finding is consistent with trace forward analyses conducted by the US Department of Agriculture's Animal and Plant Health Inspection Service. At the same time, genetic diversities in several states equaled those observed in California, Oregon, and Washington and two-thirds of multilocus genotypes exhibited limited geographic distributions, indicating that mutation was common during or subsequent to migration. Together, these data suggest that migration, rapid mutation, and genetic drift all play a role in structuring the genetic diversity of P. ramorum in US nurseries. This work demonstrates that fast-evolving genetic markers can be used to examine the evolutionary processes acting on recently introduced pathogens and to infer their putative migration patterns, thus showing promise for the application of forensics to plant pathogens.

In recent years approximate Bayesian computation (ABC) methods have become popular in population genetics as an alternative to full-likelihood methods to make inferences under complex demographic models. Most ABC methods rely on the choice of a set of summary statistics to extract information from the data. In this article we tested the use of the full allelic distribution directly in an ABC framework. Although the ABC techniques are becoming more widely used, there is still uncertainty over how they perform in comparison with full-likelihood methods. We thus conducted a simulation study and provide a detailed examination of ABC in comparison with full likelihood in the case of a model of admixture. This model assumes that two parental populations mixed at a certain time in the past, creating a hybrid population, and that the three populations then evolve under pure drift. Several aspects of ABC methodology were investigated, such as the effect of the distance metric chosen to measure the similarity between simulated and observed data sets. Results show that in general ABC provides good approximations to the posterior distributions obtained with the full-likelihood method. This suggests that it is possible to apply ABC using allele frequencies to make inferences in cases where it is difficult to select a set of suitable summary statistics and when the complexity of the model or the size of the data set makes it computationally prohibitive to use full-likelihood methods.

Large scale genomic surveys are partly motivated by the idea that the neutral genetic variation of a population may be used to reconstruct its migration history. However, our ability to trace back the colonization pathways of a species from their genetic footprints is limited by our understanding of the genetic consequences of a range expansion. Here, we study, by means of simulations and analytical methods, the neutral dynamics of gene frequencies in an asexual population undergoing a continual range expansion in one dimension. During such a colonization period, lineages can fix at the wave front by means of a "surfing" mechanism [Edmonds, C.A., Lillie, A.S., Cavalli-Sforza, L.L., 2004. Mutations arising in the wave front of an expanding population. Proc. Natl. Acad. Sci. 101, 975-979]. We quantify this phenomenon in terms of (i) the spatial distribution of lineages that reach fixation and, closely related, (ii) the continual loss of genetic diversity (heterozygosity) at the wave front, characterizing the approach to fixation. Our stochastic simulations show that an effective population size can be assigned to the wave that controls the (observable) gradient in heterozygosity left behind the colonization process. This effective population size is markedly higher in the presence of cooperation between individuals ("pushed waves") than when individuals proliferate independently ("pulled waves"), and increases only sub-linearly with deme size. To explain these and other findings, we develop a versatile analytical approach, based on the physics of reaction-diffusion systems, that yields simple predictions for any deterministic population dynamics. Our analytical theory compares well with the simulation results for pushed waves, but is less accurate in the case of pulled waves when stochastic fluctuations in the tip of the wave are important.

If an influenza pandemic struck today, borders might close, the global economy would be severely impacted, international vaccine supplies and health are systems would be overwhelmed, and some people might panic. To limit the fallout, the industrialized world must create a detailed response strategy involving the public and private sectors. Some experts feel we are overdue for a flu pandemic and the SARS pandemic of 2003 could have been the wake up call to begin preparations. Fortunately there is some assistance coming from the federal government. On January 12, U.S. Department of Health & Human Services Secretary Mike Leavitt announced funding to assist in the preparation for a pandemic flu response. $100 million is being provided initially with another $250 million due later this year to assist states in pandemic flu preparedness. Texas' initial allocation is $5,875,044. While some believe that the AI (H5N1) causing illness and deaths in Asia and Turkey will be the pandemic flu strain, there is no guarantee that will occur. Thus, without knowing which strain may lead to a pandemic, development and manufacturing of a vaccine is delayed.

A simple model of random copying among individuals, similar to the population genetic model of random drift, can predict the variability in the popularity of cultural variants. Here, we show that random drift also explains a biologically relevant cultural phenomenon--changes in the distributions of popularity of dog breeds in the United States in each of the past 50 years. There are, however, interesting deviations from the model that involve large changes in the popularity of certain breeds. By identifying meaningful departures from our null model, we show how it can serve as a foundation for studying culture change quantitatively, using the tools of population genetics.

Eastern North America is one of at least six regions of the world where agriculture is thought to have arisen wholly independently. The primary evidence for this hypothesis derives from morphological changes in the archaeobotanical record of three important crops--squash, goosefoot and sunflower--as well as an extinct minor cultigen, sumpweed. However, the geographical origins of two of the three primary domesticates--squash and goosefoot--are now debated, and until recently sunflower (Helianthus annuus L.) has been considered the only undisputed eastern North American domesticate. The discovery of 4,000-year-old domesticated sunflower remains from San Andrés, Tabasco, implies an earlier and possibly independent origin of domestication in Mexico and has stimulated a re-examination of the geographical origin of domesticated sunflower. Here we describe the genetic relationships and pattern of genetic drift between extant domesticated strains and wild populations collected from throughout the USA and Mexico. We show that extant domesticates arose in eastern North America, with a substantial genetic bottleneck occurring during domestication.

It has been proposed that free-living microorganisms exhibit ubiquitous dispersal, do not form geographically isolated populations and rarely (if ever) speciate via allopatry. We studied island-like hot spring cyanobacterial communities in which geographical isolation should be prominent and detectable if it influences the evolution of bacteria. The genetic diversity of cyanobacteria indigenous to North American, Japanese, New Zealand and Italian springs was surveyed by (i) amplification and cloning of 16S rRNA and 16S-23S internal transcribed spacer regions; (ii) lineage-specific oligonucleotide probing (used to verify the predominance of cloned sequences); and (iii) lineage-specific polymerase chain reaction (PCR) (used to search for possible rare genotypes). Phylogenetic and distribution patterns were found to be consistent with the occurrence of geographical isolation at both global and local spatial scales, although different cyanobacterial lineages were found to vary in their distribution. A lack of correspondence between biological patterning and the chemical character of springs sampled suggested that the geographical distribution of thermophilic cyanobacteria cannot be explained by the 20 potential niche-determining chemical parameters that we assayed. Thus, geographical isolation (i.e. genetic drift) must in part be responsible for driving the observed evolutionary divergences. Geographical isolation may be an important underestimated aspect of microbial evolution.

Nested cladistic analysis (NCA) is increasingly being used to infer historical population-level processes, including population fragmentation, range expansion and long-distance colonization. However, the effects on interpretation of NCA inferences of stochastic extinction of haplotypes due to genetic drift (lineage sorting), or of haplotype loss via localized biotic or climatic influences, have not been thoroughly explored. We provide empirical evidence suggesting that NCA may misinterpret population history when haplotypes or haplotype groups from one clade are replaced by those of another clade. We do so by using NCA to analyse mitochondrial sequences from the toad Bufo woodhousii from 45 locations spanning the Great Plains and southwestern USA. Portions of this region were glaciated and/or desertified in the late Pleistocene and early Holocene, and hence uninhabitable for plains-dwelling organisms. Although NCA inferences of isolation-by-distance and gradual range expansion in B. woodhousii are compatible with expectations based on climatic data and toad biology, NCA also detected several instances of long-distance movement. Such movement seems unlikely, given the low vagility of this species. We conclude that inferences of long-distance colonization likely result from extinction of haplotypes in intervening areas. We suggest using additional methods to look for congruent inferences, and amending the NCA inference key, to help avoid misinterpretations resulting from haplotype extinction.

Using properties of moment stationarity we develop exact expressions for the mean and covariance of allele frequencies at a single locus for a set of populations subject to drift, mutation, and migration. Some general results can be obtained even for arbitrary mutation and migration matrices, for example: (1) Under quite general conditions, the mean vector depends only on mutation rates, not on migration rates or the number of populations. (2) Allele frequencies covary among all pairs of populations connected by migration. As a result, the drift, mutation, migration process is not ergodic when any finite number of populations is exchanging genes. In addition, we provide closed-form expressions for the mean and covariance of allele frequencies in Wright's finite-island model of migration under several simple models of mutation, and we show that the correlation in allele frequencies among populations can be very large for realistic rates of mutation unless an enormous number of populations are exchanging genes. As a result, the traditional diffusion approximation provides a poor approximation of the stationary distribution of allele frequencies among populations. Finally, we discuss some implications of our results for measures of population structure based on Wright's F-statistics.

The presence of four lysosomal storage diseases (LSDs) at increased frequency in the Ashkenazi Jewish population has suggested to many the operation of natural selection (carrier advantage) as the driving force. We compare LSDs and nonlysosomal storage diseases (NLSDs) in terms of the number of mutations, allele-frequency distributions, and estimated coalescence dates of mutations. We also provide new data on the European geographic distribution, in the Ashkenazi population, of seven LSD and seven NLSD mutations. No differences in any of the distributions were observed between LSDs and NLSDs. Furthermore, no regular pattern of geographic distribution was observed for LSD versus NLSD mutations-with some being more common in central Europe and others being more common in eastern Europe, within each group. The most striking disparate pattern was the geographic distribution of the two primary Tay-Sachs disease mutations, with the first being more common in central Europe (and likely older) and the second being exclusive to eastern Europe (primarily Lithuania and Russia) (and likely much younger). The latter demonstrates a pattern similar to two other recently arisen Lithuanian mutations, those for torsion dystonia and familial hypercholesterolemia. These observations provide compelling support for random genetic drift (chance founder effects, one approximately 11 centuries ago that affected all Ashkenazim and another approximately 5 centuries ago that affected Lithuanians), rather than selection, as the primary determinant of disease mutations in the Ashkenazi population.