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     

Analysis of arthropod bloodmeals using molecular genetic markers

Little is known about the transmission dynamics of human malaria and other vector-borne diseases, partly because of the limited availability and distribution of appropriate tools for quantifying human-mosquito contact rates. Recent developments in molecular biology have allowed a significant increase in the efficacy and reliability of bloodmeal identification, and DNA-based molecular markers are now being harnessed for typing arthropod bloodmeals. The extent to which these markers have been used for analysis of mosquito bloodmeals and the potential they might have for the future is discussed, and the contributions that the advent of PCR has made are examined here.     

Containment of arthropod disease vectors

Effective containment of arthropod vectors of infectious diseases is necessary to prevent transmission of pathogens by released, infected vectors and to prevent vectors that escape from establishing populations that subsequently contribute to increased disease. Although rare, past releases illustrate what can go wrong and justify the need for guidelines that minimize risks. An overview of recommendations for insectary facilities, practices, and equipment is provided, and features of four recently published and increasingly rigorous arthropod containment levels (ACLs 1-4) are summarized. ACL-1 is appropriate for research that constitutes the lowest risk level, including uninfected arthropods or vectors that are infected with micro-organisms that do not cause disease in humans, domestic animals, or wildlife. ACL-2 is appropriate for indigenous and exotic arthropods that represent a moderate risk, including vectors infected or suspected of being infected with biosafety level (BSL)-2 infectious agents and arthropods that have been genetically modified in ways that do not significantly affect their fecundity, survival, host preference, or vector competence. ACL-3 is recommended for arthropods that are or may be infected with BSL-3 infectious agents. ACL-3 places greater emphasis on pathogen containment and more restricted access to the insectary than ACL-2. ACL-4 is intended for arthropods that are infected with the most dangerous BSL-4 infectious agents, which can cause life-threatening illness by aerosol or arthropod bite. Adherence to these guidelines will result in laboratory-based arthropod vector research that minimizes risks and results in important new contributions to applied and basic science.     

Polymorphism of untranscribed spacer rDNA as a molecular genetic marker in population studies of cattle]

Variable polymorphic patterns were detected using EcoRI-SalI fragment of bovine rDNA, including 3'-end of 28S rRNA gene with the adjacent portion of the transcribed spacer, as a probe for hybridization. Some features of these polymorphic patterns are similar to DNA fingerprints detected with the M13 probe. Bovine rDNA spacer polymorphism was used as a molecular genetic marker for population analysis of individual specific patterns of 4 cattle breeds with the help of the Jeffreys' method. It was supposed that the probability of identical fingerprints appearance could be the characteristics of heterogeneity of cattle populations. The observed length of polymorphic gragments ranged from 2000 to 6000 bp. The mean number of fragments per individual for all breeds was 15.05. The probability of identical patterns appearance was very high: from 1.18 x 10(-5) in ajshir's breed to 1.43 x 10(-7) in "white and black"s' breed. So, high probability seems to be dependent on the high allelic frequency and the way of breeding.     

Use of DNA fingerprinting for human population genetic studies

DNA fingerprinting techniques have been used in population genetic studies on many different kinds of organisms. Here, we present new applications for multilocus DNA fingerprint probes in population studies and demonstrate the applicability of DNA fingerprinting to human population genetics, using M13 phage DNA as a probe. The new approach, which is based on a factor method of numerical coding of non-quantitative data (factor correspondence analysis-FCA), shows good agreement between population position, as indicated by the three principal factors, and ethnogenetic proximity.     

Fundamentals of fungal molecular population genetic analyses

The last two decades have seen tremendous growth in the development and application of molecular methods in the analyses of fungal species and populations. In this paper, I provide an overview of the molecular techniques and the basic analytical tools used to address various fundamental population and evolutionary genetic questions in fungi. With increasing availability and decreasing cost, DNA sequencing is becoming a mainstream data acquisition method in fungal evolutionary genetic studies. However, other methods, especially those based on the polymerase chain reaction, remain powerful in addressing specific questions for certain groups of taxa. These developments are bringing fungal population and evolutionary genetics into mainstream ecology and evolutionary biology.     

Environmental information systems for the control of arthropod vectors of disease

Over the last decade, remote sensing technologies and geographical information systems have moved from the research arena into the hands of vector control specialists. This review explains remote sensing approaches and spatial information technologies used for investigations of arthropod pests and vectors of diseases affecting humans and livestock. Relevant applications are summarized with examples of studies on African horse sickness vector Culicoides midges (Diptera: Ceratopogonidae), malaria vector Anopheles and arbovirus vector culicine mosquitoes (Diptera: Culicidae), leishmaniasis vector Phlebotomus sandflies (Diptera: Psychodidae), trypanosomiasis vector tsetse (Diptera: Glossinidae), loaiasis vector Chrysops (Diptera: Tabanidae), Lyme disease vector Ixodes and other ticks (Acari: Ixodidae). Methods and their uses are tabulated and discussed with recommendations for efficiency, caution and progress in this burgeoning field.     

Making sense of genetic estimates of effective population size

