Genotyping by sequencing suggests overwintering of Peronospora destructor in southwestern Québec,Canada

Several Peronospora species are carried by wind over short and long distances, from warmer climates where they survive on living plants to cooler climates. In eastern Canada, this annual flow of sporangia was thought to be the main source of Peronospora destructor responsible for onion downy mildew. However, the results of a recent study showed that the increasing frequency of onion downy mildew epidemics in eastern Canada is associated with warmer autumns, milder winters, and previous year disease severity, suggesting overwintering of the inoculum in an area where the pathogen is not known to be endogenous. In this study, genotyping by sequencing was used to investigate the population structure of P. destructor at the landscape scale. The study focused on a particular region of southwestern Québec—Les Jardins de Napierville—to determine if the populations were clonal and regionally differentiated. The data were characterized by a high level of linkage disequilibrium, characteristic of clonal organisms. Consequently, the null hypothesis of random mating was rejected when tested on predefined or nonpredefined populations, indicating that linkage disequilibrium was not a function of population structure and suggesting a mixed reproduction mode. Discriminant analysis of principal components performed with predefined population assignment allowed grouping P. destructor isolates by geographical regions, while analysis of molecular variance confirmed that this genetic differentiation was significant at the regional level. Without using a priori population assignment, isolates were clustered into four genetic clusters. These results represent a baseline estimate of the genetic diversity and population structure of P. destructor.     

Effect of Mating Structure on Variation in Linkage Disequilibrium

Measurement of linkage disequilibrium involves two sampling processes. First, there is the sampling of gametes in the population to form successive generations, and this generates disequilibrium dependent on the effective population size (N_e) and the mating structure. Second, there is sampling of a finite number (n) of individuals to estimate the population disequilibrium.——Two-locus descent measures are used to describe the mating system and are transformed to disequilibrium moments at the final sampling. Approximate eigenvectors for the transition matrix of descent measures are used to obtain formulae for the variance of the observed disequilibria as a function of N_e, mating structure, n, and linkage or recombination parameter.——The variance of disequilibrium is the same for monoecious populations with or without random selfing and for dioecious populations with random pairing for each progeny. With monogamy, the variance is slightly higher, the proportional difference being greater for unlinked loci.     

Construction of a sequence-based bin map and mapping of QTLs for downy mildew resistance at four developmental stages in Chinese cabbage (<Emphasis Type="Italic">Brassica rapa</Emphasis> L. ssp. <Emphasis Type="Italic">pekinensis</Emphasis>)

Specific-locus amplified fragment sequencing is a high-resolution method for genetic mapping, genotyping, and single nucleotide polymorphism (SNP) marker discovery. Previously, a major QTL for downy mildew resistance, BraDM, was mapped to linkage group A08 in a doubled-haploid population derived from Chinese cabbage lines 91–112 and T12–19. The aim of the present study was to improve the linkage map and identify the genetic factors involved in downy mildew resistance. We detected 53,692 high quality SLAFs, of which 7230 were polymorphic, and 3482 of the polymorphic markers were used in genetic map construction. The final map included 1064 bins on ten linkage groups and was 858.98 cM in length, with an average inter-locus distance of 0.81 cM. We identified six QTLs that are involved in downy mildew resistance. The four major QTLs, sBrDM8, yBrDM8, rBrDM8, and hBrDM8, for resistance at the seedling, young plant, rosette, and heading stages were mapped to A08, and are identical to BraDM. The two minor resistance QTLs, rBrDM6 (A06) and hBrDM4 (A04), were active at the rosette and heading stages. The major QTL sBrDM8 defined a physical interval of ~228 Kb on A08, and a serine/threonine kinase family gene, Bra016457, was identified as the possible candidate gene. We report here the first high-density bin map for Chinese cabbage, which will facilitate mapping QTLs for economically important traits and SNP marker development. Our results also expand knowledge of downy mildew resistance in Chinese cabbage and provide three SNP markers (A08-709, A08-028, and A08-018) that we showed to be effective when used in MAS to breed for downy mildew resistance in B. rapa.     

The Effects of Overdominance on Linkage in a Multilocus System

