Some further results from a stochastic fertility model for women

Any realistic model of human fertility should encompass the distributions and interactions of three time intervals a fecund married woman may experience repeatedly in her childbearing period: (1) waiting time for a conception, (2) gestation period, and (3) period of postpartum amenorrhea. Perrin & Sheps (1964) presented a model in which human reproduction is viewed as a Markov renewal process with a finite number of states. Das Gupta (1973b) presented a general probability model of fertility along the lines suggested by Perrin & Sheps which removes two limitations of their model. First, it does not assume that the distributions of durations of stay in the fertility states are independent of each other. Second, it allows us to study the effect of breast-feeding on demographic characteristics, such as interval between live births or birth rate. Results derived in Das Gupta (1973b) include the distributions of time intervals and the exact probabilities of different states at a particular time. The present paper includes additional results pertaining to the same general model, such as the distribution of number of conceptions in a fixed period of time, the distribution of time needed for a fixed number of conceptions, pregnancy rate and fertility rate, and the distribution of the time elapsed since last live birth. The general results are applied to specific models to obtain some known results.     

The role of migration in the genetic structure of populations in temporally and spatially varying environments II. Island Models

The Island Model introduced by Sewall Wright (1951) has proven to be a useful construction for studying the interaction of genetic drift, population subdivision, and mutation. Interest in the model has recently increased because of its relevance to certain questions involving the rate of differentiation of sub-populations under the neutral allele hypothesis (e.g., Smith, 1970; Latter, 1973). It is perhaps the only realistic population structure in which the test for neutrality proposed by Lewontin and Krakauer (1973) is valid (Lewontin and Krakauer, 1975). If data from natural populations is to be compared to the predictions of the Island Model, it is desirable to have an alternative model with the same migration pattern but with natural selection operating. In this paper one such model will be introduced where the stochastic element comes from random fluctuations in the environment rather than from genetic drift. The model is a direct extension of the one in the previous paper in this series (Gillespie, 1975) which dealt with a population which is subdivided into two patches with restricted migration between them.     

An approach to a physico-chemical model of olfactory stimulation in vertebrates by single compounds

During recent years, numerous attempts have been made to correlate both quantitative (Davies &; Taylor, 1959; Engen, 1962; Beck, 1964; Engen, Cain &; Rovee, 1968; Cain, 1969; Dravnieks &; Laffoit, 1970; Laffort, 1969a,b) and qualitative (Davies, 1965; Amoore &; Venstrom, 1965; Döving, 1966a,b; Wright &; Michels, 1964; Leveteau &; MacLeod, 1969) odorous properties of single compounds to their molecular properties. These attempts have been only partially successful.In the present paper we will try to explain the several odorous properties of single compounds on the basis of the non-specific properties of odorants involved in solubility.This model is a first approach, and although it gives statistically highly significant relations, it is not as accurate as those advanced with respect to the physical and sensory dimensions of stimuli in the fields of vision and audition.We will first give the present definitions of the most suitable physicochemical parameters, and then advance quantitative and qualitative models for single compounds. Quantitative odorous properties are: odour threshold, rate of change of odour intensity with odorant concentration in the suprathreshold region, and the somewhat controversial upper odour intensity. Qualitative properties refer to odour character.     

Xylidine isomers. VII. Structure-hepatotoxic activity relationships of some xylidines

Hepatotoxic studies of xylidines differ with regard to the animal species and type of isomer. The aim of the present paper is to correlate electronic structure data and physical-chemical properties, studied by the authors during a few previous papers (Sahini & Weinberg, 1973; Sahini, Weinberg & Vasilescu 1976, Weinberg & Sahini, to be published; Weinberg & Sahini, 1980), with the hepatotoxic activity of xylidines. Two possible biochemical mechanisms are advanced by help of correlation equations.     

The fine structure of Ameson pulvis (Microspora,Microsporida) and its implications regarding classification and chromosome cycle

Nosema pulvisPerez, 1905, Ameson pulvis (Perez) Sprague, 1977, in muscles of the crabs Carcinus maenas and C. mediterraneus from the coast of France, was observed with the electron microscope. It was found to be structurally similar to the type species A. michaelis (Sprague, 1970). Sprague, 1977, having moniliform sporogonial plasmodia, unikaryotic sporoblasts, and hirsute sporulation stages. It is treated as distinct from A. michaelis because it has slightly smaller spores (by comparison with syntype material of A. michaelis) and appears to have fewer coils in the polar filament. The results require the removal of the genus Ameson from the family Nosematidae Labbé, 1899, where Sprague (1977) had placed it under the erroneous supposition that its sporoblasts are diplokaryotic. Ameson is transferred to family Unikaryonidae Sprague, 1977. Ameson is distinguished from PereziaLéger and Duboscq, 1909, shown by Ormieres et al. to have a similar developmental pattern, by presence of appendages on its sporulation stage. A. nelsoni (Sprague, 1950), the third, and only other species of Ameson, lacks the appendages and is transferred to genus Perezia.     

