Measuring population differentiation using GST or D? A simulation study with microsatellite DNA markers under a finite island model and nonequilibrium conditions

The genetic differentiation of populations is a key parameter in population genetic investigations. Wright's F(ST) (and its relatives such as G(ST) ) has been a standard measure of differentiation. However, the deficiencies of these indexes have been increasingly realized in recent years, leading to some new measures being proposed, such as Jost's D (Molecular Ecology, 2008; 17, 4015). The existence of these new metrics has stimulated considerable debate and induced some confusion on which statistics should be used for estimating population differentiation. Here, we report a simulation study with neutral microsatellite DNA loci under a finite island model to compare the performance of G(ST) and D, particularly under nonequilibrium conditions. Our results suggest that there exist fundamental differences between the two statistics, and neither G(ST) nor D operates satisfactorily in all situations for quantifying differentiation. D is very sensitive to mutation models but G(ST) noticeably less so, which limits D's utility in population parameter estimation and comparisons across genetic markers. Also, the initial heterozygosity of the starting populations has some important effects on both the individual behaviours of G(ST) and D and their relative behaviours in early differentiation, and this effect is much greater for D than G(ST) . In the early stages of differentiation, when initial heterozygosity is relatively low (<0.5, if the number of subpopulations is large), G(ST) increases faster than D; the opposite is true when initial heterozygosity is high. Therefore, the state of the ancestral population appears to have some lasting impacts on population differentiation. In general, G(ST) can measure differentiation fairly well when heterozygosity is low whatever the causes; however, when heterozygosity is high (e.g. as a result of either high mutation rate or high initial heterozygosity) and gene flow is moderate to strong, G(ST) fails to measure differentiation. Interestingly, when population size is not very small (e.g. N ≥ 1000), G(ST) measures differentiation quite linearly with time over a long duration when gene flow is absent or very weak even if mutation rate is not low (e.g. μ = 0.001). In contrast, D, as a differentiation measure, performs rather robustly in all these situations. In practice, both indexes should be calculated and the relative levels of heterozygosities (especially H(S) ) and gene flow taken into account. We suggest that a comparison of the two indexes can generate useful insights into the evolutionary processes that influence population differentiation.     

Allele-specific gene expression in a wild nonhuman primate population

Natural populations hold enormous potential for evolutionary genetic studies, especially when phenotypic, genetic and environmental data are all available on the same individuals. However, untangling the genotype-phenotype relationship in natural populations remains a major challenge. Here, we describe results of an investigation of one class of phenotype, allele-specific gene expression (ASGE), in the well-studied natural population of baboons of the Amboseli basin, Kenya. ASGE measurements identify cases in which one allele of a gene is overexpressed relative to the alternative allele of the same gene, within individuals, thus providing a control for background genetic and environmental effects. Here, we characterize the incidence of ASGE in the Amboseli baboon population, focusing on the genetic and environmental contributions to ASGE in a set of eleven genes involved in immunity and defence. Within this set, we identify evidence for common ASGE in four genes. We also present examples of two relationships between cis-regulatory genetic variants and the ASGE phenotype. Finally, we identify one case in which this relationship is influenced by a novel gene-environment interaction. Specifically, the dominance rank of an individual's mother during its early life (an aspect of that individual's social environment) influences the expression of the gene CCL5 via an interaction with cis-regulatory genetic variation. These results illustrate how environmental and ecological data can be integrated into evolutionary genetic studies of functional variation in natural populations. They also highlight the potential importance of early life environmental variation in shaping the genetic architecture of complex traits in wild mammals.     

Generation of two modified mouse alleles of the Hic1 tumor suppressor gene

HIC1 (hypermethylated in cancer 1) is a tumor suppressor gene located on chromosome 17p13.3, a region frequently hypermethylated or deleted in human neoplasias. In mouse, Hic1 is essential for embryonic development and exerts an antitumor role in adult animals. Since Hic1-deficient mice die perinatally, we generated a conditional Hic1 null allele by flanking the Hic1-coding region by loxP sites. When crossed to animals expressing Cre recombinase in a cell-specific manner, the Hic1 conditional mice will provide new insights into the function of Hic1 in developing and mature tissues. Additionally, we used gene targeting to replace sequence-encoding amino acids 186-893 of Hic1 by citrine fluorescent protein cDNA. We demonstrate that the distribution of Hic1-citrine fusion polypeptide corresponds to the expression pattern of wild-type Hic1. Consequently, Hic1-citrine "reporter" mice can be used to monitor the activity of the Hic1 locus using citrine fluorescence.     

csd alleles in the red dwarf honey bee (Apis florea,Hymenoptera: Apidae) show exceptionally high nucleotide diversity

Abstract The single locus complementary sex determination (sl‐csd) gene is the primary gene determining the gender of honey bees (Apis spp.). While the csd gene has been well studied in the Western honey bee (Apis mellifera), and comparable data exist in both the Eastern honey bee (Apis cerana) and the giant honey bee (Apis dorsata), no studies have been conducted in the red dwarf honey bee, Apis florea. In this study we cloned the genomic region 3 of the A. florea csd gene from 60 workers, and identified 12 csd alleles. Analysis showed that similar to A. mellifera, region 3 of the csd gene contains a RS domain at the N terminal, a proline‐rich domain at the C terminal, and a hypervariable region in the middle. However, the A. florea csd gene possessed a much higher level of nucleotide diversity, compared to A. mellifera, A. cerana and Apis dorsata. We also show that similar to the other three Apis species, in A. florea, nonsynonymous mutations in the csd gene are selectively favored in young alleles.     