Computer simulations were performed with overdominant multiple alleles among tightly linked multiple loci under a multiplicative fitness model. The quantity X²/N(n — 1) was introduced as a new measure of linkage disequilibrium which, unlike previously available measures, can be applied to multiple allele models, where N is the sample size, and n is the number of alleles at the locus possessing fewest alleles. Simulations showed that (1) With multiple (three or four) alleles, the approach to stable disequilibrium is slower and the amount of disequilibrium established is weaker than in a two allele system. (2) The number of complementary chromosomes is a function of number of alleles and of population size. (3) As population size increases, the rate of the approach to stable disequilibrium is slower. (4) There is an optimum selection coefficient which minimizes the transient fixation probability of alleles when linkage is present. (5) The absence of linkage disequilibrium is in most cases not a practical method of testing the hypothesis of balancing selection of genetic polymorphisms because it depends strongly on population size in determining linkage disequilibria.     

A Genetic Map of Lettuce (Lactuca sativa L.) with Restriction Fragment Length Polymorphism, Isozyme, Disease Resistance and Morphological Markers

A detailed linkage map of lettuce was constructed using 53 genetic markers including 41 restriction fragment length polymorphism (RFLP) loci, five downy mildew resistance genes, four isozyme loci and three morphological markers. The genetic markers were distributed into nine linkage groups and cover 404 cM which may be 25-30% of the lettuce genome. The majority (31 of 34) of the RFLP probes detected single segregating loci, although seven of these may have been homologous to further monomorphic loci. When several loci were detected by a single probe, the loci were generally linked, suggesting tandem duplications. One probe, however, detected loci in three linkage groups suggesting translocations. The five downy mildew resistance genes (Dm1, Dm3, Dm4, Dm5/8 and Dm13), segregating in the Calmar x Kordaat cross, represented each of the four resistance gene linkage groups. Dm5/8 is flanked by two cDNA loci, each located 10 cM away. These flanking markers will be used to study the source of variation in downy mildew genes and are also part our strategy to clone resistance genes.     

Genetic mapping and localization of a major QTL for seedling resistance to downy mildew in Chinese cabbage (Brassica rapa ssp. pekinensis)

Downy mildew caused by the fungus Peronospora parisitica is a serious threat to members of the Brassicaceae family. Annually, a substantial loss of yield is caused by the widespread presence of this disease in warm and humid climates. The aim of this study was to localize the genetic factors affecting downy mildew resistance in Chinese cabbage (Brassica rapa ssp. pekinensis). To achieve this goal, we improved a preexisting genetic map of a doubled-haploid population derived from a cross between two diverse Chinese cabbage lines, 91-112 and T12-19, via microspore culture. Microsatellite simple sequence repeat (SSR) markers, isozyme markers, sequence-related amplified polymorphism markers, sequence-characterized amplified region markers and sequence-tagged-site markers were integrated into the previously published map to construct a composite Chinese cabbage map. In this way, the identities of linkage groups corresponding to the Brassica A genome reference map were established. The new map contains 519 markers and covers a total length of 1,070 cM, with an average distance between markers of 2.06 cM. All markers were designated as A1–A10 through alignment and orientation using 55 markers anchored to previously published B. rapa or B. napus reference maps. Of the 89 SSR markers mapped, 15 were newly developed from express sequence tags in Genbank. The phenotypic assay indicated that a single major gene controls seedling resistance to downy mildew, and that a major QTL was detected on linkage group A8 by both interval and MQM mapping methods. The RAPD marker K14-1030 and isozyme marker PGM flanked this major QTL in a region spanning 2.9 cM, and the SSR marker Ol12G04 was linked to this QTL by a distance of 4.36 cM. This study identified a potential chromosomal segment and tightly linked markers for use in marker-assisted selection to improve downy mildew resistance in Chinese cabbage.     

Construction of a reference linkage map of Vitis amurensis and genetic mapping of Rpv8, a locus conferring resistance to grapevine downy mildew

Downy mildew, caused by the oomycete Plasmopara viticola, is one of the major threats to grapevine. All traditional cultivars of grapevine (Vitis vinifera) are susceptible to downy mildew, the control of which requires regular application of fungicides. In contrast, many sources of resistance to P. viticola have been described in the Vitis wild species, among which is V. amurensis Rupr. (Vitaceae), a species originating from East Asia. A genetic linkage map of V. amurensis, based on 122 simple sequence repeat and 6 resistance gene analogue markers, was established using S1 progeny. This map covers 975?cM on 19 linkage groups, which represent 82% of the physical coverage of the V. vinifera reference genetic map. To measure the general level of resistance, the sporulation of P. viticola and the necrosis produced in response to infection, five quantitative and semi-quantitative parameters were scored 6?days post-inoculation on the S1 progeny. A quantitative trait locus (QTL) analysis allowed us to identify on linkage group 14 a major QTL controlling the resistance to downy mildew found in V. amurensis, which explained up to 86.3% of the total phenotypic variance. This QTL was named ??Resistance to Plasmopara viticola 8?? (Rpv8).     

