Determining the source of individuals: multilocus genotyping in nonequilibrium population genetics

Integrated Discriminant Evaluation of Molecular Genetic Markers and Genetic Diversity Parameters of Endangered Balearic Dog Breeds

The genetic diversity analysis of six dog breeds, including Ca de Bestiar (CB), Ca de Bou (CBOU), Podenco Ibicenco (PI), Ca Rater (CR), Ca Mè (CM), and Ca de Conills (CC), reveals insightful findings. CB showcases the highest mean number of alleles (6.17) and heterozygosity values, with significant deviations from Hardy-Weinberg equilibrium (HWE) observed in five markers, indicating high intra-racial genetic diversity (average observed heterozygosity (Ho) = 0.754, expected heterozygosity (He) = 0.761). In contrast, CBOU presents the lowest mean number of alleles (5.05) and heterozygosity values, coupled with moderate polymorphic information content (PIC) values and a moderate level of intra-racial genetic diversity (average Ho = 0.313, He = 0.394). PI demonstrates moderate genetic diversity with an average of 5.75 alleles and highly informative PIC values, while CR displays robust genetic diversity with an average of 6.61 alleles and deviations from equilibrium, indicating potential risks of inbreeding (average Ho = 0.563, He = 0.658). CM exhibits moderate genetic diversity and deviations from equilibrium, similar to CBOU, with an average of 6.5 alleles and moderate PIC values (average Ho = 0.598, He = 0.676). Conversely, CC shows a wider range of allelic diversity and deviations from equilibrium (average Ho = 0.611, He = 0.706), suggesting a more diverse genetic background. Inter-racial analysis underscores distinct genetic differentiation between breeds, emphasizing the importance of informed breeding decisions and proactive genetic management strategies to preserve diversity, promote breed health, and ensure long-term sustainability across all breeds studied.

Breed identification using breed-informative SNPs and machine learning based on whole genome sequence data and SNP chip data

BACKGROUND

Breed identification is useful in a variety of biological contexts. Breed identification usually involves two stages, i.e., detection of breed-informative SNPs and breed assignment. For both stages, there are several methods proposed. However, what is the optimal combination of these methods remain unclear. In this study, using the whole genome sequence data available for 13 cattle breeds from Run 8 of the 1,000 Bull Genomes Project, we compared the combinations of three methods (Delta, F_ST, and I_n) for breed-informative SNP detection and five machine learning methods (KNN, SVM, RF, NB, and ANN) for breed assignment with respect to different reference population sizes and difference numbers of most breed-informative SNPs. In addition, we evaluated the accuracy of breed identification using SNP chip data of different densities.

RESULTS

We found that all combinations performed quite well with identification accuracies over 95% in all scenarios. However, there was no combination which performed the best and robust across all scenarios. We proposed to integrate the three breed-informative detection methods, named DFI, and integrate the three machine learning methods, KNN, SVM, and RF, named KSR. We found that the combination of these two integrated methods outperformed the other combinations with accuracies over 99% in most cases and was very robust in all scenarios. The accuracies from using SNP chip data were only slightly lower than that from using sequence data in most cases.

CONCLUSIONS

The current study showed that the combination of DFI and KSR was the optimal strategy. Using sequence data resulted in higher accuracies than using chip data in most cases. However, the differences were generally small. In view of the cost of genotyping, using chip data is also a good option for breed identification.

RNA-Seq Provides Novel Genomic Resources for Noug (Guizotia abyssinica) and Reveals Microsatellite Frequency and Distribution in Its Transcriptome

