Estimating Sampling Selection Bias in Human Genetics: A Phenomenological Approach

This research is the first empirical attempt to calculate the various components of the hidden bias associated with the sampling strategies routinely-used in human genetics, with special reference to surname-based strategies. We reconstructed surname distributions of 26 Italian communities with different demographic features across the last six centuries (years 1447–2001). The degree of overlapping between "reference founding core" distributions and the distributions obtained from sampling the present day communities by probabilistic and selective methods was quantified under different conditions and models. When taking into account only one individual per surname (low kinship model), the average discrepancy was 59.5%, with a peak of 84% by random sampling. When multiple individuals per surname were considered (high kinship model), the discrepancy decreased by 8–30% at the cost of a larger variance. Criteria aimed at maximizing locally-spread patrilineages and long-term residency appeared to be affected by recent gene flows much more than expected. Selection of the more frequent family names following low kinship criteria proved to be a suitable approach only for historically stable communities. In any other case true random sampling, despite its high variance, did not return more biased estimates than other selective methods. Our results indicate that the sampling of individuals bearing historically documented surnames (founders'' method) should be applied, especially when studying the male-specific genome, to prevent an over-stratification of ancient and recent genetic components that heavily biases inferences and statistics.     

Stabilizing Selection,Purifying Selection,and Mutational Bias in Finite Populations

Genomic traits such as codon usage and the lengths of noncoding sequences may be subject to stabilizing selection rather than purifying selection. Mutations affecting these traits are often biased in one direction. To investigate the potential role of stabilizing selection on genomic traits, the effects of mutational bias on the equilibrium value of a trait under stabilizing selection in a finite population were investigated, using two different mutational models. Numerical results were generated using a matrix method for calculating the probability distribution of variant frequencies at sites affecting the trait, as well as by Monte Carlo simulations. Analytical approximations were also derived, which provided useful insights into the numerical results. A novel conclusion is that the scaled intensity of selection acting on individual variants is nearly independent of the effective population size over a wide range of parameter space and is strongly determined by the logarithm of the mutational bias parameter. This is true even when there is a very small departure of the mean from the optimum, as is usually the case. This implies that studies of the frequency spectra of DNA sequence variants may be unable to distinguish between stabilizing and purifying selection. A similar investigation of purifying selection against deleterious mutations was also carried out. Contrary to previous suggestions, the scaled intensity of purifying selection with synergistic fitness effects is sensitive to population size, which is inconsistent with the general lack of sensitivity of codon usage to effective population size.     

出生在动物园的黑颈鹤

据上世纪80年代统计,全世界仅有800多只黑颈鹤,主要分布在中国。动物园显然不可以再直接从野外捕捉黑颈鹤补充自己的动物资源,可偏偏这宝贝在动物园里就是不肯繁殖。动手解决这个难题,拉开了动物园从利用自然转向保护自然的序幕。[编者按]     

Selection and Biased Gene Conversion in a Multigene Family: Consequences of Interallelic Bias and Threshold Selection

In a previous paper, I investigated the interactions in a gene family of additive selection and biased gene conversion in a finite population when conversion events are rare. Here I extend my "weak-conversion limit" model by allowing biased interallelic conversion (conversion between alleles at the same locus) of arbitrary frequency and various threshold selection schemes for rare interlocus conversion events. I suggest that it is not unreasonable for gene families to experience threshold fitness functions, and show that certain types of thresholds can greatly constrain the rate at which advantageous alleles are fixed as compared to other fitness schemes, such as additive selection. It is also shown that the double sampling process operating on a gene family in a finite population (sampling over the number of genes in the gene family and over the number of individuals in the population) can have interesting consequences. For selectively neutral alleles that experience interallelic bias, the probability of fixation at each single locus may be essentially neutral, but the cumulative effects on the entire gene family of small departures from neutrality can be significant, especially if the gene family is large. Thus, in some situations, gene families can respond to directional forces that are weak in comparison to drift at single loci.     

Self-Esteem Modulates the Time Course of Self-Positivity Bias in Explicit Self-Evaluation

Researchers have suggested that certain individuals may show a self-positivity bias, rating themselves as possessing more positive personality traits than others. Previous evidence has shown that people evaluate self-related information in such a way as to maintain or enhance self-esteem. However, whether self-esteem would modulate the time course of self-positivity bias in explicit self-evaluation has never been explored. In the present study, 21 participants completed the Rosenberg self-esteem scale and then completed a task where they were instructed to indicate to what extent positive/negative traits described themselves. Behavioral data showed that participants endorsed positive traits as higher in self-relevance compared to the negative traits. Further, participants’ self-esteem levels were positively correlated with their self-positivity bias. Electrophysiological data revealed smaller N1 amplitude and larger late positive component (LPC) amplitude to stimuli consistent with the self-positivity bias (positive-high self-relevant stimuli) when compared to stimuli that were inconsistent with the self-positivity bias (positive-low self-relevant stimuli). Moreover, only in individuals with low self-esteem, the latency of P2 was more pronounced in processing stimuli that were consistent with the self-positivity bias (negative-low self-relevant stimuli) than to stimuli that were inconsistent with the self-positivity bias (positive-low self-relevant stimuli). Overall, the present study provides additional support for the view that low self-esteem as a personality variable would affect the early attentional processing.     

