J. B. S. Haldane, Ernst Mayr and the Beanbag Genetics Dispute

Starting from the early decades of the twentieth century, evolutionary biology began to acquire mathematical overtones. This took place via the development of a set of models in which the Darwinian picture of evolution was shown to be consistent with the laws of heredity discovered by Mendel. The models, which came to be elaborated over the years, define a field of study known as population genetics. Population genetics is generally looked upon as an essential component of modern evolutionary theory. This article deals with a famous dispute between J. B. S. Haldane, one of the founders of population genetics, and Ernst Mayr, a major contributor to the way we understand evolution. The philosophical undercurrents of the dispute remain relevant today. Mayr and Haldane agreed that genetics provided a broad explanatory framework for explaining how evolution took place but differed over the relevance of the mathematical models that sought to underpin that framework. The dispute began with a fundamental issue raised by Mayr in 1959: in terms of understanding evolution, did population genetics contribute anything beyond the obvious? Haldane's response came just before his death in 1964. It contained a spirited defense, not just of population genetics, but also of the motivations that lie behind mathematical modelling in biology. While the difference of opinion persisted and was not glossed over, the two continued to maintain cordial personal relations.

Does size matter? Comparative population genetics of two butterflies with different wingspans

The dispersal ability of a species is central to its biology, affecting other processes like local adaptation, population and community dynamics, and genetic structure. Among the intrinsic, species-specific factors that affect dispersal ability in butterflies, wingspan was recently shown to explain a high amount of variance in dispersal ability. In this study, a comparative approach was adopted to test whether a difference in wingspan translates into a difference in population genetic structure. Two closely related butterfly species from subfamily Satyrinae, family Nymphalidae, which are similar with respect to all traits that affect dispersal ability except for wingspan, were studied. Melanitis leda (wingspan 60-80 mm) and Ypthima baldus (wingspan 30-40 mm) were collected from the same areas along the Western Ghats of southern India. Amplified fragment length polymorphisms were used to test whether the species with a higher wingspan (M. leda) exhibited a more homogenous population genetic structure, as compared to a species with a shorter wingspan (Y. baldus). In all analyses, Y. baldus exhibited greater degree of population genetic structuring. This study is one of the few adopting a comparative approach to establish the relationship between traits that affect dispersal ability and population genetic structure.

Transposable element-medicated evolution of sex: A population genetic model

For a population made up of individuals capable of sexual as well as asexual modes of reproduction, conditions for the spread of a transposable element are explored using a one-locus, two-haplotype model. The analysis is then extended to include the possibility that the transposable element can modulate the probability of sexual reproduction, thus casting Hickey’s (1982,Genetics 101: 519–531) suggestion in a population genetics framework. The model explicitly includes the cost of sexual reproduction, fitness disadvantage to the transposable element, probability of transposition, and the predisposition for sexual reproduction in the presence and absence of the transposable element. The model predicts several kinds of outcome, including initial frequency dependence and stable polymorphism. More importantly, it is seen that for a wide range of parameter values, the transposable element can go to fixation. Therefore it is able to convert the population from a predominantly asexual to a predominantly sexual mode of reproduction. Viewed in conjunction with recent results implicating short stretches of apparently non-coding DNA in sex determination (McCoubreyet al. 1988,Science 242: 1146–1151), the model hints at the important role this mechanism could have played in the evolution of sexuality.

Sibling rivalry between seeds within a fruit: Some population genetic models

Competition between seeds within a fruit for parental resources is described using one-locus-two-allele models. While a �normal� allele leads to an equitable distribution of resources between seeds (a situation which also corresponds to the parental optimum), the �selfish� allele is assumed to cause the seed carrying it to usurp a higher proportion of the resources. The outcome of competition between �selfish� alleles is also assumed to lead to an asymmetric distribution of resources, the �winner� being chosen randomly. Conditions for the spread of an initially rare selfish allele and the optimal resource allocation corresponding to the evolutionarily stable strategy, derived for species with n-seeded fruits, are in accordance with expectations based on Hamilton�s inclusive fitness criteria. Competition between seeds is seen to be most intense when there are only two seeds, and decreases with increasing number of seeds, suggesting that two-seeded fruits would be rarer than one-seeded or many-seeded ones. Available data from a large number of plant species are consistent with this prediction of the model.