Genomic resources and tools are essential for improving crops and conserving their genetic resources. Guizotia abyssinica (noug), an outcrossing edible oilseed crop, has highly limited genomic resources. Hence, RNA-Seq based transcriptome sequencing of 30 noug genotypes was performed to generate novel genomic resources and assess their usefulness. The genotypes include self-compatible and self-incompatible types, which differ in maturity time, photoperiod sensitivity, or oil content and quality. RNA-Seq was performed on Illumina HiSeq 2500 platform, and the transcript was reconstructed de novo, resulting in 409,309 unigenes. The unigenes were characterized for simple sequence repeats (SSRs), and served as a reference for single nucleotide polymorphism (SNP) calling. In total, 40,776 SSRs were identified in 35,639 of the 409,309 unigenes. Of these, mono, di, tri, tetra, penta and hexanucleotide repeats accounted for 55.4, 20.8, 21.1, 2.3, 0.2, and 0.2%, respectively. The average G+C content of the unigenes and their SSRs were 40 and 22.1%, respectively. The vast majority of mononucleotide repeat SSRs (97%) were of the A/T type. AG/CT and CCA/TGG were the most frequent di and trinucleotide repeat SSRs. A different number of single nucleotide polymorphism (SNP) loci were discovered in each genotype, of which 1,687 were common to all 30 genotypes and 5,531 to 28 of them. The mean observed heterozygosity of the 5,531 SNPs was 0.22; 19.4% of them had polymorphism information content above 0.30 while 17.2% deviated significantly from Hardy-Weinberg equilibrium (P < 0.05). In both cluster and principal coordinate analyses, the genotypes were grouped into four major clusters. In terms of population structure, the genotypes are best represented by three genetic populations, with significant admixture within each. Genetic similarity between self-compatible genotypes was higher, due to the narrow genetic basis, than that between self-incompatible genotypes. The genotypes that shared desirable characteristics, such as early maturity, and high oil content were found to be genetically diverse, and hence superior cultivars with multiple desirable traits can be developed through crossbreeding. The genomic resources developed in this study are vital for advancing research in noug, such as genetic linkage mapping and genome-wide association studies, which could lead to genomic-led breeding.

New Transcriptome-Based SNP Markers for Noug (Guizotia abyssinica) and Their Conversion to KASP Markers for Population Genetics Analyses

The development and use of genomic resources are essential for understanding the population genetics of crops for their efficient conservation and enhancement. Noug (Guizotia abyssinica) is an economically important oilseed crop in Ethiopia and India. The present study sought to develop new DNA markers for this crop. Transcriptome sequencing was conducted on two genotypes and 628 transcript sequences containing 959 single nucleotide polymorphisms (SNPs) were developed. A competitive allele-specific PCR (KASP) assay was developed for the SNPs and used for genotyping of 24 accessions. A total of 554 loci were successfully genotyped across the accessions, and 202 polymorphic loci were used for population genetics analyses. Polymorphism information content (PIC) of the loci varied from 0.01 to 0.37 with a mean of 0.24, and about 49% of the loci showed significant deviation from the Hardy-Weinberg equilibrium. The mean expected heterozygosity was 0.27 suggesting moderately high genetic variation within accessions. Low but significant differentiation existed among accessions (FST = 0.045, p < 0.0001). Landrace populations from isolated areas may have useful mutations and should be conserved and used in breeding this crop. The genomic resources developed in this study were shown to be useful for population genetics research and can also be used in, e.g., association genetics.

Mining whole genome sequence data to efficiently attribute individuals to source populations

Whole genome sequence (WGS) data could transform our ability to attribute individuals to source populations. However, methods that efficiently mine these data are yet to be developed. We present a minimal multilocus distance (MMD) method which rapidly deals with these large data sets as well as methods for optimally selecting loci. This was applied on WGS data to determine the source of human campylobacteriosis, the geographical origin of diverse biological species including humans and proteomic data to classify breast cancer tumours. The MMD method provides a highly accurate attribution which is computationally efficient for extended genotypes. These methods are generic, easy to implement for WGS and proteomic data and have wide application.

Testing theoretical metapopulation conditions with genotypic data from Boreal Chorus Frogs (Pseudacris maculata)

The metapopulation concept has far-reaching implications in ecology and conservation biology. Hanski’s criteria operationally define metapopulations, yet testing them is hindered by logistical and financial constraints inherent to the collection of long-term demographic data. Hence, ecologists and conservationists often assume metapopulation existence for dispersal-limited species that occupy patchy habitats. To advance understanding of metapopulation theory and improve conservation of metapopulations, we used population and landscape genetic tools to develop a methodological framework for evaluating Hanski’s criteria. We used genotypic data (11 microsatellite loci) from a purported metapopulation of Boreal Chorus Frogs (Pseudacris maculata (Agassiz, 1850)) in Colorado, U.S.A., to test Hanski’s four criteria. We found support for each criterion: (1) significant genetic differentiation between wetlands, suggesting distinct breeding populations; (2) wetlands had small effective population sizes and recent bottlenecks, suggesting populations do not experience long-term persistence; (3) population graphs provided evidence of gene flow between patches, indicating potential for recolonization; and (4) multiscale bottleneck analyses suggest asynchrony, indicating that simultaneous extinction of all populations was unlikely. Our methodological framework provides a logistically and financially feasible alternative to long-term demographic data for identifying amphibian metapopulations.