The Effect of Selection Bias in Studies of Fads and Fashions

Most studies of fashion and fads focus on objects and practices that once were popular. We argue that limiting the sample to such trajectories generates a selection bias that obscures the underlying process and generates biased estimates. Through simulations and the analysis of a data set that has previously not been used to analyze the rise and fall of cultural practices, the New York Times text archive, we show that studying a whole range of cultural objects, both popular and less popular, is essential for understanding the drivers of popularity. In particular, we show that estimates of statistical models of the drivers of popularity will be biased if researchers use only trajectories of those practices that once were popular.     

Time Perception and Depressive Realism: Judgment Type,Psychophysical Functions and Bias

The effect of mild depression on time estimation and production was investigated. Participants made both magnitude estimation and magnitude production judgments for five time intervals (specified in seconds) from 3 sec to 65 sec. The parameters of the best fitting psychophysical function (power law exponent, intercept, and threshold) were determined individually for each participant in every condition. There were no significant effects of mood (high BDI, low BDI) or judgment (estimation, production) on the mean exponent, n = .98, 95% confidence interval (.96–1.04) or on the threshold. However, the intercept showed a ‘depressive realism’ effect, where high BDI participants had a smaller deviation from accuracy and a smaller difference between estimation and judgment than low BDI participants. Accuracy bias was assessed using three measures of accuracy: difference, defined as psychological time minus physical time, ratio, defined as psychological time divided by physical time, and a new logarithmic accuracy measure defined as ln (ratio). The ln (ratio) measure was shown to have approximately normal residuals when subjected to a mixed ANOVA with mood as a between groups explanatory factor and judgment and time category as repeated measures explanatory factors. The residuals of the other two accuracy measures flagrantly violated normality. The mixed ANOVAs of accuracy also showed a strong depressive realism effect, just like the intercepts of the psychophysical functions. There was also a strong negative correlation between estimation and production judgments. Taken together these findings support a clock model of time estimation, combined with additional cognitive mechanisms to account for the depressive realism effect. The findings also suggest strong methodological recommendations.     

Mammalian Male Mutation Bias: Impacts of Generation Time and Regional Variation in Substitution Rates

In mammals, males undergo a greater number of germline cell divisions compared with females. Thus, the male germline accumulates more DNA replication errors, which result in male mutation bias—a higher mutation rate for males than for females. The phenomenon of male mutation bias has been investigated mostly for rodents and primates, however, it has not been studied in detail for other mammalian orders. Here we sequenced and analyzed five introns of three genes (DBX/DBY, UTX/UTY, and ZFX/ZFY) homologous between X and Y chromosomes in several species of perissodactyls (horses and rhinos) and of primates. Male mutation bias was evident: substitution rate was higher for a Y chromosome intron than for its X chromosome homologue for all five intron pairs studied. Substitution rates varied regionally among introns sequenced on the same chromosome and this variation influenced male mutation bias inferred from each intron pair. Interestingly, we observed a positive correlation in substitution rates between homologous X and homologous Y introns as well as between orthologous primate and perissodactyl introns. The male-to-female mutation rate ratio estimated from concatenated sequences of five perissodactyl introns was 3.88 (95% CI = 2.90–6.07). Using the data generated here and estimates available in the literature, we compared male mutation bias among several mammalian orders. We conclude that male mutation bias is significantly higher for organisms with long generation times (primates, perissodactyls, and felids) than for organisms with short generation times (e.g., rodents) since the former undergo a greater number of male germline cell divisions. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. [Reviewing Editor: Dr. Deborah Charlesworth]     

Born again at the synapse: a new function for the anaphase promoting complex/cyclosome

Currently, perhaps the most significant biological problem is to understand the mechanisms of learning and memory, and many of the answers will come from molecular explanations of synaptic plasticity. Two new papers have established a surprising connection: the Anaphase Promoting Complex/Cyclosome (APC/C) has a second function in controlling local protein stability at synapses, and hence in the control of behavior .     

The Coalescent Time in the Presence of Background Fertility Selection

Selection ultimately entails differential reproductive success over several generations. This can be measured as a correlation of the number of progeny an individual has with the number of progeny its parent had. This correlation could have a genetic or a cultural basis. The effect of such a correlation is to multiply the single generation sampling variance (Vdeltap) in the diffusion approximation for fixation time by (1-b)+bx(1+r)/(1-r), where bxrn is the correlation between the number of progeny of an individual and its ancestor n generations ago (e.g., b is the heritability and br is the resultant parent-offspring progeny number correlation if the progeny number is genetically determined). This results in a reduction of the fixation or coalescent time by division by this factor. Sex differences in this correlation have been observed, and this provides an explanation for the difference of coalescent times of y-chromosomes and mitochondria.     

Selection Intensity for Codon Bias and the Effective Population Size of Escherichia Coli