A finite population model of mlecular evolution: theory and computation

This article is concerned with the evolution of haploid organisms that reproduce asexually. In a seminal piece of work, Eigen and coauthors proposed the quasispecies model in an attempt to understand such an evolutionary process. Their work has impacted antiviral treatment and vaccine design strategies. Yet, predictions of the quasispecies model are at best viewed as a guideline, primarily because it assumes an infinite population size, whereas realistic population sizes can be quite small. In this paper we consider a population genetics-based model aimed at understanding the evolution of such organisms with finite population sizes and present a rigorous study of the convergence and computational issues that arise therein. Our first result is structural and shows that, at any time during the evolution, as the population size tends to infinity, the distribution of genomes predicted by our model converges to that predicted by the quasispecies model. This justifies the continued use of the quasispecies model to derive guidelines for intervention. While the stationary state in the quasispecies model is readily obtained, due to the explosion of the state space in our model, exact computations are prohibitive. Our second set of results are computational in nature and address this issue. We derive conditions on the parameters of evolution under which our stochastic model mixes rapidly. Further, for a class of widely used fitness landscapes we give a fast deterministic algorithm which computes the stationary distribution of our model. These computational tools are expected to serve as a framework for the modeling of strategies for the deployment of mutagenic drugs.

Variation at Diabetes- and Obesity-Associated Loci May Mirror Neutral Patterns of Human Population Diversity and Diabetes Prevalence in India

South Asian populations harbor a high degree of genetic diversity, due in part to demographic history. Two studies on genome-wide variation in Indian populations have shown that most Indian populations show varying degrees of admixture between ancestral north Indian and ancestral south Indian components. As a result of this structure, genetic variation in India appears to follow a geographic cline. Similarly, Indian populations seem to show detectable differences in diabetes and obesity prevalence between different geographic regions of the country. We tested the hypothesis that genetic variation at diabetes-and obesity-associated loci may be potentially related to different genetic ancestries. We genotyped 2977 individuals from 61 populations across India for 18 SNPs in genes implicated in T2D and obesity. We examined patterns of variation in allele frequency across different geographical gradients and considered state of origin and language affiliation. Our results show that most of the 18 SNPs show no significant correlation with latitude, the geographic cline reported in previous studies, or by language family. Exceptions include KCNQ1 with latitude and THADA and JAK1 with language, which suggests that genetic variation at previously ascertained diabetes-associated loci may only partly mirror geographic patterns of genome-wide diversity in Indian populations.

Effects of Withania somnifera and Tinospora cordifolia Extracts on the Side Population Phenotype of Human Epithelial Cancer Cells: Toward Targeting Multidrug Resistance in Cancer

Recent reports suggest the existence of a subpopulation of stem-like cancer cells, termed as cancer stem cells (CSCs), which bear functional and phenotypic resemblance with the adult, tissue-resident stem cells. Side population (SP) assay based on differential efflux of Hoechst 33342 has been effectively used for the isolation of CSCs. The drug resistance properties of SP cells are typically due to the increased expression of ABC transporters leading to drug efflux. Conventionally used chemotherapeutic drugs may often leads to an enrichment of SP, revealing their inability to target the drug-resistant SP and CSCs. Thus, identification of agents that can reduce the SP phenotype is currently in vogue in cancer therapeutics. Withania somnifera (WS) and Tinospora cordifolia (TC) have been used in Ayurveda for treating various diseases, including cancer. In the current study, we have investigated the effects of ethanolic (ET) extracts of WS and TC on the cancer SP phenotype. Interestingly, we found significant decrease in SP on treatment with TC-ET, but not with WS-ET. The SP-inhibitory TC-ET was further fractionated into petroleum ether (TC-PET), dichloromethane (TC-DCM), and n-butyl alcohol (TC-nBT) fractions using bioactivity-guided fractionation. Our data revealed that TC-PET and TC-DCM, but not TC-nBT, significantly inhibited SP in a dose-dependent manner. Furthermore, flow cytometry-based functional assays revealed that TC-PET and TC-DCM significantly inhibited ABC-B1 and ABC-G2 transporters and sensitized cancer cells toward chemotherapeutic drug-mediated cytotoxicity. Thus, the TC-PET and TC-DCM may harbor phytochemicals with the potential to reverse the drug-resistant phenotype, thus improving the efficacy of cancer chemotherapy.