The last decade has seen an explosion of interest in use of genetic markers to estimate effective population size, N_e. Effective population size is important both theoretically (N_e is a key parameter in almost every aspect of evolutionary biology) and for practical application (N_e determines rates of genetic drift and loss of genetic variability and modulates the effectiveness of selection, so it is crucial to consider in conservation). As documented by Palstra & Fraser ( 2012 ), most of the recent growth in N_e estimation can be attributed to development or refinement of methods that can use a single sample of individuals (the older temporal method requires at least two samples separated in time). As with other population genetic methods, performance of new N_e estimators is typically evaluated with simulated data for a few scenarios selected by the author(s). Inevitably, these initial evaluations fail to fully consider the consequences of violating simplifying assumptions, such as discrete generations, closed populations of constant size and selective neutrality. Subsequently, many researchers studying natural or captive populations have reported estimates of N_e for multiple methods; often these estimates are congruent, but that is not always the case. Because true N_e is rarely known in these empirical studies, it is difficult to make sense of the results when estimates differ substantially among methods. What is needed is a rigorous, comparative analysis under realistic scenarios for which true N_e is known. Recently, Gilbert & Whitlock ( 2015 ) did just that for both single‐sample and temporal methods under a wide range of migration schemes. In the current issue of Molecular Ecology, Wang ( 2016 ) uses simulations to evaluate performance of four single‐sample N_e estimators. In addition to assessing effects of true N_e, sample size, and number of loci, Wang also evaluated performance under changing abundance, physical linkage and genotyping errors, as well as for some alternative life histories (high rates of selfing; haplodiploids). Wang showed that the sibship frequency (SF) and linkage disequilibrium (LD) methods perform dramatically better than the heterozygote excess and molecular coancestry methods under most scenarios (see Fig. 1, modified from figure 2 in Wang 2016 ), and he also concluded that SF is generally more versatile than LD. This article represents a truly Herculean effort, and results should be of considerable value to researchers interested in applying these methods to real‐world situations.     

Laboratory containment practices for arthropod vectors of human and animal pathogens

Arthropod-borne pathogens have an impact on the health and well-being of humans and animals throughout the world. Research involving arthropod vectors of disease is often dependent on the ability to maintain the specific arthropod species in laboratory colonies. The author reviews current arthropod containment practices and discusses their importance from public health and ecological perspectives.     

The population genetic structure of vectors and our understanding of disease epidemiology

Understanding and predicting disease epidemiology relies on clear knowledge about the basic biology of the organisms involved. Despite the key role that arthropod vectors play in disease dynamics and detailed mechanistic work on the vector-pathogen interface, little information is often available about how these populations function under natural conditions. Population genetic studies can help fill this void by providing information about the taxonomic status of species, the spatial limits of populations, and the nature of gene flow among populations. Here, I briefly review different types of population genetic structure and some recent examples of where this information has provided key elements for understanding pathogen transmission in tick-borne systems.     

Molecular bases of proliferation of Francisella tularensis in arthropod vectors

Arthropod vectors are important vehicles for transmission of Francisella tularensis between mammals, but very little is known about the F. tularensis–arthropod vector interaction. Drosophila melanogaster has been recently developed as an arthropod vector model for F. tularensis. We have shown that intracellular trafficking of F. tularensis within human monocytes‐derived macrophages and D. melanogaster‐derived S2 cells is very similar. Within both evolutionarily distant host cells, the Francisella‐containing phagosome matures to a late endosome‐like phagosome with limited fusion to lysosomes followed by rapid bacterial escape into the cytosol where the bacterial proliferate. To decipher the molecular bases of intracellular proliferation of F. tularensis within arthropod‐derived cells, we screened a comprehensive library of mutants of F. tularensis ssp. novicida for their defect in intracellular proliferation within D. melanogaster‐derived S2 cells. Our data show that 394 genes, representing 22% of the genome, are required for intracellular proliferation within D. melanogaster‐derived S2 cells, including many of the Francisella Pathogenicity Island (FPI) genes that are also required for proliferation within mammalian macrophages. Functional gene classes that exhibit growth defect include metabolic (25%), FPI (2%), type IV pili (1%), transport (16%) and DNA modification (5%). Among 168 most defective mutants in intracellular proliferation in S2 cells, 80 are defective in lethality and proliferation within adult D. melanogaster. The observation that only 135 of the 394 mutants that are defective in S2 cells are also defective in human macrophages indicates that F. tularensis utilize common as well as distinct mechanisms to proliferate within mammalian and arthropod cells. Our studies will facilitate deciphering the molecular aspects of F. tularensis–arthropod vector interaction and its patho‐adaptation to infect mammals.     

Effects of population structure on genetic association studies

Population-based case-control association is a promising approach for unravelling the genetic basis of complex diseases. One potential problem of this approach is the presence of population structure in the samples. Using the Collaborative Study on the Genetics of Alcoholism (COGA) single-nucleotide polymorphism (SNP) datasets, we addressed three questions: How can the degree of population structure be quantified, and how does the population structure affect association studies? How accurate and efficient is the genomic control method in correcting for population structure? The amount of population structure in the COGA SNP data was found to inflate the p-value in association tests. Genomic control was found to be effective only when the appropriate number of markers was used in the control group in order to correctly calibrate the test. The approach presented in this paper could be used to select the appropriate number of markers for use in the genomic control method of correcting population structure.     

Faeces as a source of DNA for molecular studies in a threatened population of great bustards

With recent advances in molecular biology, it is now possible to use the trace amounts of DNA in faeces to non-invasively sample endangered species for genetic studies. A highly vulnerable population of approximately 100 great bustards (Otis tarda) exists in Morocco necessitating the use of non-invasive protocols to study their genetic structure. Here we report a reliable silica-based method to extract DNA from great bustard faeces. We found that successful extraction and amplification correlated strongly with faeces freshness and composition. We could not extract amplifiable DNA from 30% of our samples as they were dry or contained insect material. However 100% of our fresh faecal samples containing no obvious insect material worked, allowing us to assess the levels of genetic variation among 25 individuals using a 542 bp control region sequence. We were able to extract DNA from four out of five other avian species, demonstrating that faeces represents a suitable source of DNA for population genetics studies in a broad range of species. This revised version was published online in July 2006 with corrections to the Cover Date.