Theoretical Study of near Neutrality. I. Heterozygosity and Rate of Mutant Substitution

In order to clarify the nature of "near neutrality" in molecular evolution and polymorphism, extensive simulation studies were performed. Selection coefficients of new mutations are assumed to be small so that both random genetic drift and selection contribute to determining the behavior of mutants. The model also incorporates normally distributed spatial fluctuation of selection coefficients. If the system starts from "average neutrality," it will move to a better adapted state, and most new mutations will become "slightly deleterious." Monte Carlo simulations have indicated that such adaptation is attained, but that the rate of such "progress" is very low for weak selection. In general, the larger the population size, the more effective the selection becomes. Also, as selection becomes weaker, the behavior of the mutants approaches that of completely neutral genes. Thus, the weaker the selection, the smaller is the effect of population size on mutant dynamics. Increase of heterozygosity with population size is very pronounced for subdivided populations. The significance of these results is discussed in relation to various observed facts on molecular evolution and polymorphism, such as generation-time dependency and overdispersion of the molecular clock, or contrasting patterns of DNA and protein polymorphism among some closely related species.     

Effective Size and F-Statistics of Subdivided Populations. II. Dioecious Species

Assuming discrete generations and autosomal inheritance involving genes that do not affect viability or reproductive ability, we have derived recurrence equations for the inbreeding coefficient and coancestry between individuals within and among subpopulations for a subdivided monoecious population with arbitrary distributions of male and female gametes per family, variable pollen and seed migration rates, and partial selfing. From the equations, formulas for effective size and expressions for F-statistics are obtained. For the special case of a single unsubdivided population, our equations reduce to the simple expressions derived by previous authors. It is shown that population structure (subdivision and migration) is important in determining the inbreeding coefficient and effective size. Failure to recognize internal structures of populations may lead to considerable bias in predicting effective size. Inbreeding coefficient, coancestry between individuals within and among subpopulations accrue at different and variable rates over initial generations before they converge to the same asymptotic rate of increase. For a given population, the smaller the pollen and seed migration rates, the more generations are required to attain the asymptotic rate and the larger the asymptotic effective size. The equations presented herein can be used for the study of evolutionary biology and conservation genetics.     

Population Genetics of Euphydryas Butterflies. I. Genetic Variation and the Neutrality Hypothesis

Twenty-one populations of the checkerspot butterfly, Euphydryas editha, and ten populations of Euphydryas chalcedona were sampled for genetic variation at eight polymorphic enzyme loci. Both species possessed loci that were highly variable from population to population and loci that were virtually identical across all populations sampled. Our data indicate that the neutrality hypothesis is untenable for the loci studied, and therefore selection is indicated as the major factor responsible for producing these patterns. Thorough ecological work allowed gene flow to be ruled out (in almost all instances) as a factor maintaining similar gene frequencies across populations. The Lewontin-Krakauer test indicated magnitudes of heterogeneity among standardized variances of gene frequencies inconsistent with the neutrality hypothesis. The question of whether or not to correct this statistic for sample size is discussed. Observed equitability of gene frequencies of multiple allelic loci was found to be greater than that predicted under the neutrality hypothesis. Genetic differentiation persisting through two generations was found between the one pair of populations known to exchange significant numbers of individuals per generation. Two matrices of genetic distance between populations, based on the eight loci sampled, were found to be significantly correlated with a matrix of environmental distance, based on measures of fourteen environmental parameters. Correlations between gene frequencies and environmental parameters, results of multiple regression analysis, and results of principle component analysis showed strong patterns of association and of "explained" variation. The correlation analyses suggest which factors might be further investigated as proximate selective agents.     

Decline and Local Extinction of Caribbean Eusocial Shrimp

The tropical shrimp genus Synalpheus includes the only eusocial marine animals. In much of the Caribbean, eusocial species have dominated the diverse fauna of sponge-dwelling shrimp in coral rubble for at least the past two decades. Here we document a recent, dramatic decline and apparent local extinction of eusocial shrimp species on the Belize Barrier Reef. Our collections from shallow reefs in central Belize in 2012 failed to locate three of the four eusocial species formerly abundant in the area, and showed steep declines in colony size and increases in frequency of queenless colonies prior to their disappearance. Concordant with these declines, several nonsocial, pair-forming Synalpheus species increased in frequency. The decline in eusocial shrimp is explained in part by disappearance of two sponge species on which they specialize. Eusocial shrimp collections from Jamaica in 2012 showed similar patterns of decline in colony size and increased queenlessness compared with prior Jamaican collections. The decline and local extinction of eusocial shrimp happened against a backdrop of changes in coral assemblages during recent decades, and may reflect changes in abundance and quality of dead coral substratum and succession of the diverse cryptic organisms living within it. These changes document potentially worrisome declines in a unique taxon of eusocial marine animals.     