Historical baselines in marine bioinvasions: Implications for policy and management

The human-mediated introduction of marine non-indigenous species is a centuries- if not millennia-old phenomenon, but was only recently acknowledged as a potent driver of change in the sea. We provide a synopsis of key historical milestones for marine bioinvasions, including timelines of (a) discovery and understanding of the invasion process, focusing on transfer mechanisms and outcomes, (b) methodologies used for detection and monitoring, (c) approaches to ecological impacts research, and (d) management and policy responses. Early (until the mid-1900s) marine bioinvasions were given little attention, and in a number of cases actively and routinely facilitated. Beginning in the second half of the 20th century, several conspicuous non-indigenous species outbreaks with strong environmental, economic, and public health impacts raised widespread concerns and initiated shifts in public and scientific perceptions. These high-profile invasions led to policy documents and strategies to reduce the introduction and spread of non-indigenous species, although with significant time lags and limited success and focused on only a subset of transfer mechanisms. Integrated, multi-vector management within an ecosystem-based marine management context is urgently needed to address the complex interactions of natural and human pressures that drive invasions in marine ecosystems.

Rapid evolution in insect pests: the importance of space and time in population genomics studies

Pest species in agroecosystems often exhibit patterns of rapid evolution to environmental and human-imposed selection pressures. Although the role of adaptive processes is well accepted, few insect pests have been studied in detail and most research has focused on selection at insecticide resistance candidate genes. Emerging genomic datasets provide opportunities to detect and quantify selection in insect pest populations, and address long-standing questions about mechanisms underlying rapid evolutionary change. We examine the strengths of recent studies that stratify population samples both in space (along environmental gradients and comparing ancestral vs. derived populations) and in time (using chronological sampling, museum specimens and comparative phylogenomics), resulting in critical insights on evolutionary processes, and providing new directions for studying pests in agroecosystems.

Microsatellite marker development from next-generation sequencing in the New England cottontail (Sylvilagus transitionalis) and cross-amplification in the eastern cottontail (S. floridanus)

Objective

The New England cottontail (Sylvilagus transitionalis) is a species of high conservation priority in the Northeastern United States, and was a candidate for federal listing under the Endangered Species Act until a recent decision determined that conservation actions were sufficient to preclude listing. The aim of this study was to develop a suite of microsatellite loci to guide future research efforts such as the analysis of population genetic structure, genetic variation, dispersal, and genetic mark-recapture population estimation.

Results

Thirty-five microsatellite markers containing tri- and tetranucleotide sequences were developed from shotgun genomic sequencing of tissue from S. transitionalis, S. obscurus, and S. floridanus. These loci were screened in n = 33 wild S. transitionalis sampled from a population in eastern Massachusetts, USA. Thirty-two of the 35 loci were polymorphic with 2–6 alleles, and observed heterozygosities of 0.06–0.82. All loci conformed to Hardy–Weinberg Equilibrium proportions and there was no evidence of linkage disequilibrium or null alleles. Primers for 33 of the 35 loci amplified DNA extracted from n = 6 eastern cottontail (S. floridanus) samples, of which nine revealed putative species-diagnostic alleles. These loci will provide a useful tool for conservation genetics investigations of S. transitionalis and a potential diagnostic species assay for differentiating sympatric eastern and New England cottontails.

Microsatellite markers for Cryptolestes ferrugineus (Coleoptera: Laemophloeidae) and other Cryptolestes species

Cryptolestes ferrugineus (Stephens, 1831) is an important insect pest of stored products. Due to its broad host range, short life cycle, and high reproductive capacity, this species has rapidly colonized temperate and tropical regions around the world. In this study, we isolated 18 novel polymorphic microsatellite loci from an enriched genomic library based on a biotin/streptavidin capture protocol. These loci will be useful tool to better understand the genetic structure and migration patterns of C. ferrugineus throughout the world. The genetic parameters were estimated based on 80 individual C. ferrugineus from two natural populations. The results revealed that 18 loci were different polymorphic levels. The numbers of alleles ranged from 3 to 12, and eleven loci demonstrated polymorphic information contents greater than 0.5. The observed (H O) and expected (H E) heterozygosities ranged from 0.051 to 0.883 and 0.173 to 0.815, respectively. Five locus/population combinations significantly deviated from Hardy-Weinberg equilibrium. We also demonstrated the potential utility of the C. ferrugineus microsatellites as population and species markers for four additional Cryptolestes species.

