Metapopulation structure of the specialized herbivore Macrosiphoniella tanacetaria (Homoptera,Aphididae)

We investigated population dynamics, genetic diversity and spatial structure in the aphid species Macrosiphoniella tanacetaria, a specialist herbivore feeding on tansy, Tanacetum vulgare. Tansy plants (genets) consist of many shoots (ramets), and genets are grouped in sites. Thus, aphids feeding on tansy can cluster at the level of ramets, genets and sites. We studied aphid population dynamics in 1997 and 2001 and found that within sites: (i). at any time, aphids used only a fraction of the available ramets and genets; (ii). at the level of ramets, most aphid colonies survived only one week; (iii). at the level of genets, mean survival time was less than 4 weeks; and (iv). colonization and extinction events occurred throughout the season. We sampled aphids in seven sites in the Alsace region, France (4-45 km apart) and two sites in Germany in 1999 to study genetic structure within and between populations. Genetic analyses using nine microsatellite loci showed that: (i). genotypic variability was high, (ii). none of the populations was in Hardy-Weinberg equilibrium, (iii). heterozygote deficits and linkage disequilibria were frequent, and (iv). all populations were genetically differentiated, even at a small geographical scale. Renewed sampling of the Alsace sites in 2001 showed that three populations had become extinct and significant genetic changes had occurred in the remaining four populations. The frequencies of extinction and colonization events at several spatial scales suggest a hierarchical metapopulation structure for M. tanacetaria. Frequent population turnover and drift are likely causes for the genetic differentiation of M. tanacetaria populations.     

Statistical analysis of amplified fragment length polymorphism data: a toolbox for molecular ecologists and evolutionists

Recently, the amplified fragment length polymorphism (AFLP) technique has gained a lot of popularity, and is now frequently applied to a wide variety of organisms. Technical specificities of the AFLP procedure have been well documented over the years, but there is on the contrary little or scattered information about the statistical analysis of AFLPs. In this review, we describe the various methods available to handle AFLP data, focusing on four research topics at the population or individual level of analysis: (i) assessment of genetic diversity; (ii) identification of population structure; (iii) identification of hybrid individuals; and (iv) detection of markers associated with phenotypes. Two kinds of analysis methods can be distinguished, depending on whether they are based on the direct study of band presences or absences in AFLP profiles ('band-based' methods), or on allelic frequencies estimated at each locus from these profiles ('allele frequency-based' methods). We investigate the characteristics and limitations of these statistical tools; finally, we appeal for a wider adoption of methodologies borrowed from other research fields, like for example those especially designed to deal with binary data.     

The multiple-locus symmetric fertility model

Selection due to variation in the fecundity among matings of genotypes with respect to many loci each with two alleles is studied. The fitness of a mating depends only on the genotypic distinction between homozygote and heterozygote at each locus in the two individuals, and differences among loci are allowed. This symmetric fertility model is therefore a generalization of the multiple-locus symmetric viability model. The phenomena seen in the two-locus symmetric fertility model generalize—e.g., the possibility of joint stability of equilibria with linkage equilibrium and with linkage disequilibrium, and the existence of different types of totally polymorphic equilibria with the gametic proportions in linkage equilibrium. The central equilibrium with genotypic frequencies in Hardy-Weinberg proportions and gametic frequencies in Robbins proportions exists for all symmetric fertility models. For some symmetric fertility regimes additional equilibria exist with gametic frequencies in linkage equilibrium and with genotypic frequencies in Hardy-Weinberg proportions at all except one locus. These equilibria may exist in the dioecious symmetric viability model, and then they will be locally stable. For free recombination the stable equilibria show linkage equilibrium, but several of these with different numbers of polymorphic loci may be stable simultaneously.     

A new theory of MHC evolution: beyond selection on the immune genes