Genetic Diversity of Plasmodium falciparum in Haiti: Insights from Microsatellite Markers

Hispaniola, comprising Haiti and the Dominican Republic, has been identified as a candidate for malaria elimination. However, incomplete surveillance data in Haiti hamper efforts to assess the impact of ongoing malaria control interventions. Characteristics of the genetic diversity of Plasmodium falciparum populations can be used to assess parasite transmission, which is information vital to evaluating malaria elimination efforts. Here we characterize the genetic diversity of P. falciparum samples collected from patients at seven sites in Haiti using 12 microsatellite markers previously employed in population genetic analyses of global P. falciparum populations. We measured multiplicity of infections, level of genetic diversity, degree of population geographic substructure, and linkage disequilibrium (defined as non-random association of alleles from different loci). For low transmission populations like Haiti, we expect to see few multiple infections, low levels of genetic diversity, high degree of population structure, and high linkage disequilibrium. In Haiti, we found low levels of multiple infections (12.9%), moderate to high levels of genetic diversity (mean number of alleles per locus = 4.9, heterozygosity = 0.61), low levels of population structure (highest pairwise F_st = 0.09 and no clustering in principal components analysis), and moderate linkage disequilibrium (ISA = 0.05, P<0.0001). In addition, population bottleneck analysis revealed no evidence for a reduction in the P. falciparum population size in Haiti. We conclude that the high level of genetic diversity and lack of evidence for a population bottleneck may suggest that Haiti’s P. falciparum population has been stable and discuss the implications of our results for understanding the impact of malaria control interventions. We also discuss the relevance of parasite population history and other host and vector factors when assessing transmission intensity from genetic diversity data.     

Seed Transmission of Pseudoperonospora cubensis

Pseudoperonospora cubensis, an obligate biotrophic oomycete causing devastating foliar disease in species of the Cucurbitaceae family, was never reported in seeds or transmitted by seeds. We now show that P. cubensis occurs in fruits and seeds of downy mildew-infected plants but not in fruits or seeds of healthy plants. About 6.7% of the fruits collected during 2012–2014 have developed downy mildew when homogenized and inoculated onto detached leaves and 0.9% of the seeds collected developed downy mildew when grown to the seedling stage. This is the first report showing that P. cubensis has become seed-transmitted in cucurbits. Species-specific PCR assays showed that P. cubensis occurs in ovaries, fruit seed cavity and seed embryos of cucurbits. We propose that international trade of fruits or seeds of cucurbits might be associated with the recent global change in the population structure of P. cubensis.     

Development of SSR markers linked to QTL reducing leaf hair density and grapevine downy mildew resistance in <Emphasis Type="Italic">Vitis vinifera</Emphasis>

Dense leaf hairs of grapevines have been known to act as a preexisting defense structure for preventing the incidence of grapevine downy mildew, because the pathogen, Plasmopara viticola, needs water to invade grapevines, and water is repelled by a hydrophobic surface due to dense leaf hairs. In the present study, we performed regression analyses of downy mildew resistance with leaf hair density using hybrids of Vitis labrusca origin and confirmed the effect of leaf hairs. Reducing the repelling effect of leaf hairs by detergent application canceled the effect of leaf hair, which confirmed the hypothesis. Thereafter, based on QTL analyses in the population of V. vinifera ‘Muscat of Alexandria’ × the interspecific hybrid ‘Campbell Early,’ we identified a major locus in linkage group (LG) 5 of ‘Muscat of Alexandria’ controlling leaf hair density. This locus was previously reported as a small effect QTL for downy mildew resistance, however we found that the locus had high LOD scores explaining 71.9%–78.5% of the phenotypic variance of leaf hairs. Moreover, this locus was detected as a QTL for downy mildew resistance. The effect of this locus was confirmed in two other hybrid populations. Finally, we could successfully identify three traditional V. vinifera table grapes ‘Muscat of Alexandria,’ ‘Katta Kurgan,’ and ‘Parkent’ as the origin of the allele linked to hairlessness by defining the SSR haplotypes. The use of this locus for marker-assisted selections would be a promising strategy for producing grapevines with resistance by preexisting defense structure.     