Inferring ancestry from population genomic data and its applications

Ancestry inference is a frequently encountered problem and has many applications such as forensic analyses, genetic association studies, and personal genomics. The main goal of ancestry inference is to identify an individual’s population of origin based on our knowledge of natural populations. Because both self-reported ancestry in humans or the sampling location of an organism can be inaccurate for this purpose, the use of genetic markers can facilitate accurate and reliable inference of an individual’s ancestral origins. At a higher level, there are two different paradigms in ancestry inference: global ancestry inference which tries to compute the genome-wide average of the population contributions and local ancestry inference which tries to identify the regional ancestry of a genomic segment. In this mini review, I describe the numerous approaches that are currently available for both kinds of ancestry inference from population genomic datasets. I first describe the general ideas underlying such inference methods and their relationship to one another. Then, I describe practical applications in which inference of ancestry has proven useful. Lastly, I discuss challenges and directions for future research work in this area.

Genetic structure in the coral, Montastraea cavernosa: assessing genetic differentiation among and within Mesophotic reefs

Mesophotic coral reefs (30–150 m) have recently received increased attention as a potential source of larvae (e.g., the refugia hypothesis) to repopulate a select subset of the shallow water (<30 m) coral fauna. To test the refugia hypothesis we used highly polymorphic Amplified Fragment Length Polymorphism (AFLP) markers as a means to assess small-scale genetic heterogeneity between geographic locations and across depth clines in the Caribbean coral, Montastraea cavernosa. Zooxanthellae-free DNA extracts of coral samples (N = 105) were analyzed from four depths, shallow (3–10 m), medium (15–25 m), deep (30–50 m) and very deep (60–90 m) from Little Cayman Island (LCI), Lee Stocking Island (LSI), Bahamas and San Salvador (SS), Bahamas which range in distance from 170 to 1,600 km apart. Using AMOVA analysis there were significant differences in Φ_ST values in pair wise comparisons between LCI and LSI. Among depths at LCI, there was significant genetic differentiation between shallow and medium versus deep and very deep depths in contrast there were no significant differences in Φ_ST values among depths at LSI. The assignment program AFLPOP, however, correctly assigned 95.7% of the LCI and LSI samples to the depths from which they were collected, differentiating among populations as little as 10 to 20 m in depth from one another. Discriminant function analysis of the data showed significant differentiation among samples when categorized by collection site as well as collection depth. F_ST outlier analyses identified 2 loci under positive selection and 3 under balancing selection at LCI. At LSI 2 loci were identified, both showing balancing selection. This data shows that adult populations of M. cavernosa separated by depths of tens of meters exhibits significant genetic structure, indicative of low population connectivity among and within sites and are not supplying successful recruits to adjacent coral reefs less than 30 m in depth.

Patterns of mitochondrial haplotype diversity in the invasive pest Epiphyas postvittana (Lepidoptera: Tortricidae)

The light brown apple moth, Epiphyas postvittana (Walker) (Lepidoptera: Tortricidae), is a horticultural pest of Australia and New Zealand that has more recently invaded Hawaii, Europe, and California. A 2,216-bp region of the mitochondrial genome containing the cytochrome oxidase I and II genes was sequenced from 752 individuals. Haplotype network analyses revealed a major split between a predominantly Western Australian clade and all other samples, suggestive of either a deep genetic divergence or a cryptic species. Nucleotide and haplotype diversity were highest in the country of origin, Australia, and in New Zealand populations, with evidence of haplotype sharing between New Zealand and Tasmania. Nucleotide and haplotype diversity were higher in California than within the British Isles or Hawaii. From the total of 96 haplotypes, seven were found in California, of which four were private. Within California, there have been at least two introductions; based on genetic diversity we were unable to assign a likely source for a single moth found and eradicated in Los Angeles in 2007; however, our data suggest it is unlikely that Hawaii and the British Isles are sources of the major E. postvittana population found throughout the rest of the state since 2006.

Evaluation of approaches for identifying population informative markers from high density SNP chips