The selection intensity for codon bias and the synonymous diversity have been used in the recent literature to estimate the effective population size of Escherichia coli. The results have varied between 10(5) and 10(8). It is suggested here that most of this disparity can be explained by a model that accounts for the population structure of the species. Thus it is assumed that weakly selected characters, like synonymous substitutions, are selectively fixed within individual lines or colonies but spread throughout the population in an essentially neutral way when colonies replace one another. In this way, the effective population size that enters expressions for the codon bias will be that of an individual colony, which, if hitchhiking effects are considered, can be a very small number. The effective population size that appears together with the mutation rate in expressions for the synonymous diversity, on the other hand, will be related to the total number of colonies that make up the species and can be a very large number.     

Mutational Bias and Translational Selection Shaping the Codon Usage Pattern of Tissue-Specific Genes in Rice

The regulatory mechanisms of determining which genes specifically expressed in which tissues are still not fully elucidated, especially in plants. Using internal correspondence analysis, I first establish that tissue-specific genes exhibit significantly different synonymous codon usage in rice, although this effect is weak. The variability of synonymous codon usage between tissues accounts for 5.62% of the total codon usage variability, which has mainly arisen from the neutral evolutionary forces, such as GC content variation among tissues. Moreover, tissue-specific genes are under differential selective constraints, inferring that natural selection also contributes to the codon usage divergence between tissues. These findings may add further evidence in understanding the differentiation and regulation of tissue-specific gene products in plants.     

The Effects of Mutation and Natural Selection on Codon Bias in the Genes of Drosophila

Codon bias varies widely among the loci of Drosophila melanogaster, and some of this diversity has been explained by variation in the strength of natural selection. A study of correlations between intron and coding region base composition shows that variation in mutation pattern also contributes to codon bias variation. This finding is corroborated by an analysis of variance (ANOVA), which shows a tendency for introns from the same gene to be similar in base composition. The strength of base composition correlations between introns and codon third positions is greater for genes with low codon bias than for genes with high codon bias. This pattern can be explained by an overwhelming effect of natural selection, relative to mutation, in highly biased loci. In particular, this correlation is absent when examining fourfold degenerate sites of highly biased genes. In general, it appears that selection acts more strongly in choosing among fourfold degenerate codons than among twofold degenerate codons. Although the results indicate regional variation in mutational bias, no evidence is found for large scale regions of compositional homogeneity.     

Identifying Signatures of Selection in Genetic Time Series

Both genetic drift and natural selection cause the frequencies of alleles in a population to vary over time. Discriminating between these two evolutionary forces, based on a time series of samples from a population, remains an outstanding problem with increasing relevance to modern data sets. Even in the idealized situation when the sampled locus is independent of all other loci, this problem is difficult to solve, especially when the size of the population from which the samples are drawn is unknown. A standard χ²-based likelihood-ratio test was previously proposed to address this problem. Here we show that the χ²-test of selection substantially underestimates the probability of type I error, leading to more false positives than indicated by its P-value, especially at stringent P-values. We introduce two methods to correct this bias. The empirical likelihood-ratio test (ELRT) rejects neutrality when the likelihood-ratio statistic falls in the tail of the empirical distribution obtained under the most likely neutral population size. The frequency increment test (FIT) rejects neutrality if the distribution of normalized allele-frequency increments exhibits a mean that deviates significantly from zero. We characterize the statistical power of these two tests for selection, and we apply them to three experimental data sets. We demonstrate that both ELRT and FIT have power to detect selection in practical parameter regimes, such as those encountered in microbial evolution experiments. Our analysis applies to a single diallelic locus, assumed independent of all other loci, which is most relevant to full-genome selection scans in sexual organisms, and also to evolution experiments in asexual organisms as long as clonal interference is weak. Different techniques will be required to detect selection in time series of cosegregating linked loci.     

Tracking Visual Events in Time in the Absence of Time Perception: Implicit Processing at the ms Level

Previous studies have suggested that even if subjects deem two visual stimuli less than 20 ms apart to be simultaneous, implicitly they are nonetheless distinguished in time. It is unclear, however, how information is encoded within this short timescale. We used a priming paradigm to demonstrate how successive visual stimuli are processed over time intervals of less than 20 ms. The primers were two empty square frames displayed either simultaneously or with a 17ms asynchrony. The primers were followed by the target information after a delay of 25 ms to 100 ms. The two square frames were filled in one after another with a delay of 100 ms between them, and subjects had to decide on the location of the first of the frames to be filled in. In a second version of the paradigm, only one square frame was filled in, and subjects had to decide where it was positioned. The influence of the primers is revealed through faster response times depending on the location of the first and second primers. Experiment 1 replicates earlier results, with a bias towards the side of the second primer, but only when there is a delay of 75 to 100 ms between primers and targets. The following experiments suggest this effect to be relatively independent of the task context, except for a slight effect on the time course of the biases. For the temporal order judgment task, identical results were observed when subjects have to answer to the side of the second rather than the first target, showing the effect to be independent of the hand response, and suggesting it might be related to a displacement of attention. All in all the results suggest the flow of events is followed more efficiently than suggested by explicit asynchrony judgment studies. We discuss the possible impact of these results on our understanding of the sense of time continuity.