Trehalose-6-phosphate synthetase activity in extracts of Dictyostelium discoideum.

Trehalose-6-phosphate (T-6-P) synthetase activity in extracts of Dictyostelium discoideum has been reexamined in an effort to resolve discrepancies between the results of previous studies (R. Roth and M. Sussman (1966). Biochim. Biophys. Acta, 122, 225; K. A. Killick and B. E. Wright (1972). J. Biol. Chem., 247, 2967). We find that T-6-P synthetase is not cold sensitive as reported by Killick and Wright (1972), is not present in bacterial-grown vegetative cells (though subject to some modulation by other nutritional conditions), and is not in our hands unmasked or activated by ammonium sulfate fractionation. We conclude that the pattern of T-6-P synthetase accumulation and disappearance during fruiting body construction in D. discoideum is as originally described by R. Roth and M. Sussman (1968). J. Biol. Chem., 243, 5081) and confirmed elsewhere (P. C. Newell et al. (1972). J. Mol. Biol., 63, 373; R. W. Brackenbury et al. (1974). J. Mol. Biol., 90, 529; B. D. Hames and J. M. Ashworth (1974). Biochem. J., 142, 301).     

Check-list of catarrhina primate hybrids

An up-to-date check-list of nonhuman Primate hybrids has important implications in the areas of taxonomy and phylogeny. A list of this nature could, moreover, provide valuable information for further research into the comparative cytogenetics and genetics of Primates (Chiarelli, 1966).About twelve years ago a first check-list of data on Primate hybrids was published (Chiarelli, 1961). The following is a revised list including more recent data. The Table shows the species hybridized, the vitality (survival period) and fertility of the hybrids, the number of hybrids of each type, the location of the hybrid births and the source of the information.     

Kinship and covariance

Price's (1970) covariance theorem can be used to derive an expression for gene frequency change in kin selection models in which the fitness effect of an act is independent of the genotype of the recipient. This expression defines a coefficient of relatedness which subsumes r(Wright, 1922), b(Hamilton, 1972), ρ (Orlove &; Wood, 1978), and R(Michod &; Hamilton, 1980). The new coefficient extends the domain of Hamilton's rule to models in which the average gene frequency of actors differs from that of recipients.     

An autoradiographic analysis of the time of appearance of neurons in the developing chick neural retina

The pattern of incorporation of [³H]thymidine into the chick neural retina has been used to establish the time and order in which different classes of neuroepithelial cells withdraw from the cell cycle and initiate migration and differentiation.The posterior pole of the retina is the first to form during development. In this region most neuroepithelial cells complete mitotic activity between the third and sixth day of incubation. Presumptive ganglion cells initiate the withdrawal process, and they are soon followed by the neuroepithelial precursors of amacrine, horizontal, and receptor cells. Bipolar cell precursors are the last to begin and the last to complete cell cycle activity. It is worthy of note, however, that, in any given region of the retina, neuroepithelial cells of all types cease mitosis in close, overlapping succession.These results are in reasonable agreement with those previously published on the chick retina by Fujita and Horii (1963), and other investigators on the mouse (Mus), killifish (Fundulus), and toad (Xenopus). The present data are also consistent with those proposals of Angevine (1970), Jacobson, 1968a, Jacobson, 1968b, Jacobson, 1970, and others that relate the cessation of mitotic activity of neuroepithelial cells to the determination of neuronal size, axon length, and the specification of neuronal connections.     

Does a functional relationship exist between the polymerization and catalytic activity of beef liver glutamate dehydrogenase?

The recent work of Cohen &; Benedek (1976) and Cohen et al. (1975, 1976) on the apparent interdependence of beef liver glutamate dohydrogenase catalytic activity and degree of polymerization is examined in the light of previously published equilibrium and kinetic results. It is shown that some of the hypotheses central to the Cohen &; Benedek (1976) model are in contradiction with existent data. Consideration of all available information leads to the conclusion that effector-induced depolymerization may simply be an incidental side reaction in the events leading to inhibition.     

A stochastic model of animal growth