Background

Genetic markers can be used to identify and verify the origin of individuals. Motivation for the inference of ancestry ranges from conservation genetics to forensic analysis. High density assays featuring Single Nucleotide Polymorphism (SNP) markers can be exploited to create a reduced panel containing the most informative markers for these purposes. The objectives of this study were to evaluate methods of marker selection and determine the minimum number of markers from the BovineSNP50 BeadChip required to verify the origin of individuals in European cattle breeds. Delta, Wright's F_ST, Weir & Cockerham's F_ST and PCA methods for population differentiation were compared. The level of informativeness of each SNP was estimated from the breed specific allele frequencies. Individual assignment analysis was performed using the ranked informative markers. Stringency levels were applied by log-likelihood ratio to assess the confidence of the assignment test.

Results

A 95% assignment success rate for the 384 individually genotyped animals was achieved with < 80, < 100, < 140 and < 200 SNP markers (with increasing stringency threshold levels) across all the examined methods for marker selection. No further gain in power of assignment was achieved by sampling in excess of 200 SNP markers. The marker selection method that required the lowest number of SNP markers to verify the animal's breed origin was Wright's F_ST (60 to 140 SNPs depending on the chosen degree of confidence). Certain breeds required fewer markers (< 100) to achieve 100% assignment success. In contrast, closely related breeds require more markers (~200) to achieve > 95% assignment success. The power of assignment success, and therefore the number of SNP markers required, is dependent on the levels of genetic heterogeneity and pool of samples considered.

Conclusions

While all SNP selection methods produced marker panels capable of breed identification, the power of assignment varied markedly among analysis methods. Thus, with effective exploration of available high density genetic markers, a diagnostic panel of highly informative markers can be produced.

Shifting dispersal modes at an expanding species' range margin

While it is generally recognized that noncontiguous (long-distance) dispersal of small numbers of individuals is important for range expansion over large geographic areas, it is often assumed that colonization on more local scales proceeds by population expansion and diffusion dispersal (larger numbers of individuals colonizing adjacent sites). There are few empirical studies of dispersal modes at the front of expanding ranges, and very little information is available on dispersal dynamics at smaller geographic scales where we expect contiguous (diffusion) dispersal to be prevalent. We used highly polymorphic genetic markers to characterize dispersal modes at a local geographic scale for populations at the edge of the range of a newly invasive grass species (Brachypodium sylvaticum) that is undergoing rapid range expansion in the Pacific Northwest of North America. Comparisons of Bayesian clustering of populations, patterns of genetic diversity, and gametic disequilibrium indicate that new populations are colonized ahead of the invasion front by noncontiguous dispersal from source populations, with admixture occurring as populations age. This pattern of noncontiguous colonization was maintained even at a local scale. Absence of evidence for dispersal among adjacent pioneer sites at the edge of the expanding range of this species suggests that pioneer populations undergo an establishment phase during which they do not contribute emigrants for colonization of neighbouring sites. Our data indicate that dispersal modes change as the invasion matures: initial colonization processes appear to be dominated by noncontiguous dispersal from only a few sources, while contiguous dispersal may play a greater role once populations become established.

Detecting populations in the 'ambiguous' zone: kinship-based estimation of population structure at low genetic divergence

Identifying population structure is one of the most common and important objectives of spatial analyses using population genetic data. Population structure is detected either by rejecting the null hypothesis of a homogenous distribution of genetic variation, or by estimating low migration rates. Issues arise with most current population genetic inference methods when the genetic divergence is low among putative populations. Low levels of genetic divergence may be as a result of either high ongoing migration or historic high migration but no current, ongoing migration. We direct attention to recent developments in the use of the tempo-spatial distribution of closely related individuals to detect population structure or estimate current migration rates. These 'kinship-based' approaches complement more traditional population-based genetic inference methods by providing a means to detect population structure and estimate current migration rates when genetic divergence is low. However, for kinship-based methods to become widely adopted, formal estimation procedures applicable to a range of species life histories are needed.

Single nucleotide polymorphisms for assessing genetic diversity in castor bean (Ricinus communis)

BACKGROUND