Relationship between the Reproductive Capability of a Subdivided Population and Its Genetic Structure: Computer Simulation,the Simplest Multilocus Model

An original computer model, simulating joint genetic and demographic dynamics of subdivided populations, is proposed. The model accounts for the reverse effect of the genetic structure on the reproductive capability of a population, which is based on a postulated limited set of biallelic loci, controlling variation in an adaptive quantitative trait. The model allows to simulate spreading of the population, originating from a single small colony, resulting in establishment of involves genetic and demographic equilibrium (a normal population-genetic process) and reorganization of the genetic structure of the subdivided population under anthropogenic pressure, associated with a decrease in its reproductive capability (an adverse population-genetic process).     

Minimum Area Required for Local Populations of Japanese Macaques Estimated from the Relationship Between Habitat Area and Population Extinction

I estimated the minimum area required (MAR) for local populations of Japanese macaques (Macaca fuscata) from empirical data on habitat area and population viability. I used logistic regression analysis to examine the relationship between habitat area and survival/extinction among 50 populations over 50 yr. Estimated habitat areas satisfying 95–99% probability of a population persisting for 100–1000 yr range from 525 to 975 km². However, confidence limits of parameters in the logistic regression equation are very large. Moreover, the number of extinct population might be underestimated in the empirical data. Consequently, a much wider habitat area (>1000 km²) should be considered for actual conservation planning for local populations of Japanese macaques. The method involves fewer variables and assumptions than previous methods of MAR estimation, and therefore may be a more useful way to estimate MAR for various species and regions.     

一个受到开发的扩散种群的持续与灭绝条件

本文讨论了一类服从Logistic增长规律,且受到开发或捕食的扩散种群动态模型：的渐近性质,得到了具有重要生物学意义的结果．即当开发或捕食率时,种群稳定持续,否则趋于灭绝．     

具有毒素影响的二维Kolmogorov模型的持续生存和绝灭

研究了具有毒素影响的二维Kolmogorov模型,给出了该系统持续生存与绝灭的充分条件.     

Small but Powerful: Top Predator Local Extinction Affects Ecosystem Structure and Function in an Intermittent Stream

Top predator loss is a major global problem, with a current trend in biodiversity loss towards high trophic levels that modifies most ecosystems worldwide. Most research in this area is focused on large-bodied predators, despite the high extinction risk of small-bodied freshwater fish that often act as apex consumers. Consequently, it remains unknown if intermittent streams are affected by the consequences of top-predators’ extirpations. The aim of our research was to determine how this global problem affects intermittent streams and, in particular, if the loss of a small-bodied top predator (1) leads to a ‘mesopredator release’, affects primary consumers and changes whole community structures, and (2) triggers a cascade effect modifying the ecosystem function. To address these questions, we studied the top-down effects of a small endangered fish species, Barbus meridionalis (the Mediterranean barbel), conducting an enclosure/exclosure mesocosm experiment in an intermittent stream where B. meridionalis became locally extinct following a wildfire. We found that top predator absence led to ‘mesopredator release’, and also to ‘prey release’ despite intraguild predation, which contrasts with traditional food web theory. In addition, B. meridionalis extirpation changed whole macroinvertebrate community composition and increased total macroinvertebrate density. Regarding ecosystem function, periphyton primary production decreased in apex consumer absence. In this study, the apex consumer was functionally irreplaceable; its local extinction led to the loss of an important functional role that resulted in major changes to the ecosystem’s structure and function. This study evidences that intermittent streams can be affected by the consequences of apex consumers’ extinctions, and that the loss of small-bodied top predators can lead to large ecosystem changes. We recommend the reintroduction of small-bodied apex consumers to systems where they have been extirpated, to restore ecosystem structure and function.     

The Genetic Structure of Natural Populations of DROSOPHILA MELANOGASTER. Xi. Genetic Variability in a Local Population