Cryptosporidium,Giardia, Cryptococcus,Pneumocystis Genetic Variability: Cryptic Biological Species or Clonal Near-Clades?

An abundant literature dealing with the population genetics and taxonomy of Giardia duodenalis, Cryptosporidium spp., Pneumocystis spp., and Cryptococcus spp., pathogens of high medical and veterinary relevance, has been produced in recent years. We have analyzed these data in the light of new population genetic concepts dealing with predominant clonal evolution (PCE) recently proposed by us. In spite of the considerable phylogenetic diversity that exists among these pathogens, we have found striking similarities among them. The two main PCE features described by us, namely highly significant linkage disequilibrium and near-clading (stable phylogenetic clustering clouded by occasional recombination), are clearly observed in Cryptococcus and Giardia, and more limited indication of them is also present in Cryptosporidium and Pneumocystis. Moreover, in several cases, these features still obtain when the near-clades that subdivide the species are analyzed separately (“Russian doll pattern”). Lastly, several sets of data undermine the notion that certain microbes form clonal lineages simply owing to a lack of opportunity to outcross due to low transmission rates leading to lack of multiclonal infections (“starving sex hypothesis”). We propose that the divergent taxonomic and population genetic inferences advanced by various authors about these pathogens may not correspond to true evolutionary differences and could be, rather, the reflection of idiosyncratic practices among compartmentalized scientific communities. The PCE model provides an opportunity to revise the taxonomy and applied research dealing with these pathogens and others, such as viruses, bacteria, parasitic protozoa, and fungi.     

Identity of the downy mildew pathogens of basil, coleus, and sage with implications for quarantine measures

The downy mildew pathogen of basil (Ocimum spp.) has caused considerable damage throughout the past five years, and an end to the epidemics is not in sight. The downy mildew of coleus (Solenostemon spp.) is just emerging and here we report that it was very recently introduced into Germany. Although it has been recognised that these pathogens are a major threat, the identity of the pathogens is still unresolved, and so it is difficult to devise quarantine measures against them. Using morphological comparison and molecular phylogenetic reconstructions we confirmed in this study that the downy mildews of basil and coleus are unrelated to Peronospora lamii, which is a common pathogen of the weed Lamium purpureum. In addition, we conclude by the investigation of the type specimen of P. swingleii and downy mildew specimens on Salvia officinalis that the newly occurring pathogens are not identical to P. swingleii on Salvia reflexa. The taxonomy of the downy mildew pathogens of hosts from the Lamiaceae and, in particular, from the tribes Mentheae and Elsholtzieae, is discussed, and a new species is described to accommodate the downy mildew pathogen of basil and coleus, which is the first downy mildew pathogen known to be parasitic to hosts of the tribe Ocimeae.     

Mapping and validation of powdery mildew resistance loci from spring wheat cv. Naxos with SNP markers

Powdery mildew, caused by Blumeria graminis f.sp. tritici, is a major wheat disease in maritime and temperate climates. Breeding for race-non-specific or partial resistance is a cost-effective and environmentally friendly disease control strategy. The German spring wheat cultivar Naxos has proven to be a good source for partial resistance to powdery mildew. The objectives of the present study were to map the resistance loci in Naxos with use of high-density SNP markers in the Shanghai3/Catbird x Naxos inbred line population and validate the results in a different genetic background; Soru#1 x Naxos. Both populations were genotyped with the Illumina iSelect 90K wheat chip, and integrated linkage maps developed by inclusion of previously genotyped SSR and DArT markers. With the new linkage maps, we detected a total of 12 QTL for powdery mildew resistance in Shanghai3/Catbird x Naxos, of which eight were derived from Naxos. Previously reported QTL on chromosome arms 1AS and 2BL were more precisely mapped and the SNP markers enabled discovery of new QTL on 1AL, 2AL, 5AS and 5AL. In the Soru#1 x Naxos population, four QTL for powdery mildew resistance were detected, of which three had resistance from Naxos. This mapping verified the 1AS and 2AL QTL detected in Shanghai3/Catbird x Naxos, and identified a new QTL from Naxos on 2BL. In conclusion, the improved linkage maps with SNP markers enabled discovery of new resistance QTL and more precise mapping of previously known QTL. Moreover, the results were validated in an independent genetic background.     

Daytime Solar Heating Controls Downy Mildew Peronospora belbahrii in Sweet Basil