Castor bean (Ricinus communis) is an agricultural crop and garden ornamental that is widely cultivated and has been introduced worldwide. Understanding population structure and the distribution of castor bean cultivars has been challenging because of limited genetic variability. We analyzed the population genetics of R. communis in a worldwide collection of plants from germplasm and from naturalized populations in Florida, U.S. To assess genetic diversity we conducted survey sequencing of the genomes of seven diverse cultivars and compared the data to a reference genome assembly of a widespread cultivar (Hale). We determined the population genetic structure of 676 samples using single nucleotide polymorphisms (SNPs) at 48 loci.

RESULTS

Bayesian clustering indicated five main groups worldwide and a repeated pattern of mixed genotypes in most countries. High levels of population differentiation occurred between most populations but this structure was not geographically based. Most molecular variance occurred within populations (74%) followed by 22% among populations, and 4% among continents. Samples from naturalized populations in Florida indicated significant population structuring consistent with local demes. There was significant population differentiation for 56 of 78 comparisons in Florida (pairwise population phiPT values, p < 0.01).

CONCLUSION

Low levels of genetic diversity and mixing of genotypes have led to minimal geographic structuring of castor bean populations worldwide. Relatively few lineages occur and these are widely distributed. Our approach of determining population genetic structure using SNPs from genome-wide comparisons constitutes a framework for high-throughput analyses of genetic diversity in plants, particularly in species with limited genetic diversity.

Genetic perspectives on marine biological invasions

The extent to which the geographic distributions of marine organisms have been reshaped by human activities remains underappreciated, and so does, consequently, the impact of invasive species on marine ecosystems. The application of molecular genetic data in fields such as population genetics, phylogeography, and evolutionary biology have improved our ability to make inferences regarding invasion histories. Genetic methods have helped to resolve longstanding questions regarding the cryptogenic status of marine species, facilitated recognition of cryptic marine biodiversity, and provided means to determine the sources of introduced marine populations and to begin to recover the patterns of anthropogenic reshuffling of the ocean's biota. These approaches stand to aid materially in the development of effective management strategies and sustainable science-based policies. Continued advancements in the statistical analysis of genetic data promise to overcome some existing limitations of current approaches. Still other limitations will be best addressed by concerted collaborative and multidisciplinary efforts that recognize the important synergy between understanding the extent of biological invasions and coming to a more complete picture of both modern-day and historical marine biogeography.

Genetic variation in the invasive avian parasite, Philornis downsi (Diptera, Muscidae) on the Galápagos archipelago

BACKGROUND

Understanding the dispersal and genetic structure of invasive insects across islands is important for designing management plans that are appropriate at spatial and temporal scales. For invasive parasites, population dynamics are largely determined by the distribution and density of their host species. The introduced parasitic fly, Philornis downsi, parasitises nestlings of endemic birds on all major islands of the Galápagos archipelago. The fly's high mortality and fitness impacts are of conservation concern for vulnerable and declining species of Darwin's finches. Using microsatellite data in Bayesian clustering and landscape genetic analyses, we examine gene flow and dispersal in P. downsi between three islands and across habitats (highlands, lowlands) and examine for the presence of population bottlenecks. We also examine variation at the mitochondrial gene CO1 across islands to establish if cryptic species were present.

RESULTS

Both the mitochondrial and microsatellite data were consistent with there being a single species across islands. We found low genetic differentiation between islands and strong evidence for inter-island gene flow, or shared recent ancestry among individuals. Landscape genetic analysis identified two genetic clusters: one encompassing Santa Cruz and Isabela, and one on Floreana Island. There was no evidence of genetic differentiation between habitats and molecular variance was mainly attributable to within individuals. The combined P. downsi population was found to have undergone a population bottleneck.

CONCLUSION

Philornis downsi populations have high connectivity within and between islands, with low levels of genetic differentiation between Floreana and the other two islands examined. The genetic bottleneck found across islands suggests there was a small founding population or few introduction events of P. downsi. The high dispersal capacity and wide habitat use of P. downsi highlights the significant threat that this parasite poses to the Galápagos avifauna. Our findings are relevant for assessing the viability of methods to control P. downsi on Galápagos, such as the sterile insect technique.

Genetic evidence that culling increases badger movement: implications for the spread of bovine tuberculosis