Investigation of serum protein systems in Chinese pigs

Polymorphisms of transferrin (Tf), pre-albumin (Pa), haemopexin (Hpx), ceruloplasmin (Cp) and amylase (Am) of Duroc pigs and Hunan indigenous pigs were investigated using horizontal starch gel electrophoresis. Allele frequencies of Durocs determined in 2004 were compared with frequencies presented in the paper by Baker L.N. (1968) Serum protein variation in Duroc and Hampshire pigs. Vox Sanguinis15, 154-8. The number of serum protein alleles decreased over time and allele frequencies aggregated across certain alleles, including TfB, PaA, Hpx3, CpB and AmB. Differences in allele frequencies, average heterozygosities and standard genetic distances between the Duroc pigs and Hunan indigenous pig populations were examined. The relationship between Durocs and Hunan indigenous pigs was found to be more distant than those among the three Hunan indigenous pig populations.     

Non-African origin of a local beneficial mutation in D. melanogaster

It is well understood that the out-of-Africa habitat expansion of D. melanogaster was associated with the fixation of many beneficial mutations. Nevertheless, it is not clear yet whether these beneficial mutations segregated already in Africa or originated outside of Africa. In this article, we describe an ongoing selective sweep specific to one European population. One microsatellite allele has increased in a population from The Netherlands to a frequency of 18%, whereas it is virtually absent in 12 other European populations. The selective sweep resulted in a genomic region of more than 600 kb that is identical by descent. This is probably the first evidence of a beneficial mutation that has arisen outside of Africa and has resulted in a selective sweep localized in a population from The Netherlands.     

Genetic structure in the coral-reef-associated Banggai cardinalfish, Pterapogon kauderni

In this study, we used 11 polymorphic microsatellite loci to show that oceanic distances as small as 2-5 km are sufficient to produce high levels of population genetic structure (multilocus F(ST) as high as 0.22) in the Banggai cardinalfish (Pterapogon kauderni), a heavily exploited reef fish lacking a pelagic larval dispersal phase. Global F(ST) among all populations, separated by a maximum distance of 203 km, was 0.18 (R(ST) = 0.35). Moreover, two lines of evidence suggest that estimates of F(ST) may actually underestimate the true level of genetic structure. First, within-locus F(ST) values were consistently close to the theoretical maximum set by the average within-population heterozygosity. Second, the allele size permutation test showed that R(ST) values were significantly larger than F(ST) values, indicating that populations have been isolated long enough for mutation to have played a role in generating allelic variation among populations. The high level of microspatial structure observed in this marine fish indicates that life history traits such as lack of pelagic larval phase and a good homing ability do indeed play a role in shaping population genetic structure in the marine realm.     

Interactome modeling

A long-term goal of the field of interactome modeling is to understand how global and local properties of complex macromolecular networks impact on observable biological properties, and how changes in such properties can lead to human diseases. The information available at this stage of development of the field provides strong evidence for the existence of such interesting global and local properties, but also demonstrates that many more datasets will be needed to provide accurate models with increasingly predictive capacity. This review focuses on an early attempt at mapping a multicellular interactome network and on the lessons learned from that attempt.     

Non-Mendelian transmission of alleles at microsatellite loci: an example in Ixodes ricinus, the vector of Lyme disease

Microsatellite loci are usually considered to be neutral co-dominant and Mendelian markers. We undertook to study the inheritance of five microsatellite loci in the European Lyme disease vector, the tick Ixodes ricinus. Only two loci appeared fully Mendelian while the three others displayed non-Mendelian patterns that highly frequent null alleles could not fully explain. At one locus, IR27, some phenomenon seems to hinder the PCR amplification of one allele, depending on its origin (maternal imprinting) and/or its size (short allele dominance). DNA methylation, which appeared to be a possible explanation of this amplification bias, was rejected by a specific test comparing the amplification efficiency that did not differ between unmethylated and experimentally methylated DNA. The role of allele size in heterozygous individuals was then revealed from the data available on field collected ticks and consistent with the results of a theoretical approach. These observations highlight the need for prudence while inferring reproductive systems (selfing rates), parentage or even allelic frequencies from microsatellite markers, in particular for parasitic organisms for which molecular approaches often represent the only way for population biology inferences.     

Allelic structure and distribution of 103 STR loci in a Southern Tunisian population

Genotypes of 103 short tandem repeat (STR) markers distributed at an average of 40 cM intervals throughout the genome were determined for 40 individuals from the village of BirEl Hfai (BEH). This village of approximately 31.000 individuals is localized in the south-west of Tunisia. The allele frequency distributions in BEH were compared with those obtained for individuals in the CEPH (Centre d’Etude du Polymorphisme Humain) data using a Kolmo-gorov-Smirnov two-sample test. Fourteen out of the 103 markers (13.2%) showed significant differences (P<0.05) in distribution between the two populations. Population heterogeneity in BEH was indicated by an excess of observed homozygosity deviations from Hardy-Weinberg equilibrium at 3 loci (P<0.0005). No evidence for genotypic disequilibrium was found for any of the marker pairs. This demonstrated that in spite of a high inbreeding level in the population, few markers showed evidence for a different pattern of allelic distribution compared to CEPH.