The biotrophic oomycete Peronospora belbahrii causes a devastating downy mildew disease in sweet basil. Due to the lack of resistant cultivars current control measures rely heavily on fungicides. However, resistance to fungicides and strict regulation on their deployment greatly restrict their use. Here we report on a ‘green’ method to control this disease. Growth chamber studies showed that P. belbahrii could hardly withstand exposure to high temperatures; exposure of spores, infected leaves, or infected plants to 35-45°C for 6-9 hours suppressed its survival. Therefore, daytime solar heating was employed in the field to control the downy mildew disease it causes in basil. Covering growth houses of sweet basil already infected with downy mildew with transparent infra-red-impermeable, transparent polyethylene sheets raised the daily maximal temperature during sunny hours by 11-22°C reaching 40-58°C (greenhouse effect). Such coverage, applied for a few hours during 1-3 consecutive days, had a detrimental effect on the survival of P. belbahrii: killing the pathogen and/or suppressing disease progress while enhancing growth of the host basil plants.     

Species delimitation in downy mildews: the case of Hyaloperonospora in the light of nuclear ribosomal ITS and LSU sequences

Species definitions for plant pathogens have considerable practical impact for measures such as plant protection or biological control, and are also important for comparative studies involving model organisms. However, in many groups, the delimitation of species is a notoriously difficult taxonomic problem. This is particularly evident in the obligate biotrophic downy mildew genera (Peronosporaceae, Peronosporales, Oomycetes), which display a considerable diversity with respect to genetic distances and host plants, but are, for the most part, morphologically rather uniform. The recently established genus Hyaloperonospora is of particular biological interest because it shows an impressive radiation on virtually a single host family, Brassicaceae, and it contains the downy mildew parasite, Arabidopsis thaliana, of importance as a model organism. Based on the most comprehensive molecular sampling of specimens from a downy mildew genus to date, including various collections from different host species and geographic locations, we investigate the phylogenetic relationships of Hyaloperonospora by molecular analysis of the nuclear ribosomal ITS and LSU sequences. Phylogenetic trees were inferred with ML and MP from the combined dataset; partitioned Bremer support (PBrS) was used to assess potential conflict between data partitions. As in other downy mildew groups, the molecular data clearly corroborate earlier results that supported the use of narrow species delimitations and host ranges as taxonomic markers. With few exceptions, suggested species boundaries are supported without conflict between different data partitions. The results indicate that a combination of molecular and host features is a reliable means to discriminate downy mildew species for which morphological differences are unknown.     

Identification of Genetic Loci Associated with Quality Traits in Almond via Association Mapping

To design an appropriate association study, we need to understand population structure and the structure of linkage disequilibrium within and among populations as well as in different regions of the genome in an organism. In this study, we have used a total of 98 almond accessions, from five continents located and maintained at the Centro de Investigación y Tecnología Agroalimentaria de Aragón (CITA; Spain), and 40 microsatellite markers. Population structure analysis performed in ‘Structure’ grouped the accessions into two principal groups; the Mediterranean (Western-Europe) and the non-Mediterranean, with K = 3, being the best fit for our data. There was a strong subpopulation structure with linkage disequilibrium decaying with increasing genetic distance resulting in lower levels of linkage disequilibrium between more distant markers. A significant impact of population structure on linkage disequilibrium in the almond cultivar groups was observed. The mean r² value for all intra-chromosomal loci pairs was 0.040, whereas, the r² for the inter-chromosomal loci pairs was 0.036. For analysis of association between the markers and phenotypic traits, five models comprising both general linear models and mixed linear models were selected to test the marker trait associations. The mixed linear model (MLM) approach using co-ancestry values from population structure and kinship estimates (K model) as covariates identified a maximum of 16 significant associations for chemical traits and 12 for physical traits. This study reports for the first time the use of association mapping for determining marker-locus trait associations in a world-wide almond germplasm collection. It is likely that association mapping will have the most immediate and largest impact on the tier of crops such as almond with the greatest economic value.     

Density- and Frequency-Dependent Selection at the Mdh-2 Locus in DROSOPHILA PSEUDOOBSCURA