Six hundred and ninety-one second chromosomes were extracted from a Raleigh, North Carolina population, and the following experimental results were obtained: (1) Salivary gland chromosomes of all lines were observed and the number of inversion-carrying chromosomes was 130, among which 76 carried In(2R)NS, 36 carried In(2L)t, 4 carried In(2L)t and In(2R)NS, and 14 carried different kinds of rare inversions. (2) Viabilities of homozygotes and heterozygotes were examined. The frequency of lethal-carrying chromosomes was 275/691 (or 0.398):70/130 (or 0.538) in inversion-carrying chromosomes and 205/561 (or 0.365) in inversion-free chromosomes. The former is significantly higher than the latter. The average homozygote viability was 0.4342 including lethal lines and 0.7163 excluding those, the average heterozygote viability being 1.0000. The detrimental load to lethal load ratio (D:L ratio) was 0.334/0.501 = 0.67. The average viability of lethal heterozygotes was less than that of lethal-free heterozygotes, significantly in inversion-free individuals but not significantly so in inversion-carrying individuals. Inversion heterozygotes seem to have slightly better viability than the inversion-free heterozygotes on the average, but not significantly so. (3) The average degree of dominance of viability polygenes was estimated to be 0.293 +/- 0.071 for all heterozygotes whose component chromosomes had better viabilities than 0.6 of the average heterozygote viability, 0.177 +/- 0.077 for inversion-free heterozygotes and 0.489 +/- 0.082 for inversion heterozygotes. (4) Mutation rates of viability polygenes and lethal genes were estimated on the basis of genetic loads and average degrees of dominance of lethal genes and viability polygenes. Estimates were very close to those obtained by direct estimation. (5) Possible overdominance and epistasis were detected, but the magnitude must be very small. (6) The effective size of the population was estimated to be much greater than 10,000 by using the allelism rate of lethal-carrying chromosomes (0.0040) and their frequency.-On the basis of these findings and the comparison with the predicted result (Mukai and Maruyama 1971), the mechanisms of the maintenance of genetic variability in the population are discussed.     

A Population Study of Spergula arvensis: II. Genetics and Breeding Behaviour

Genetical work has shown that one locus and two alleles areinvolved in the inheritance of a seedcoat character of Spergulaarvensis and that the hetero-zygous plants are intermediatebetween the two types of homozygous plants. It is thus possibleto translate directly the observed geographical distributionof plant frequencies into distribution of gene frequencies.The inheritance of a hairiness character is less straightforward,but in this case also there is a close relation between phenotypefrequency distribution and gene frequency distribution. Using the fact that plants heterozygous for the seedcoat characterare easily recognizable, it is possible to make an estimateof the amount of outbreeding in nature. The values obtainedrange from o to 3 per cent.     

The Allelic Correlation Structure of Gainj- and Kalam-Speaking People. II. the Genetic Distance between Population Subdivisions

The patterning of allele frequency variability among 18 local groups of Gainj and Kalam speakers of highland Papua New Guinea is investigated using new genetic distance methods. The genetic distances proposed here are obtained by decomposing Sewall Wright's coefficient FST into a set of coefficients corresponding to all pairs of population subdivisions. Two statistical methods are given to estimate these quantities. One method provides estimates weighted by sample sizes, while the other method does not use sample size weighting. Both methods correct for the within-individual and between-individual-within-groups sums of squares. Genetic distances among the Gainj and Kalam subdivisions are analyzed with respect to demographic, geographic, and linguistic variables. We find that a demographic feature, group size, has the greatest demonstrable association with the patterning of genetic distances. The pattern of geographic distances among groups displays a weak congruence with the pattern of genetic distances, and the association of genetic and linguistic diversity is very low. An effect of differences in group size on genetic distances is not surprising, from basic theoretical considerations, but genetic distances have not often been analyzed with respect to these variables in the past. The lack of correspondence between genetic distances and linguistic and geographic differences is an unusual feature that distinguishes the Gainj and Kalam from most other tribal populations.     

Theoretical chemistry of gold. II

Gold is an element whose unique properties are strongly influenced by relativistic effects. A large body of appropriate calculations now exists and its main conclusions are summarized. The present paper completes the recent review by [Pyykkö, Angew. Chem. Int. Ed., 43 (2004) 4412].     

用SSR研究栲树群体遗传结构

利用微卫星（SSR）分子标记对福建省内4个栲树（Castanopsis fargesii Franch.）群体遗传结构进行了研究。SSR标记揭示了栲树群体丰富的遗传变异：平均等位基因数A＝9.0,平均有效等位基因数Ne=4.8,平均期望杂合度He=0.65,而群 具有较低的Fst值（Fst=0.031）。SSR每个位点的等位基因频率分布在栲树群体间都存在显或极显差异,表明根据SSR等位基因频率分布亦能了解各体的分化。SSR标记使栲树群体中一些稀有等位基因得以表现,54个SSR等位基因中有15个等位基因仅出现在1个或2个群体中,且频率较低,在遗传多样性保护中更应注重保护这些稀有的等位变异。     

The Genetic Structure of Natural Populations of DROSOPHILA MELANOGASTER. Xii. Linkage Disequilibrium in a Large Local Population