DNA methylation analysis of phenotype specific stratified Indian population

Background: DNA methylation and its perturbations are an established attribute to a wide spectrum of phenotypic variations and disease conditions. Indian traditional system practices personalized medicine through indigenous concept of distinctly descriptive physiological, psychological and anatomical features known as prakriti. Here we attempted to establish DNA methylation differences in these three prakriti phenotypes. Methods: Following structured and objective measurement of 3416 subjects, whole blood DNA of 147 healthy male individuals belonging to defined prakriti (Vata, Pitta and Kapha) between the age group of 20-30years were subjected to methylated DNA immunoprecipitation (MeDIP) and microarray analysis. After data analysis, prakriti specific signatures were validated through bisulfite DNA sequencing. Results: Differentially methylated regions in CpG islands and shores were significantly enriched in promoters/UTRs and gene body regions. Phenotypes characterized by higher metabolism (Pitta prakriti) in individuals showed distinct promoter (34) and gene body methylation (204), followed by Vata prakriti which correlates to motion showed DNA methylation in 52 promoters and 139 CpG islands and finally individuals with structural attributes (Kapha prakriti) with 23 and 19 promoters and CpG islands respectively. Bisulfite DNA sequencing of prakriti specific multiple CpG sites in promoters and 5'-UTR such as; LHX1 (Vata prakriti), SOX11 (Pitta prakriti) and CDH22 (Kapha prakriti) were validated. Kapha prakriti specific CDH22 5'-UTR CpG methylation was also found to be associated with higher body mass index (BMI). Conclusion: Differential DNA methylation signatures in three distinct prakriti phenotypes demonstrate the epigenetic basis of Indian traditional human classification which may have relevance to personalized medicine.

A Ca2+-dependent early functional heterogeneity in amoebae of Dictyostelium discoideum, revealed by flow cytometry

When freshly starved amoebae of Dictyostelium discoideum are loaded with the Ca2+-specific dye indo-1/AM and analyzed in a fluorescence-activated cell sorter, they exhibit a quasi-bimodal distribution of fluorescence. This permits a separation of the population into two classes: H, or ''high Ca2+-indo-1 fluorescence,'' and L, or ''low Ca2+-indo-1 fluorescence.'' Simultaneous monitoring of Ca2+-indo-1 and Ca2+-chlortetracycline fluorescence shows that by and large the same cells tend to have high (or low) levels of both cytoplasmic and sequestered Ca2+. Next we label H cells with tetramethylrhodamine isothiocyanate (TRITC) and mix them in a 1:4 ratio with L cells, In the slugs that result, TRITC fluorescence is confined mainly to the anterior prestalk region. This implies that amoebae with relatively high Ca2+ at the vegetative stage tend to develop into prestalk cells and those with low Ca2+ into prespores. Polysphondylium violaceum, a cellular slime mold that does not possess prestalk and prespore cells, also does not display a Ca2+-dependent heterogeneity at the vegetative stage or in slugs. Finally, confirming earlier findings with the fluorophore fura-2 (Azhar ef al., Curr. Sci. 68, 337-342 (1995)), a prestalk-prespore difference in cellular Ca2+ is present in the cells of the slug in vivo. These findings are discussed in light of the possible roles of Ca2+ for cell differentiation in D. discoideum.

Precipitation In Small Systems--II. Mean Field Equations More Effective Than Population Balance