We have studied differences in the number of Drosophila pseudoobscura produced in a culture when the flies differ with respect to two alleles (F and S) at the Mdh-2 locus, which codes for a malate dehydrogenase enzyme. The studies were done at low and at high density in two- and three-genotype combinations (S/S, F/F and S/F), with one-genotype cultures as controls.——Density affects the fitness of the Mdh-2 genotypes. Different genotypes are differently affected, and the genotype of the competitors also makes a difference on the fitness of a given genotype. When three genotypes are present in a culture, particularly at high density, intergenotypic competition is less intense than intragenotypic competition at several frequency combinations. That is, there is "overcompensation": the three genotypes together exploit the environmental resources better than one genotype alone.—The fitness of the genotypes is frequency dependent in both two-genotype and three-genotype combinations. An inverse relationship between frequency and fitness is observed at high density. This may lead to a stable polymorphism, because the fitness of a genotype increases as its frequency decreases.—Forty independent strains, sampled from a natural population, were used in the experiments. This ensures that more than 95% of the variation present in the genome in the natural population is also present is the experimental cultures. It also ensures that the genetic background of the Mdh-2 alleles is randomized in the same way as it is in nature. However, the possibility remains that Mdh-2 alleles in nature are nonrandomly associated with alleles at closely linked loci. If linkage disequilibrium is present in the experiments because it exists in nature, then the observed effects (such as frequency-dependent selection) would affect the Mdh-2 locus in nature as well.     

SNP detection and population structure evaluation of Salix gordejevii Y. L. Chang et Skv. in Hunshandake Sandland,Inner Mongolia,China

Single nucleotide polymorphisms (SNPs) are the most abundant and richest form of genomic polymorphism and, hence, are highly favorable markers for genetic map construction and genome-wide association studies. Based on the DNA specific-locus amplified fragment sequencing (SLAF-seq) for large-scale SNP detection, the genetic diversity and population structure of Salix gordejevii Y. L. Chang et Skv., a valuable sand-fixing shrub, was assessed in 199 accessions from 20 populations in Hunshandake Sandland of northern China. A total of 623.15 M reads resulted in 30.49 × sequencing depth on average and a mean Q30 of 95.70%, and 2,287,715 SNPs in 178,509 polymorphic SLAF tags were obtained. By discarding minor allele frequency > 0.05 and integrity > 0.8, a total of 93,600 SNPs were retained for population genetic analyses, which revealed that 199 individuals could be divided into six groups based on cross-validation errors. However, this grouping pattern did not match the geographical distribution, indicating that there is no apparent geographic barrier in the blank areas where S. gordejevii was not distributed in Hunshandake Sandland. In addition, the physical distance of linkage disequilibrium decay in the analyzed S. gordejevii individuals was 18.5 kb when r² = 0.1. The linkage disequilibrium decay distances for different chromosomes varied from 4.6 kb (chromosome 16) to 37.8 kb (chromosome 3). The obtained SNPs offer suitable marker resources for further genetic and genomic studies and will benefit S. gordejevii breeding programs.     

Simulating a base population in honey bee for molecular genetic studies

Background

Over the past years, reports have indicated that honey bee populations are declining and that infestation by an ecto-parasitic mite (Varroa destructor) is one of the main causes. Selective breeding of resistant bees can help to prevent losses due to the parasite, but it requires that a robust breeding program and genetic evaluation are implemented. Genomic selection has emerged as an important tool in animal breeding programs and simulation studies have shown that it yields more accurate breeding value estimates, higher genetic gain and low rates of inbreeding. Since genomic selection relies on marker data, simulations conducted on a genomic dataset are a pre-requisite before selection can be implemented. Although genomic datasets have been simulated in other species undergoing genetic evaluation, simulation of a genomic dataset specific to the honey bee is required since this species has a distinct genetic and reproductive biology. Our software program was aimed at constructing a base population by simulating a random mating honey bee population. A forward-time population simulation approach was applied since it allows modeling of genetic characteristics and reproductive behavior specific to the honey bee.

Results

Our software program yielded a genomic dataset for a base population in linkage disequilibrium. In addition, information was obtained on (1) the position of markers on each chromosome, (2) allele frequency, (3) χ² statistics for Hardy-Weinberg equilibrium, (4) a sorted list of markers with a minor allele frequency less than or equal to the input value, (5) average r² values of linkage disequilibrium between all simulated marker loci pair for all generations and (6) average r² value of linkage disequilibrium in the last generation for selected markers with the highest minor allele frequency.

Conclusion

We developed a software program that takes into account the genetic and reproductive biology specific to the honey bee and that can be used to constitute a genomic dataset compatible with the simulation studies necessary to optimize breeding programs. The source code together with an instruction file is freely accessible at http://msproteomics.org/Research/Misc/honeybeepopulationsimulator.html