Seven hundred and three second chromosomes were extracted from a Raleigh, North Carolina population of Drosophila melanogaster in 1970. Additionally, four hundred and eighty-nine third chromosomes were extracted from a large cage population founded from the flies in the 1970 Raleigh collection. The alpha glycerol-3-phosphate dehydrogenase-1, malate dehydrogenase-1, alcohol dehydrogenase, and alpha amylase loci were studied from the second chromosomes, and the esterase-6, esterase-C, and octanol dehydrogenase loci were analyzed from the third chromosomes. Inversions, relative viability and fecundity were studied for both classes of chromosomes. The following significant findings were obtained: (1) All loci examined were polymorphic or had at least two alleles at appreciable frequencies. Analysis of the combined data from this experiment with that of Mukai, Mettler and Chigusa (1971) revealed that the frequencies of the genes in the second chromosomes collected in early August were approximately the same over three years. (2) Linkage disequilibria between and among isozyme genes inter se were not detected except in a few cases which can be considered due to non-random sampling. (3) Linkage disequilibria between isozyme genes and polymorphic inversions were detected when the recombination values between the breakage points of the inversions and the genes in question were small. In only a few cases, were second and third order linkage disequilibria including polymorphic inversions detected. (4) Evidence for either variation among genotypes within loci or cumulative effects of heterozygosity was found for viability and fecundity. As a result of these findings, it was tentatively concluded that although selection might be perceptibly operating on some polymorphic isozyme loci, most of the polymorphic isozyme genes are selectively neutral or near-neutral in the populations studied.     

Effective Size and F-Statistics of Subdivided Populations. I. Monoecious Species with Partial Selfing

Assuming discrete generations and autosomal inheritance involving genes that do not affect viability or reproductive ability, we have derived recurrence equations for the inbreeding coefficient and coancestry between individuals within and among subpopulations for a subdivided monoecious population with arbitrary distributions of male and female gametes per family, variable pollen and seed migration rates, and partial selfing. From the equations, formulas for effective size and expressions for F-statistics are obtained. For the special case of a single unsubdivided population, our equations reduce to the simple expressions derived by previous authors. It is shown that population structure (subdivision and migration) is important in determining the inbreeding coefficient and effective size. Failure to recognize internal structures of populations may lead to considerable bias in predicting effective size. Inbreeding coefficient, coancestry between individuals within and among subpopulations accrue at different and variable rates over initial generations before they converge to the same asymptotic rate of increase. For a given population, the smaller the pollen and seed migration rates, the more generations are required to attain the asymptotic rate and the larger the asymptotic effective size. The equations presented herein can be used for the study of evolutionary biology and conservation genetics.     

The Role of Density Regulation in Extinction Processes and Population Viability Analysis

We review the role of density dependence in the stochastic extinction of populations and the role density dependence has played in population viability analysis (PVA) case studies. In total, 32 approaches have been used to model density regulation in theoretical or applied extinction models, 29 of them are mathematical functions of density dependence, and one approach uses empirical relationships between density and survival, reproduction, or growth rates. In addition, quasi-extinction levels are sometimes applied as a substitute for density dependence at low population size. Density dependence further has been modelled via explicit individual spacing behaviour and/or dispersal. We briefly summarise the features of density dependence available in standard PVA software, provide summary statistics about the use of density dependence in PVA case studies, and discuss the effects of density dependence on extinction probability. The introduction of an upper limit for population size has the effect that the probability of ultimate extinction becomes 1. Mean time to extinction increases with carrying capacity if populations start at high density, but carrying capacity often does not have any effect if populations start at low numbers. In contrast, the Allee effect is usually strong when populations start at low densities but has only a limited influence on persistence when populations start at high numbers. Contrary to previous opinions, other forms of density dependence may lead to increased or decreased persistence, depending on the type and strength of density dependence, the degree of environmental variability, and the growth rate. Furthermore, effects may be reversed for different quasi-extinction levels, making the use of arbitrary quasi-extinction levels problematic. Few systematic comparisons of the effects on persistence between different models of density dependence are available. These effects can be strikingly different among models. Our understanding of the effects of density dependence on extinction of metapopulations is rudimentary, but even opposite effects of density dependence can occur when metapopulations and single populations are contrasted. We argue that spatially explicit models hold particular promise for analysing the effects of density dependence on population viability provided a good knowledge of the biology of the species under consideration exists. Since the results of PVAs may critically depend on the way density dependence is modelled, combined efforts to advance statistical methods, field sampling, and modelling are urgently needed to elucidate the relationships between density, vital rates, and extinction probability.