Part I (Manjunath et al., 1994, Chem. Engng Sci. 49, 1451-1463) of this paper showed that the random particle numbers and size distributions in precipitation processes in very small drops obtained by stochastic simulation techniques deviate substantially from the predictions of conventional population balance. The foregoing problem is considered in this paper in terms of a mean field approximation obtained by applying a first-order closure to an unclosed set of mean field equations presented in Part I. The mean field approximation consists of two mutually coupled partial differential equations featuring (i) the probability distribution for residual supersaturation and (ii) the mean number density of particles for each size and supersaturation from which all average properties and fluctuations can be calculated. The mean field equations have been solved by finite difference methods for (i) crystallization and (ii) precipitation of a metal hydroxide both occurring in a single drop of specified initial supersaturation. The results for the average number of particles, average residual supersaturation, the average size distribution, and fluctuations about the average values have been compared with those obtained by stochastic simulation techniques and by population balance. This comparison shows that the mean field predictions are substantially superior to those of population balance as judged by the close proximity of results from the former to those from stochastic simulations. The agreement is excellent for broad initial supersaturations at short times but deteriorates progressively at larger times. For steep initial supersaturation distributions, predictions of the mean field theory are not satisfactory thus calling for higher-order approximations. The merit of the mean field approximation over stochastic simulation lies in its potential to reduce expensive computation times involved in simulation. More effective computational techniques could not only enhance this advantage of the mean field approximation but also make it possible to use higher-order approximations eliminating the constraints under which the stochastic dynamics of the process can be predicted accurately.

On the moulding of population viscosity by natural selection

In this paper, we explore the conjoint evolution of dispersal and social behaviour. The model investigated is of a population distributed over a number of sites each with a carrying capacity of two adults and an episode of dispersal in the juvenile stage. The fertilities are governed by whether an individual and its neighbour are selfish or co-operative. It is shown that the best dispersal strategy for the co-operative genotype always involves lower levels of dispersal; and further that ecological conditions favouring low levels of dispersal increase the selective advantage of a co-operative genotype. Given this positive feedback, we suggest that in any taxon viscosity and co-operativity will tend to be correlated and bimodally distributed. Hence we predict the existence of two kinds of animal societies; viscous and co-operative (e.g. quasi-social wasps such as Mischocyttarus), and non-viscous and selfish (e.g. communal sphecid wasps such as Cerceris), and relatively few social groups with intermediate levels of co-operativity and viscosity. We also suggest that when one of the two sexes disperses, it will be the sex with lower potential for co-operative behaviour.

Inbreeding and the incidence of childhood genetic disorders in Karnataka, South India

Consanguineous marriages are strongly favoured among the populations of South India. In a study conducted on 407 infants and children, a total of 35 genetic diseases was diagnosed in 63 persons: 44 with single gene defects, 12 with polygenic disorders, and seven with Down's syndrome. The coefficient of inbreeding of the total study group, F = 0.0414, was significantly higher than that previously calculated for the general population, F = 0.0271, and autosomal recessive disorders formed the largest single disease category diagnosed. The results suggest that long term inbreeding may not have resulted in appreciable elimination of recessive lethals and sub-lethals from the gene pool.

Evolution of cooperation by reciprocation within structured demes

The iterative two-person Prisoners’ Dilemma game has been generalised to theN-person case. The evolution of cooperation is explored by matching the Tit For Tat (TFT) strategy (Axelrod and Hamilton 1981) against the selfish strategy. Extension of TFT toN-person situations yields a graded set of strategies from the softest TFT, which continues cooperation even if only one of the opponents reciprocates it, to the hardest, which would do so only when all the remaining opponents cooperate. The hardest TFT can go to fixation against the selfish strategy provided it crosses a threshold frequencypc. All the other TFT are invadable by the selfish (D) or the pure defector strategy, while none can invadeD. Yet, provided a thresholdpc is crossed, they can coexist stably withD. AsN, the size of the group increases, the threshold pc also increases, indicating that the evolution of cooperation is more difficult for larger groups. Under certain conditions, only the soft TFT can coexist stably against the selfish strategyD, while the harder ones cannot. An interesting possibility of a complete takeover of the selfish population by successive invasions by harder and harder TFT strategies is also presented.