The plasticity of growth of animals in time, due to the resilience of their response to the ways they can be fed, suggests the difficulty of describing growth by a stochastic model in the time domain. A model is presented which avoids this difficulty by describing growth as a Markov process in the food-consumed domain, assuming that, at conception, (1) the maximum mature weight as a number α of biomass units of mass μ, and (2) the probability B of production of a biomass unit per unit of food consumed, are specified. Constancy of α, μ and B, as the animal feeds, is the basis of the proposed Markov process. The mean growth from infancy to maturity in the food-consumed domain is then the old law of diminishing returns empirically formulated first by Spillman (1924) for cattle and swine, and confirmed by Titus, Jull &; Hendricks (1934) for fowl, and by Parks (1972) across species from mice to steers. The solution also leads to the possibility that the distribution of weights in a population of growing animals of the same species, is related to the distribution of mature weights among the individuals. An experiment by Lister &; McCance (1967) with well-fed and severely undernourished pigs, shows the stability of growth in the foodconsumed domain compared to the plasticity in the time domain. Other implications of the model are discussed.     

Variance of genetic distance and correlation of heterozygosity between populations under the pressure of stepwise mutation

Using the stepwise mutation model of Ohta and Kimura (1973), formulas are developed for the correlation of heterozygosity and the variance of genetic distance between two finite populations. Studied in detail is the case where the sizes of the two descendant populations are equal to that of the ancestral population and the mutation rate is the same for all loci. Numerical computations are carried out by using the present formulas and those of Li and Nei (1975Genet. Res.25) for the infinite-allele model. The results are as follows: The correlation of heterozygosity decreases with time faster for the stepwise mutation model than for the infinite-allele model. However, the relationships between the correlation of heterozygosity and the normalized genetic identity for the two models are very similar, if the average heterozygosities of the two populations are around 0.20 or less. On the other hand, the variance of genetic distance for the stepwise mutation model may become considerably smaller than that for the infinite-allele model, if the average heterozygosities of the two populations are larger than 0.05. The ratio of the standard deviation to the mean is, however, very large for the stepwise mutation model as well as the infinite-allele model.     

Unified theory of 1f and conductance noise in nerve membrane

A theory of $\text{1}\text{f}$ and conductance noise is given for ionic channels in nerve membrane. The theory is based on the assumption that the channels are in constant, stochastically independent, rotational motion within a fluid bilayer membrane. The resulting expression for the current noise power density S contains a conduction noise term consistent withStevens (1972) and Hill & Chen (1972) and a $\text{1}\text{f}$ noise term consistent with Lundstrom & McQueen (1974) and Clay & Shlesinger (1976). The expression for S also contains a third term which is the spectrum of the product of the single channel conduction noise and $\text{1}\text{f}$ noise correlation functions. This term is independent of the number of channels in the membrane, R. Consequently, the expression for S effectively reduces to a sum of $\text{1}\text{f}$ and conduction noise for R 10–100 which is in agreement with noise measurements on squid axon. The theory is applied in detail to potassium squid noise measurements of Conti, DeFelice & Wanke (1975) using the stochastic analysis of single file ion motion developed in our previous paper (Clay & Shlesinger (1976)).     

The Plasmodium vaughani--Complex

Plasmodium vaughaniNovy and MacNeal, 1904, Plasmodium tenueLaveran and Marullaz, 1914, and Plasmodium merulaeCorradetti and Scanga, 1972 are shown to differ. It is suggested that P. tenue and P. merulae should be considered as subspecies belonging to Plasmodium vaughani-complex.More investigations are needed for a sufficient knowledge of the complex, particularly because at least 36 species of birds harbor P. vaughani-like parasites and cover an immense geographical area in all the parts of the world.     

Local mate competition and the sex ratio

Hamilton (1967) pointed out that Fisher's (1930) argument predicting an equality of the sex ratio may break down when there is local competition for mates. He considered in particular a model in which the environment consists of a number of isolated patches, each of which is colonized by a number of inseminated females; the offspring breed within the patch before dispersal. The present paper provides a careful derivation of the equilibrium sex ratio under this model in both diploid and haplo-diploid populations, and extends the model to consider the effects of having a finite number of patches.We suggest that the equilibrium sex ratio is not simply a function of the amount of inbreeding or sib-mating, as suggested by Maynard Smith (1978), but that the detailed breeding structure of the population must be taken into account.     

A note on the sampling theory for infinite alleles and infinite sites models