The Eurasian badger (Meles meles) has been implicated in the transmission of bovine tuberculosis (TB, caused by Mycobacterium bovis) to cattle. However, evidence suggests that attempts to reduce the spread of TB among cattle in Britain by culling badgers have mixed effects. A large-scale field experiment (the randomized badger culling trial, RBCT) showed that widespread proactive badger culling reduced the incidence of TB in cattle within culled areas but that TB incidence increased in adjoining areas. Additionally, localized reactive badger culling increased the incidence of TB in cattle. It has been suggested that culling-induced perturbation of badger social structure may increase individual movements and elevate the risk of disease transmission between badgers and cattle. Field studies support this hypothesis, by demonstrating increases in badger group ranges and the prevalence of TB infection in badgers following culling. However, more evidence on the effect of culling on badger movements is needed in order to predict the epidemiological consequences of this control strategy. Here, analysis of the genetic signatures of badger populations in the RBCT revealed increased dispersal following culling. While standard tests provided evidence for greater dispersal after culling, a novel method indicated that this was due to medium- and long-distance dispersal, in addition to previously reported increases in home-range size. Our results also indicated that, on average, badgers infected with M. bovis moved significantly farther than did uninfected badgers. A disease control strategy that included culling would need to take account of the potentially negative epidemiological consequences of increased badger dispersal.

Introduction history of Drosophila subobscura in the New World: a microsatellite-based survey using ABC methods

Drosophila subobscura is a Palearctic species that was first observed in South and North America in the early 1980s, and that rapidly invaded broad latitudinal ranges on both continents. To trace the source and history of this invasion, we obtained genotypic data on nine microsatellite loci from two South American, two North American and five European populations of D. subobscura. We analysed these data with traditional statistics as well as with an approximate Bayesian computation (ABC) framework. ABC methods yielded the strongest support for the scenario involving a serial introduction with founder events from Europe into South America, and then from South America into North America. Stable effective population size of the source population was very large (around one million individuals), and the propagule size was notably smaller for the introduction into South America (i.e. high bottleneck severity index with only a few effective founders) but considerably larger for the subsequent introduction into North America (i.e. low bottleneck severity index with around 100-150 effective founders). Finally, the Mediterranean region of Europe (and most likely Barcelona from the localities so far analysed) is proposed as the source of the New World flies, based on mean individual assignment statistics.

High genetic diversity and population differentiation in Boechera fecunda, a rare relative of Arabidopsis

Conservation of endangered species becomes a critical issue with the increasing rates of extinction. In this study, we use 13 microsatellite loci and 27 single-copy nuclear loci to investigate the population genetics of Boechera fecunda, a rare relative of Arabidopsis thaliana, known from only 21 populations in Montana. We investigated levels of genetic diversity and population structure in comparison to its widespread congener, Boechera stricta, which shares similar life history and mating system. Despite its rarity, B. fecunda had levels of genetic diversity similar to B. stricta for both microsatellites and nucleotide polymorphism. Populations of B. fecunda are highly differentiated, with a majority of genetic diversity existing among populations (F(ST) = 0.57). Differences in molecular diversity and allele frequencies between western and eastern population groups suggest they experienced very different evolutionary histories.

GENETIC DIVERSITY AND EXPANDING NONINDIGENOUS RANGE OF THE RHIZOCEPHALAN LOXOTHYLACUS PANOPAEI PARASITIZING MUD CRABS IN THE WESTERN NORTH ATLANTIC

Nonindigenous parasite introductions and range expansions have become a major concern because of their potential to restructure communities and impact fisheries. Molecular markers provide an important tool for reconstructing the pattern of introduction. The parasitic castrator Loxothylacus panopaei, a rhizocephalan barnacle, infects estuarine mud crabs in the Gulf of Mexico and southeastern Florida. A similar parasite introduced into Chesapeake Bay before 1964, presumably via infected crabs associated with oysters from the Gulf of Mexico, was identified as L. panopaei. Our samples of this species during 2004 and 2005 show that the introduced range has expanded as far south as Edgewater, Florida, just north of the northern endemic range limit. The nonindigenous range expanded southward at a rate of up to 165 km/yr with relatively high prevalence, ranging from 30 to 93%. Mitochondrial DNA sequences from the cytochrome oxidase I gene showed that these nonindigenous L. panopaei are genetically distinct from the endemic parasites in southeastern Florida and the eastern Gulf of Mexico. The genetic difference was also associated with distinct host spectra. These results are incompatible with an eastern Gulf source population, but suggest that unrecognized genetic and phenotypic population structure may occur among Gulf of Mexico populations of Loxothvlacus.

