Statistical methods for examining genetic influences of resistance to anti-epileptic drugs

It is generally accepted that genetics may be an important factor in explaining the variation between patients’ responses to certain drugs. However, identification and confirmation of the responsible genetic variants is proving to be a challenge in many cases. A number of difficulties that maybe encountered in pursuit of these variants, such as non-replication of a true effect, population structure and selection bias, can be mitigated or at least reduced by appropriate statistical methodology. Another major statistical challenge facing pharmacogenetics studies is trying to detect possibly small polygenic effects using large volumes of genetic data, while controlling the number of false positive signals. Here we review statistical design and analysis options available for investigations of genetic resistance to anti-epileptic drugs.

A novel parallel approach to the likelihood-based estimation of admixture in population genetics

Inferring population admixture from genetic data and quantifying it is a difficult but crucial task in evolutionary and conservation biology. Unfortunately state-of-the-art probabilistic approaches are computationally demanding. Effectively exploiting the computational power of modern multiprocessor systems can thus have a positive impact to Monte Carlo-based simulation of admixture modeling. A novel parallel approach is briefly described and promising results on its message passing interface (MPI)-based C implementation are reported.

Determination of mitochondrial genetic diversity in mammals

Mitochondrial DNA (mtDNA) is one of the most Popular population genetic markers. Its relevance as an indicator Of Population size and history has recently been questioned by several large-scale studies in animals reporting evidence for recurrent adaptive evolution, at least in invertebrates. Here we focus on mammals, a more restricted taxonomic group for which the issue of mtDNA near neutrality is crucial. By analyzing the distribution of mtDNA diversity across species and relating 4 to allozyme diversity, life-history traits, and taxonomy, we show that (i) mtDNA in mammals (toes not reject the nearly neutral model; (ii) mtDNA diversity, however, is unrelated to any of the 14 life-history and ecological variables that we analyzed, including body mass, geographic range, and The World Conservation Union (IUCN) categorization; (iii) mtDNA diversity is highly variable between mammalian orders and families; (iv) this taxonomic effect is most likely explained by variations of mutation rate between lineages. These results are indicative of a strong stochasticity of effective population size in mammalian species. They Suggest that, even in the absence of selection, mtDNA genetic diversity is essentially unpredictable, knowing species biology, and probably uncorrelated to species abundance.

Simple sequence repeats reveal uneven distribution of genetic diversity in chloroplast genomes of Brassica oleracea L. and (n=9) wild relatives

Diversity in the chloroplast genome of 171 accessions representing the Brassica 'C' (n = 9) genome, including domesticated and wild B. oleracea and nine inter-fertile related wild species, was investigated using six chloroplast SSR (microsatellite) markers. The lack of diversity detected among 105 cultivated and wild accessions of B. oleracea contrasted starkly with that found within its wild relatives. The vast majority of B. oleracea accessions shared a single haplotype, whereas as many as six haplotypes were detected in two wild species, B. villosa Biv. and B. cretica Lam.. The SSRs proved to be highly polymorphic across haplotypes, with calculated genetic diversity values (H) of 0.23-0.87. In total, 23 different haplotypes were detected in C genome species, with an additional five haplotypes detected in B. rapa L. (A genome n = 10) and another in B. nigra L. (B genome, n = 8). The low chloroplast diversity of B. oleracea is not suggestive of multiple domestication events. The predominant B. oleracea haplotype was also common in B. incana Ten. and present in low frequencies in B. villosa, B. macrocarpa Guss, B. rupestris Raf. and B. cretica. The chloroplast SSRs reveal a wealth of diversity within wild Brassica species that will facilitate further evolutionary and phylogeographic studies of this important crop genus.

Population genetic models can be used to study the evolution of the interacting behaviours of parents and their progeny

An example of the evolution of the interacting behaviours of parents and progeny is studied using iterative equations linking the frequencies of the gametes produced by the progeny to the frequencies of the gametes in the parental generation. This population genetics approach shows that a model in which both behaviours are determined by a single locus can lead to a stable equilibrium in which the two behaviours continue to segregate. A model in which the behaviours are determined by genes at two separate loci leads eventually to fixation of the alleles at both loci but this can take many generations of selection. Models of the type described in this paper will be needed to understand the evolution of complex behaviour when genomic or experimental information is available about the genetic determinants of behaviour and the selective values of different genomes. (c) 2007 Elsevier Inc. All rights reserved.

Domestication of plants in the americas: Insights from mendelian and molecular genetics

Background Plant domestication occurred independently in four different regions of the Americas. In general, different species were domesticated in each area, though a few species were domesticated independently in more than one area. The changes resulting from human selection conform to the familiar domestication syndrome, though different traits making up this syndrome, for example loss of dispersal, are achieved by different routes in crops belonging to different families. Genetic and Molecular Analyses of Domestication Understanding of the genetic control of elements of the domestication syndrome is improving as a result of the development of saturated linkage maps for major crops, identification and mapping of quantitative trait loci, cloning and sequencing of genes or parts of genes, and discoveries of widespread orthologies in genes and linkage groups within and between families. As the modes of action of the genes involved in domestication and the metabolic pathways leading to particular phenotypes become better understood, it should be possible to determine whether similar phenotypes have similar underlying genetic controls, or whether human selection in genetically related but independently domesticated taxa has fixed different mutants with similar phenotypic effects. Conclusions Such studies will permit more critical analysis of possible examples of multiple domestications and of the origin(s) and spread of distinctive variants within crops. They also offer the possibility of improving existing crops, not only major food staples but also minor crops that are potential export crops for developing countries or alternative crops for marginal areas.

The population genetic structure of clonal organisms generated by exponentially bounded and fat-tailed dispersal

Long distance dispersal (LDD) plays an important role in many population processes like colonization, range expansion, and epidemics. LDD of small particles like fungal spores is often a result of turbulent wind dispersal and is best described by functions with power-law behavior in the tails ("fat tailed"). The influence of fat-tailed LDD on population genetic structure is reported in this article. In computer simulations, the population structure generated by power-law dispersal with exponents in the range of -2 to -1, in distinct contrast to that generated by exponential dispersal, has a fractal structure. As the power-law exponent becomes smaller, the distribution of individual genotypes becomes more self-similar at different scales. Common statistics like G(ST) are not well suited to summarizing differences between the population genetic structures. Instead, fractal and self-similarity statistics demonstrated differences in structure arising from fat-tailed and exponential dispersal. When dispersal is fat tailed, a log-log plot of the Simpson index against distance between subpopulations has an approximately constant gradient over a large range of spatial scales. The fractal dimension D-2 is linearly inversely related to the power-law exponent, with a slope of similar to -2. In a large simulation arena, fat-tailed LDD allows colonization of the entire space by all genotypes whereas exponentially bounded dispersal eventually confines all descendants of a single clonal lineage to a relatively small area.

The genetic basis of resistance to anticoagulants in rodents

Anticoagulant compounds, i.e., derivatives of either 4-hydroxycoumarin (e.g., warfarin, bromadiolone) or indane-1,3-dione (e.g., diphacinone, chlorophacinone), have been in worldwide use as rodenticides for > 50 years. These compounds inhibit blood coagulation by repression of the vitamin K reductase reaction (VKOR). Anticoagulant-resistant rodent populations have been reported from many countries and pose a considerable problem for pest control. Resistance is transmitted as an autosomal dominant trait although, until recently, the basic genetic mutation was unknown. Here, we report on the identification of eight different mutations in the VKORC1 gene in resistant laboratory strains of brown rats and house mice and in wild-caught brown rats from various locations in Europe with five of these mutations affecting only two amino acids (Tyr139Cys, Tyr139Ser, Tyr139Phe and Leu128Gln, Leu128Ser). By recombinant expression of VKORC1 constructs in HEK293 cells we demonstrate that mutations at Tyr139 confer resistance to warlarin at variable degrees while the other mutations, in addition, dramatically reduce VKOR activity. Our data strongly argue for at least seven independent mutation events in brown rats and two in mice. They suggest that mutations in VKORC1 are the genetic basis of anticoagulant resistance in wild populations of rodents, although the mutations alone do not explain all aspects of resistance that have been reported. We hypothesize that these mutations, apart from generating structural changes in the VKORC1 protein, may induce compensatory mechanisms to maintain blood clotting. Our findings provide the basis for a DNA-based field monitoring of anticoagulant resistance in rodents.

Recent developments in genetic data analysis: what can they tell us about human demographic history?

Over the last decade, a number of new methods of population genetic analysis based on likelihood have been introduced. This review describes and explains the general statistical techniques that have recently been used, and discusses the underlying population genetic models. Experimental papers that use these methods to infer human demographic and phylogeographic history are reviewed. It appears that the use of likelihood has hitherto had little impact in the field of human population genetics, which is still primarily driven by more traditional approaches. However, with the current uncertainty about the effects of natural selection, population structure and ascertainment of single-nucleotide polymorphism markers, it is suggested that likelihood-based methods may have a greater impact in the future.

The Bayesian revolution in genetics

Bayesian statistics allow scientists to easily incorporate prior knowledge into their data analysis. Nonetheless, the sheer amount of computational power that is required for Bayesian statistical analyses has previously limited their use in genetics. These computational constraints have now largely been overcome and the underlying advantages of Bayesian approaches are putting them at the forefront of genetic data analysis in an increasing number of areas.

Identifying adaptive genetic divergence among populations from genome scans

The identification of signatures of natural selection in genomic surveys has become an area of intense research, stimulated by the increasing ease with which genetic markers can be typed. Loci identified as subject to selection may be functionally important, and hence (weak) candidates for involvement in disease causation. They can also be useful in determining the adaptive differentiation of populations, and exploring hypotheses about speciation. Adaptive differentiation has traditionally been identified from differences in allele frequencies among different populations, summarised by an estimate of F-ST. Low outliers relative to an appropriate neutral population-genetics model indicate loci subject to balancing selection, whereas high outliers suggest adaptive (directional) selection. However, the problem of identifying statistically significant departures from neutrality is complicated by confounding effects on the distribution of F-ST estimates, and current methods have not yet been tested in large-scale simulation experiments. Here, we simulate data from a structured population at many unlinked, diallelic loci that are predominantly neutral but with some loci subject to adaptive or balancing selection. We develop a hierarchical-Bayesian method, implemented via Markov chain Monte Carlo (MCMC), and assess its performance in distinguishing the loci simulated under selection from the neutral loci. We also compare this performance with that of a frequentist method, based on moment-based estimates of F-ST. We find that both methods can identify loci subject to adaptive selection when the selection coefficient is at least five times the migration rate. Neither method could reliably distinguish loci under balancing selection in our simulations, even when the selection coefficient is twenty times the migration rate.

A genetic linkage map of microsatellite, gene-specific and morphological markers in diploid Fragaria

Diploid Fragaria provide a potential model for genomic studies in the Rosaceae. To develop a genetic linkage map of diploid Fragaria, we scored 78 markers (68 microsatellites, one sequence-characterised amplified region, six gene-specific markers and three morphological traits) in an interspecific F2 population of 94 plants generated from a cross of F.vesca f. semperflorens × F. nubicola. Co-segregation analysis arranged 76 markers into seven discrete linkage groups covering 448 cM, with linkage group sizes ranging from 100.3 cM to 22.9 cM. Marker coverage was generally good; however some clustering of markers was observed on six of the seven linkage groups. Segregation distortion was observed at a high proportion of loci (54%), which could reflect the interspecific nature of the progeny and, in some cases, the self-incompatibility of F. nubicola. Such distortion may also account for some of the marker clustering observed in the map. One of the morphological markers, pale-green leaf (pg) has not previously been mapped in Fragaria and was located to the mid-point of linkage group VI. The transferable nature of the markers used in this study means that the map will be ideal for use as a framework for additional marker incorporation aimed at enhancing and resolving map coverage of the diploid Fragaria genome. The map also provides a sound basis for linkage map transfer to the cultivated octoploid strawberry.

The role played by the interaction between genetic factors and attachment in the stress response in infancy

Background: The importance of understanding which environmental and biological factors are involved in determining individual differences in physiological response to stress is widely recognized, given the impact that stress has on physical and mental health. Methods: The child-mother attachment relationship and some genetic polymorphisms (5-HTTLPR, COMT and GABRA6) were tested as predictors of salivary cortisol and alpha amylase concentrations, two biomarkers of hypothalamic-pituitary-adrenocortical (HPA) axis and sympathetic adrenomedullary (SAM) system activity, during the Strange Situation (SS) procedure in a sample of more than 100 healthy infants, aged 12 to 18 months. Results: Individual differences in alpha amylase response to separation were predicted by security of attachment in interaction with 5-HTTLPR and GABRA6 genetic polymorphisms, whereas alpha amylase basal levels were predicted by COMT x attachment interaction. No significant effect of attachment, genetics and their interaction on cortisol activity emerged. Conclusions: These results help to disentangle the role played by both genetic and environmental factors in determining individual differences in stress response in infancy. The results also shed light on the suggestion that HPA and SAM systems are likely to have different characteristic responses to stress.