This note considers sampling theory for a selectively neutral locus where it is supposed that the data provide nucleotide sequences for the genes sampled. It thus anticipates that technical advances will soon provide data of this form in volume approaching that currently obtained from electrophoresis. The assumption made on the nature of the data will require us to use, in the terminology ofKimura (Theor. Pop. Biol.2, 174–208 (1971)), the “infinite sites” model of Karlin and McGregor (Proc. Fifth Berkeley Symp. Math. Statist. Prob.4, 415–438 (1967)) rather that the “infinite alleles” model of Kimura and Crow (Genetics49, 174–738 (1964)). We emphasize that these two models refer not to two different real-world circumstances, but rather to two different assumptions concerning our capacity to investigate the real world. We compare our results where appropriate with corresponding sampling theory of Ewens (Theor. Pop. Biol.3, 87–112 (1972)) for the “infinite alleles” model. Note finally that some of our results depend on an assumption of independence of behavior at individual sites; a parallel paper byWatterson (submitted for publication (1974)) assumes no recombination between sites. Real-world behavior will lie between these two assumptions, closer to the situation assumed by Watterson than in this note. Our analysis provides upper bounds for increased efficiency in using complete nucleotide sequences.     

Co-operative response of chemically excitable membrane: I. Formulation: Unified theory of co-operativity

The lattice-model of Changeux, Thiery, Tung & Kittel (1966) was extended in order to examine the co-operative response of chemically excitable membrane and the exact mathematical correspondence to the Ising model was shown. In this model, two conformational states S and R with different affinities for the ligand are assumed to be accessible to each protomer, which is interacting with the nearest-neighbor protomers. The model is applicable to any kind of symmetrically interacting system consisting of oligomers and lattices and is an extension of previously proposed models of allosteric protein. It includes the model of Monod, Wyman, & Changeux (1965) and that of Koshland, Némethy & Filmer (1966) as the extreme cases of the oligomer. By assuming that a state-transition from S to R in a protomer is accompanied by a unit increase in conductance, the characteristics of dose-response curves of chemically excitable membrane are examined. The Hill's coefficient n_H of dose-response curve, the measure of the co-operativity, is shown to be proportional to the square of the mean fluctuation of the state function, the fraction of protomers in R state.     

One- and two-locus models of the origin of worker behavior in hymenoptera

Hamilton, 1964, Hamilton, 1972 suggested that, if females are singly inseminated, the asymmetries in intra-familial relationships inherent in haplo-diploidy should predispose the Hymenoptera to evolve female sex-limited worker behavior, if either a worker raises sibs as reproductives but raises a bias of sisters (pathway 1), or she raises a combination of sibs and sons as reproductives but substitutes sons for brothers preferentially (pathway 2). To investigate these hypotheses when mating is random, we present mathematical and numerical analyses of the one-locus two-allele model of Charnov (1978a) which we generalize. Furthermore, we extend the Charnov model to a more realistic two-locus situation where we postulate that the ability of a worker to raise a bias of sisters or to substitute sons for brothers preferentially is controlled by a locus separate from the locus controlling worker behavior, rather than being a pleiotropic effect of that locus. Considerations of the relative stringencies of the conditions for substitution of the allele for worker behavior and the constraints on linkage suggest that pathway 2 is more likely than pathway 1 for the origin of worker behavior. For pathway 1, we compare the Charnov model with its diplo-diploid analogue. This comparison reveals that, depending on the penetrance of the allele for worker behavior in the heterozygous state, both models permit the evolution of worker behavior even when workers raise the ratio of sisters to brothers laid by the queens, and that an advantage to haplo-diploidy becomes apparent only when workers raise a bias of sisters. For pathway 1, we also consider a generalization of the Craig model (Craig, 1979) which is quite similar to the Charnov model.     

The myosin filament: VII. Changes in internal structure along the length of the filament

Improved fixation procedures have enabled substructure to be observed by electron microscopy in transverse sections of vertebrate skeletal muscle thick filaments as thin as 140 nm. Optical diffraction combined with digital autocorrelation analysis, focal series and tilting experiments have confirmed the presence of a regular substructure having a repeat near 4 nm and shown that it is highly unlikely to be an artifact associated with the electron microscope imaging system. The results obtained strongly suggest that the thick filament is constructed from a bundle of rod-like subfilaments arranged parallel to the thick filament axis to within less than a degree. This cannot easily be reconciled with the general theory of thick filament structure proposed by Squire (1973), but it is consistent with the model proposed by Pepe, 1966, Pepe, 1967. Optical diffraction of 140 nm thick serial transverse sections has also suggested a structural change along the length of the filament that is manifest by a variation in the proportion of filaments showing strong diffraction maxima in one, two or three directions.