The major histocompatibility complex (MHC) is a dense region of immune genes with high levels of polymorphism, which are arranged in haplotype blocks. Traditional models of balancing selection (i.e. overdominance and negative frequency dependence) were developed to study the population genetics of single genes. However, the MHC is a multigene family surrounded by linked (non-neutral) polymorphisms, and not all of its features are well explained by these models. For example, (i) the high levels of polymorphism in small populations, (ii) the unexpectedly large genetic differentiation between populations, (iii) the shape of the allelic genealogy associated with trans-species evolution, and (iv) the close associations between particular MHC (human leucocyte antigen, HLA) haplotypes and the approximately 100 pathologies in humans. Here, I propose a new model of MHC evolution named Associative Balancing Complex evolution that can explain these phenomena. The model proposes that recessive deleterious mutations accumulate as a 'sheltered load' nearby MHC genes. These mutations can accumulate because (i) they are rarely expressed as homozygotes given the high MHC gene diversity and (ii) purifying selection is inefficient with low recombination rates (cf. Muller's ratchet). Once fixed, these mutations add to balancing selection and further reinforce linkage through epistatic selection against recombinants.     

Accuracy of haplotype frequency estimation for biallelic loci, via the expectation-maximization algorithm for unphased diploid genotype data

Haplotype analyses have become increasingly common in genetic studies of human disease because of their ability to identify unique chromosomal segments likely to harbor disease-predisposing genes. The study of haplotypes is also used to investigate many population processes, such as migration and immigration rates, linkage-disequilibrium strength, and the relatedness of populations. Unfortunately, many haplotype-analysis methods require phase information that can be difficult to obtain from samples of nonhaploid species. There are, however, strategies for estimating haplotype frequencies from unphased diploid genotype data collected on a sample of individuals that make use of the expectation-maximization (EM) algorithm to overcome the missing phase information. The accuracy of such strategies, compared with other phase-determination methods, must be assessed before their use can be advocated. In this study, we consider and explore sources of error between EM-derived haplotype frequency estimates and their population parameters, noting that much of this error is due to sampling error, which is inherent in all studies, even when phase can be determined. In light of this, we focus on the additional error between haplotype frequencies within a sample data set and EM-derived haplotype frequency estimates incurred by the estimation procedure. We assess the accuracy of haplotype frequency estimation as a function of a number of factors, including sample size, number of loci studied, allele frequencies, and locus-specific allelic departures from Hardy-Weinberg and linkage equilibrium. We point out the relative impacts of sampling error and estimation error, calling attention to the pronounced accuracy of EM estimates once sampling error has been accounted for. We also suggest that many factors that may influence accuracy can be assessed empirically within a data set-a fact that can be used to create "diagnostics" that a user can turn to for assessing potential inaccuracies in estimation.     

STORM: software for testing hypotheses of relatedness and mating patterns

Storm is a software package that allows users to test a variety of hypotheses regarding patterns of relatedness and patterns of mate choice and/or mate compatibility within a population. These functions are based on four main calculations that can be conducted either independently or in the hypothesis-testing framework: internal relatedness; homozygosity by loci; pairwise relatedness; and a new metric called allele inheritance, which calculates the proportion of loci at which an offspring inherits a paternal allele different from that inherited from its mother. STORM allows users to test four hypotheses based on these calculations and Monte Carlo simulations: (i) are individuals within observed associations or groupings more/less related than expected; (ii) do observed offspring have more/less genetic variability (based on internal relatedness or homozygosity by loci) than expected from the gene pool; (iii) are observed mating pairs more/less related than expected if mating is random with respect to relatedness; and (iv) do observed offspring inherit paternal alleles different from those inherited from the mother more/less often than expected based on Mendelian inheritance.     

Genomics of natural populations: Evolutionary forces that establish and maintain gene arrangements in Drosophila pseudoobscura

The evolution of complex traits in heterogeneous environments may shape the order of genes within chromosomes. Drosophila pseudoobscura has a rich gene arrangement polymorphism that allows one to test evolutionary genetic hypotheses about how chromosomal inversions are established in populations. D. pseudoobscura has >30 gene arrangements on a single chromosome that were generated through a series of overlapping inversion mutations with >10 inversions with appreciable frequencies and wide geographic distributions. This study analyses the genomic sequences of 54 strains of Drosophila pseudoobscura that carry one of six different chromosomal arrangements to test whether (i) genetic drift, (ii) hitchhiking with an adaptive allele, (iii) direct effects of inversions to create gene disruptions caused by breakpoints, or (iv) indirect effects of inversions in limiting the formation of recombinant gametes are responsible for the establishment of new gene arrangements. We found that the inversion events do not disrupt the structure of protein coding genes at the breakpoints. Population genetic analyses of 2,669 protein coding genes identified 277 outlier loci harbouring elevated frequencies of arrangement‐specific derived alleles. Significant linkage disequilibrium occurs among distant loci interspersed between regions with low levels of association indicating that distant allelic combinations are held together despite shared polymorphism among arrangements. Outlier genes showing evidence of genetic differentiation between arrangements are enriched for sensory perception and detoxification genes. The data presented here support the indirect effect of inversion hypothesis where chromosomal inversions are favoured because they maintain linked associations among multilocus allelic combinations among different arrangements.     

Endocrine disruptors: present issues, future directions

A variety of natural products and synthetic chemicals, known collectively as endocrine-disrupting compounds (EDCs), mimic or interfere with the mechanisms that govern vertebrate reproductive development and function. At present, research has focused on (i) the morphological and functional consequences of EDCs; (ii) identifying and determining the relative potencies of synthetic and steroidal compounds that have endocrine-disrupting effects; (iii) the mechanism of action of EDCs at the molecular level; and (iv) the recognition that in "real life," contamination usually reflects mixtures of EDCs. Future research must examine (i) the interactive nature of EDCs, particularly whether the threshold concept as developed in traditional toxicological research applies to these chemicals; (ii) when and how EDCs act at the physiological level, particularly how they may organize the neural substrates of reproductive physiology and behavior; (iii) the various effects these compounds have on different species, individuals, and even tissues; and (iv) how adaptations may evolve in natural populations with continued exposure to EDCs. Several predictions are offered that reflect these new perspectives. Specifically, (i) the threshold assumption will be found not to apply to EDCs because they mimic the actions of endogenous molecules (e.g., estrogen) critical to development; hence, the threshold is automatically exceeded with exposure. (ii) Behavior can compound and magnify the effects of EDCs over successive generations; that is, bioaccumulated EDCs inherited from the mother not only influence the morphological and physiological development of the offspring but also the offsprings' reproductive behavior as adults. This adult behavior, in turn, can have further consequences on the sexual development of their own young. (iii) The sensitivity of a species or an individual to a compound is related to species (individual)-typical concentrations of circulating gonadal steroid hormones. Related to this is the recent finding that alternate forms of the putative receptors are differentially distributed, thereby contributing to the different effects that have been observed. (iv) Except in extraordinary situations, populations often continue to exist in contaminated sites. One possible explanation for this observation that needs to be considered is that animals can rapidly adapt to the nature and level of contamination in their environment. It is unlikely that successive generations coincidentally become insensitive to gonadal steroid hormones fundamentally important as biological regulators of development and reproduction. Rather, adaptive alterations in the genes that encode steroid receptors may occur with chronic exposure to EDCs, allowing the sex hormone receptor to discriminate natural steroids from EDCs.     

An adaptive explanation for geographically structured allozyme variation inDichroplus elongatus (Orthoptera: Acrididae)

Previous studies on allozyme variation in five populations of the grasshopperDichroplus elongatus along a geographical gradient in Argentina revealed a significant degree of population structuring and a significant association between one of the loci (Aat-1) and latitude. As this could not be entirely explained by historical factors, the possible adaptive significance of this locus or loci in linkage disequilibrium was investigated in the present study, with an emphasis on the role of environmental variables correlated with latitude. The present paper reports a study of the geographical organization of allozyme diversity over a wider range of the distribution ofD. elongatus. The relation of allelic frequencies with geographic climatic variables was analysed. We have found (i) that different loci (Aat-1 andPep-1) covary significantly with different variables, (ii) discordance between genetic and geographical distances, (iii) copious gene flow that would mask allelic frequency differences due solely to genetic drift, (iv) and temporal stability of the gradient found. The results suggest that causes other than drift and migration may explain the observed directional patterns of variation inD. elongatus.     

Identification of kin structure among Guam rail founders: a comparison of pedigrees and DNA profiles

Kin structure among founders can have a significant effect on subsequent population structure. Here we use the correlation between DNA profile similarity and relatedness calculated from pedigrees to test hypotheses regarding kin structure among founders to the captive Guam rail (Rallus owstoni) population. Five different pedigrees were generated under the following hypotheses: (i) founders are unrelated; (ii) founders are unrelated except for same-nest chicks; (iii) founders from the same major site are siblings; (iv) founders from the same local site are siblings; and (v) founders are related as defined by a UPGMA cluster analysis of DNA similarity data. Relatedness values from pedigrees 1, 2 and 5 had the highest correlation with DNA similarity but the correlation between relatedness and similarity were not significantly different among pedigrees. Pedigree 5 resulted in the highest correlation overall when using only relatedness values that changed as a result of different founder hypotheses. Thus, founders were assigned relatedness based on pedigree 5 because it had the highest correlations with DNA similarity, was the most conservative approach, and incorporated all field data. The analyses indicated that estimating relatedness using DNA profiles remains problematic, therefore we compared mean kinship, a measure of genetic importance, with mean DNA profile similarity to determine if genetic importance among individuals could be determined via use of DNA profiles alone. The significant correlation suggests this method may provide more information about population structure than was previously thought. Thus, DNA profiles can provide a reasonable explanation for founder relatedness and mean DNA profile similarity may be helpful in determining relative genetic importance of individuals when detailed pedigrees are absent.     

Survival quantitative trait locus fine mapping by measuring and testing for Hardy-Weinberg and linkage disequilibrium

I show that fine-scale localization of a survival-related locus can be accomplished on the basis of deviations from Hardy-Weinberg equilibrium and linkage disequilibrium at closely linked marker loci. The method is based on chi(2)-tests and they can be performed for age-specific samples of alive (or dead) individuals, as for combined samples of alive and dead individuals.     

Cooperation,clumping and the evolution of multicellularity

The evolution of multicellular organisms represents one of the major evolutionary transitions in the history of life. A potential advantage of forming multicellular clumps is that it provides an efficiency benefit to pre-existing cooperation, such as the production of extracellular ‘public goods’. However, this is complicated by the fact that cooperation could jointly evolve with clumping, and clumping could have multiple consequences for the evolution of cooperation. We model the evolution of clumping and a cooperative public good, showing that (i) when considered separately, both clumping and public goods production gradually increase with increasing genetic relatedness; (ii) in contrast, when the traits evolve jointly, a small increase in relatedness can lead to a major shift in evolutionary outcome—from a non-clumping state with low public goods production to a cooperative clumping state with high values of both traits; (iii) high relatedness makes it easier to get to the cooperative clumping state and (iv) clumping can be inhibited when it increases the number of cells that the benefits of cooperation must be shared with, but promoted when it increases relatedness between those cells. Overall, our results suggest that public goods sharing can facilitate the formation of well-integrated cooperative clumps as a first step in the evolution of multicellularity.     

6R instrumented spatial linkages for anatomical joint motion measurement--Part 1: Design.

Six-revolute-joint instrumented spatial linkages (6R ISLs) have become often-used devices to measure the complete six-degree-of-freedom motion of anatomical joints. Accuracy of motion measurement depends on ISL design and calibration technique. In this paper, a design process is outlined that uses computer graphics and numerical methods as aids in developing 6R ISLs that (i) physically assemble within the desired range of motion of the joint; (ii) do not collide with either the experimental apparatus or the subject joint; (iii) avoid singular linkage configurations that can cause forces to be applied to the joint; and (iv) measure selected anatomical motions most accurately. It is found that a certain subgroup of 6R linkages are suitable for accurate measurement of specific motions, and can be the basis for new ISL designs. General guidelines are developed that can assist in the generation of unique linkage designs for different anatomical joints. The design process is demonstrated in the creation of an ISL to measure knee motion.     

spagedi: a versatile computer program to analyse spatial genetic structure at the individual or population levels

spag e d i version 1.0 is a software primarily designed to characterize the spatial genetic structure of mapped individuals or populations using genotype data of codominant markers. It computes various statistics describing genetic relatedness or differentiation between individuals or populations by pairwise comparisons and tests their significance by appropriate numerical resampling. spag e d i is useful for: (i) detecting isolation by distance within or among populations and estimating gene dispersal parameters; (ii) assessing genetic relatedness between individuals and its actual variance, a parameter of interest for marker based inferences of quantitative inheritance; (iii) assessing genetic differentiation among populations, including the case of haploids or autopolyploids.     

PERMANENT GENETIC RESOURCES: Isolation and characterization of polymorphic microsatellite loci in common evening primrose (Oenothera biennis)

We developed nine polymorphic microsatellite loci for evening primrose (Oenothera biennis). These loci have two to 18 alleles per locus and observed heterozygosities ranging from 0 to 0.879 in a sample of 34 individuals. In a pattern consistent with the functionally asexual reproductive system of this species, 17/36 pairs of loci revealed significant linkage disequilibrium and three loci showed significant deviations from Hardy-Weinberg equilibrium. The loci will be informative in identifying genotypes in multigenerational field studies to assess changes in genotype frequencies.     

Evaluation of protective potentials of a potentized homeopathic drug, Chelidonium majus, during azo dye induced hepatocarcinogenesis in mice

Several cytogenetical and enzymatic protocols were used to test if two microdoses of Chelidonium majus, namely Chelidonium-30 (Ch-30) and Chelidonium-200 (Ch-200), used as homeopathic drugs, showed anti-tumor activity and also favorably modulated genotoxic damages produced by an azo dye in mice at several intervals of fixation. Different sets of healthy mice were fed: (i) hepatocarcinogen, p-dimethylaminoazobenzene (p-DAB, initiator) + phenobarbital (PB, promoter), (ii) only p-DAB, (iii) only PB, and (iv) neither p-DAB nor PB (normal control). Mice fed with p-DAB + PB were divided into different sets that were also fed either Ch-30 (v) or Ch-200 (vi) or diluted alcohol (vii), the "vehicle" of the microdoses of Chelidonium. All mice of group (i), a few of group (ii) and group (vii) and none of groups (iii) and (iv) developed tumors in liver at the longer intervals of fixation. The frequencies of chromosome aberrations (CA), micronucleated erythrocytes (MN), mitotic index (MI) and sperm head abnormality (SHA) were much higher in groups (i) and (vii) mice than in groups (ii), (iii) and (iv) mice at all fixation intervals. However, in mice of both groups (v) and (vi), the frequencies of CA, MN, SHA were strikingly less than those of groups (i) and (vii), and moderately less than those of groups (ii) and (iii). Both Ch-30 and Ch-200 also modulated favourably some toxicity marker enzymes like acid and alkaline phosphatases, peroxidases, glutamate oxaloacetate and glutamate pyruvate transaminases in liver, kidney and spleen tissues of the carcinogen fed mice. The microdoses of Chelidonium having no visible ill effects of their own, may be strong candidates for use in delaying/protecting liver cancer.     

Genetic Mapping in Escherichia coli K-12 by Radiation-Induced Crossing-Over

Conventional methods for chromosomal mapping in Escherichia coli are (i) interruption of matings to obtain minimum marker entry times, (ii) linkage analysis of recombinants, and (iii) cotransduction. Method (i) has a resolution of about 0.5 min (5 x 10(4) nucleotides) and is not useful for distances less than about 1 min; methods (ii) and (iii) are capable of better resolution but are generally not very reproducible and no general theory is available for translating crossing-over and cotransduction frequencies into physical chromosomal distances. We found that when merozygotes are irradiated (X rays or ultraviolet light) soon after marker transfer, high linkage values (0.8 to 1.0) between nearby marker pairs decrease with radiation dose to 0.5. Our results are quantitatively consistent with the idea that radiations induce crossing-over lesions proportional to dose, and the number of such lesions between two markers is proportional to the physical separation of the markers in the range that can also be measured by interruption of mating (0.5 to 4.0 min). Additivity relations among markers are also satisfied. We used this technique to measure the distances (0.1 to 1.0 min) between several pairs of closely linked markers.     

Polymorphic microsatellite loci for the cardinal fish (<Emphasis Type="Italic">Apogon imberbis</Emphasis>)

Eight polymorphic microsatellite loci were isolated and characterized for the Cardinal fish (Apogon imberbis), a coastal-reef fish endemic to the Mediterranean Sea. Characterization of 30 Cardinal fish individuals form the western Mediterranean showed moderate to high allelic diversity ranging from 6 to 19 alleles per locus. Two loci showed significant departures from Hardy-Weinberg equilibrium presumably due to null alleles. No evidence of linkage disequilibrium was found for any locus pairwise comparasions. This microsatellite set could be useful for any basic population genetic studies of this species.     

Assessing population genetic structure via the maximisation of genetic distance

Background

The inference of the hidden structure of a population is an essential issue in population genetics. Recently, several methods have been proposed to infer population structure in population genetics.

Methods

In this study, a new method to infer the number of clusters and to assign individuals to the inferred populations is proposed. This approach does not make any assumption on Hardy-Weinberg and linkage equilibrium. The implemented criterion is the maximisation (via a simulated annealing algorithm) of the averaged genetic distance between a predefined number of clusters. The performance of this method is compared with two Bayesian approaches: STRUCTURE and BAPS, using simulated data and also a real human data set.

Results

The simulations show that with a reduced number of markers, BAPS overestimates the number of clusters and presents a reduced proportion of correct groupings. The accuracy of the new method is approximately the same as for STRUCTURE. Also, in Hardy-Weinberg and linkage disequilibrium cases, BAPS performs incorrectly. In these situations, STRUCTURE and the new method show an equivalent behaviour with respect to the number of inferred clusters, although the proportion of correct groupings is slightly better with the new method. Re-establishing equilibrium with the randomisation procedures improves the precision of the Bayesian approaches. All methods have a good precision for F_ST ≥ 0.03, but only STRUCTURE estimates the correct number of clusters for F_ST as low as 0.01. In situations with a high number of clusters or a more complex population structure, MGD performs better than STRUCTURE and BAPS. The results for a human data set analysed with the new method are congruent with the geographical regions previously found.

Conclusion

This new method used to infer the hidden structure in a population, based on the maximisation of the genetic distance and not taking into consideration any assumption about Hardy-Weinberg and linkage equilibrium, performs well under different simulated scenarios and with real data. Therefore, it could be a useful tool to determine genetically homogeneous groups, especially in those situations where the number of clusters is high, with complex population structure and where Hardy-Weinberg and/or linkage equilibrium are present.     

"Fundamental theorem" for two Loci

THE CHANGE IN MEAN FITNESS FOR TWO LOCI UNDER SELECTION CAN BE DESCRIBED BY FOUR TERMS: (i) the variance of fitness, (ii) a weighted between-gamete covariance, (iii) a function of recombination, linkage disequilibrium and the slope of the surface of mean fitness on disequilibrium, and (iv) a function of these two parameters and the curvature of the surface. Independent derivations of this equation by different methods, although thought at one time to be in disagreement, give algebraically identical results.