Historical range expansion determines the phylogenetic diversity introduced during contemporary species invasion

For a species rapidly expanding its geographic range, such as during biological invasion, most alleles in the introduced range will have their evolutionary origins in the native range. Yet, the way in which historical processes occurring over evolutionary time in the native range contribute to the diversity sampled during contemporary invasion is largely unknown. We used chloroplast DNA (cpDNA) gene genealogies and coalescent methods to study two congeneric plants, Silene latifolia and S. vulgaris. We examined how phylogenetic diversity was shaped by demographic growth and historical range expansions in the native European range, and how this history affected the diversity sampled during their recent invasion of North America. Genealogies from both species depart from neutrality, likely as a result of demographic expansion in the ancestral range, the timing of which corresponds to shortly after each species originated. However, the species differ in the spatial distribution of cpDNA lineages across the native range. Silene latifolia shows a highly significant phylogeographic structure that most likely reflects different avenues of the post-glacial expansion into northern Europe from Mediterranean refugia. By contrast, cpDNA lineages in S. vulgaris have been widely scattered across Europe during, or since, the most recent post-glacial expansion. These different evolutionary histories resulted in dramatic differences in how phylogenetic diversity was sampled during invasion of North America. In S. latifolia, relatively few, discrete invasion events from a structured native range resulted in a rather severe genetic bottleneck, but also opportunities for admixture among previously isolated lineages. In S. vulgaris, lack of genetic structure was accompanied by more representative sampling of phylogenetic diversity during invasion, and reduced potential for admixture. Our results provide clear insights into how historical processes may feed forward to influence the phylogenetic diversity of species invading new geographic ranges.

Detecting female precise natal philopatry in green turtles using assignment methods

It is well established that sea turtles return to natal rookeries to mate and lay their eggs, and that individual females are faithful to particular nesting sites within the rookery. Less certain is whether females are precisely returning to their natal beach. Attempts to demonstrate such precise natal philopatry with genetic data have had mixed success. Here we focused on the green turtles of three nesting sites in the Ascension Island rookery, separated by 5-15 km. Our approach differed from previous work in two key areas. First, we used male microsatellite data (five loci) reconstructed from samples collected from their offspring (N = 17) in addition to data for samples taken directly from females (N = 139). Second, we employed assignment methods in addition to the more traditional F-statistics. No significant genetic structure could be demonstrated with F(ST). However, when average assignment probabilities of females were examined, those for nesting populations in which they were sampled were indeed significantly higher than their probabilities for other populations (Mann-Whitney U-test: P < 0.001). Further evidence was provided by a significant result for the mAI(C) test (P < 0.001), supporting greater natal philopatry for females compared with males. The results suggest that female natal site fidelity was not sufficient for significant genetic differentiation among the nesting populations within the rookery, but detectable with assignment tests.

Genetic Structure and Differentiation of Three Chinese Indigenous Cattle Populations

Levels of genetic differentiation, gene flow, and genetic structure of three indigenous cattle populations (Luxi, Bohai, and Minnan) and two reference cattle populations (Chinese Holstein and Qinhai yak) in China were estimated using the information from 12 microsatellites, and 141 microsatellite alleles were identified. The mean number of alleles per locus ranged from 2.9005 in yak to 4.9722 in Holstein. The observed heterozygosity ranged from 0.5325 (yak) to 0.7719 (Holstein); 29 private alleles were detected. The global heterozygote deficit across all populations amounted to 58.5% (p < 0.001). The overall significant (p < 0.001) deficit of heterozygotes because of inbreeding within breeds amounted to 43.2%. The five cattle populations were highly differentiated (F (st) = 26.9%, p < 0.001) at all loci. The heterozygote deficit within the population was highest in Luxi cattle and lowest in yak. The average number of effective migrants exchanged per generation was highest (1.149) between Luxi and Holstein, and lowest (0.509) between Luxi and yak. With the application of prior population information, cluster analysis achieved posterior probabilities from 91% to 98% of correctly assigning individuals to populations. Combining the information of cluster analysis, gene flow, and Structure analysis, the five cattle populations belong to three genetic clusters, a taurine (Luxi and Chinese Holstein), a zebu (Bohai and Minnan), and a yak cluster. This indicates that Bohai black is closer to Bos indicus than Luxi cattle. The evolution and development of three indigenous cattle populations are discussed.