Revisiting Offspring Maxima in Branching Processes

2000 Mathematics Subject Classification: 60J80; 60G70.

Extremes of Bivariate Geometric Variables with Application to Bisexual Branching Processes

2000 Mathematics Subject Classification: 60J80, 60G70.

Conservation of dual-targeted proteins in Arabidopsis and rice points to a similar pattern of gene-family evolution

Gene duplication followed by acquisition of specific targeting information and dual targeting were evolutionary strategies enabling organelles to cope with overlapping functions. We examined the evolutionary trend of dual-targeted single-gene products in Arabidopsis and rice genomes. The number of paralogous proteins encoded by gene families and the dual-targeted orthologous proteins were analysed. The number of dual-targeted proteins and the corresponding gene-family sizes were similar in Arabidopsis and rice irrespective of genome sizes. We show that dual targeting of methionine aminopeptidase, monodehydroascorbate reductase, glutamyl-tRNA synthetase, and tyrosyl-tRNA synthetase was maintained despite occurrence of whole-genome duplications in Arabidopsis and rice as well as a polyploidization followed by a diploidization event (gene loss) in the latter.

Family size, economics and child gender preference: A case study in the Nyeri district of Kenya

Kenyan women have more children, especially in rural areas, than in most developing nations. This is widely believed to be an impediment to Kenya’s economic development. Thus, factors influencing family size in the Kenyan context are important for its future. A brief review of economic theories of fertility leads to the conclusion that both economics and social/cultural factors must be considered simultaneously when examining factors that determine the number of children in a family. The need to do this is borne out in Kenya’s situation by utilising responses from a random sample of rural households in the Nyeri district of Kenya. Economic and social/cultural factors intertwine to influence family sizes in this district. After providing a summary of the main statistical results from the survey, we use multiple regression analysis to explore the influences of a woman’s age, level of education, whether she has outside employment, whether the family keeps livestock, whether she expresses a preference for more boys than girls, whether the family uses only family labour (including child labour) and the size of the farm, which is used as a proxy for family income. It was found that preference for male children has an important positive influence on family size in this district. Women were found to have greater preference for male children than their male counterparts possibly because of their fear of being disinherited if they do not produce an heir for their husbands. Preference for sons was also found in allocation of human capital resources at the household level in that the female respondents were found to have lower levels of education than their male counterparts. Various long-term policies are outlined that may help to reduce the number of offspring of women in Kenya.

BadiRate: Estimating Family Turnover Rates by Likelihood-Based Methods

Motivation: The comparative analysis of gene gain and loss rates is critical for understanding the role of natural selection and adaptation in shaping gene family sizes. Studying complete genome data from closely related species allows accurate estimation of gene family turnover rates. Current methods and software tools, however, are not well designed for dealing with certain kinds of functional elements, such as microRNAs or transcription factor binding sites. Results: Here, we describe BadiRate, a new software tool to estimate family turnover rates, as well as the number of elements in internal phylogenetic nodes, by likelihood-based methods and parsimony. It implements two stochastic population models, which provide the appropriate statistical framework for testing hypothesis, such as lineage-specific gene family expansions or contractions. We have assessed the accuracy of BadiRate by computer simulations, and have also illustrated its functionality by analyzing a representative empirical dataset.

BadiRate: estimating family turnover rates by likelihood-based methods

Motivation: The comparative analysis of gene gain and loss rates is critical for understanding the role of natural selection and adaptation in shaping gene family sizes. Studying complete genome data from closely related species allows accurate estimation of gene family turnover rates. Current methods and software tools, however, are not well designed for dealing with certain kinds of functional elements, such as microRNAs or transcription factor binding sites. Results: Here, we describe BadiRate, a new software tool to estimate family turnover rates, as well as the number of elements in internal phylogenetic nodes, by likelihood-based methods and parsimony. It implements two stochastic population models, which provide the appropriate statistical framework for testing hypothesis, such as lineage-specific gene family expansions or contractions. We have assessed the accuracy of BadiRate by computer simulations, and have also illustrated its functionality by analyzing a representative empirical dataset.

Family Size Evolution in Drosophila Chemosensory Gene Families: A Comparative Analysis with a Critical Appraisal of Methods

Gene turnover rates and the evolution of gene family sizes are important aspects of genome evolution. Here, we use curated sequence data of the major chemosensory gene families from Drosophila-the gustatory receptor, odorant receptor, ionotropic receptor, and odorant-binding protein families-to conduct a comparative analysis among families, exploring different methods to estimate gene birth and death rates, including an ad hoc simulation study. Remarkably, we found that the state-of-the-art methods may produce very different rate estimates, which may lead to disparate conclusions regarding the evolution of chemosensory gene family sizes in Drosophila. Among biological factors, we found that a peculiarity of D. sechellia's gene turnover rates was a major source of bias in global estimates, whereas gene conversion had negligible effects for the families analyzed herein. Turnover rates vary considerably among families, subfamilies, and ortholog groups although all analyzed families were quite dynamic in terms of gene turnover. Computer simulations showed that the methods that use ortholog group information appear to be the most accurate for the Drosophila chemosensory families. Most importantly, these results reveal the potential of rate heterogeneity among lineages to severely bias some turnover rate estimation methods and the need of further evaluating the performance of these methods in a more diverse sampling of gene families and phylogenetic contexts. Using branch-specific codon substitution models, we find further evidence of positive selection in recently duplicated genes, which attests to a nonneutral aspect of the gene birth-and-death process.

Downhole variations in grain sizes at DSDP/ODP Hole 504B

The maximum grain sizes of plagioclase and magnetite in the groundmass of the sheeted dike complex drilled at Hole 504B have been measured. Downhole variations through a 440-m-long section show a crude zig-zag pattern consisting of a gradual decrease or increase followed by an abrupt jump. The gradual decrease or increase in grain size extends over many lithologic units, and hence, does not reflect variations in grain size within a single dike. Such a zig-zag pattern is well explained by grain-size variations through multiple dikes. By using the observed inclination of sheeted dikes of 81° ± 2.5°, thickness of the multiple dikes varies from 0.7 to 8.5 m and averages to 4 ± 1 m. The average thickness of individual dikes forming multiple dikes is 0.8 m. We expect such multiple dikes to be formed during rifting events beneath mid-oceanic spreading ridges. If the average expansion at rifting episodes is twice as wide as the average width of the multiple dike units, the full spreading rate of 7.2 cm/yr of Cocos Ridge gives 112 ± 33 yr for a time interval of the rifting. A simple one-dimensional conductive cooling model is applied to solidification of multiple dikes. Numerical simulations show that the grain-size variations observed through the drill hole are more consistent with a model where a new injection of a dike occurs periodically with a constant time interval rather than one where the next dike intrudes just after the solidification of the previous one. Grain-size variations within simple dikes from Iritono, Japan, and those for Makaopuhi lava lake, Hawaii, show that square root of crystallization time is linearly correlated with the logarithm of plagioclase size. By using an empirically derived relationship between these two variables, the variations of plagioclase size through Hole 504B are directly compared with the calculated times for crystallization. Each rifting episode at the Costa Rica Rift lasts for several years, and periodic injection of a new dike occurs into the center of a previously solidified multiple dike at time intervals varying from 1 to 12 months.

Use of the EM algorithm to detect QTL affecting multiple-traits in an across half-sib family analysis

QTL detection experiments in livestock species commonly use the half-sib design. Each male is mated to a number of females, each female producing a limited number of progeny. Analysis consists of attempting to detect associations between phenotype and genotype measured on the progeny. When family sizes are limiting experimenters may wish to incorporate as much information as possible into a single analysis. However, combining information across sires is problematic because of incomplete linkage disequilibrium between the markers and the QTL in the population. This study describes formulae for obtaining MLEs via the expectation maximization (EM) algorithm for use in a multiple-trait, multiple-family analysis. A model specifying a QTL with only two alleles, and a common within sire error variance is assumed. Compared to single-family analyses, power can be improved up to fourfold with multi-family analyses. The accuracy and precision of QTL location estimates are also substantially improved. With small family sizes, the multi-family, multi-trait analyses reduce substantially, but not totally remove, biases in QTL effect estimates. In situations where multiple QTL alleles are segregating the multi-family analysis will average out the effects of the different QTL alleles.

Impact of different pack sizes of paracetamol in the United Kingdom and Ireland on intentional overdoses: a comparative study

BACKGROUND: In order to reduce fatal self-poisoning legislation was introduced in the UK in 1998 to restrict pack sizes of paracetamol sold in pharmacies (maximum 32 tablets) and non-pharmacy outlets (maximum 16 tablets), and in Ireland in 2001, but with smaller maximum pack sizes (24 and 12 tablets). Our aim was to determine whether this resulted in smaller overdoses of paracetamol in Ireland compared with the UK. METHODS: We used data on general hospital presentations for non-fatal self-harm for 2002-2007 from the Multicentre Study of Self-harm in England (six hospitals), and from the National Registry of Deliberate Self-harm in Ireland. We compared sizes of overdoses of paracetamol in the two settings. RESULTS: There were clear peaks in numbers of non-fatal overdoses, associated with maximum pack sizes of paracetamol in pharmacy and non-pharmacy outlets in both England and Ireland. Significantly more pack equivalents (based on maximum non-pharmacy pack sizes) were used in overdoses in Ireland (mean 2.63, 95% CI 2.57-2.69) compared with England (2.07, 95% CI 2.03-2.10). The overall size of overdoses did not differ significantly between England (median 22, interquartile range (IQR) 15-32) and Ireland (median 24, IQR 12-36). CONCLUSIONS: The difference in paracetamol pack size legislation between England and Ireland does not appear to have resulted in a major difference in sizes of overdoses. This is because more pack equivalents are taken in overdoses in Ireland, possibly reflecting differing enforcement of sales advice. Differences in access to clinical services may also be relevant.

Response to Heterologous Leishmanins in Cutaneous Leishmaniasis in Nigeria: Discovery of a New Focus

A pilot study was undertaken to preliminary illustrate the leishmanin skin test (LST) positivity to distinct antigen preparations (derived from promastigote of either Leishmania major or L. amazonensis, or pooled L. mexicana, L. amazonensis and L. guyanensis) in cutaneous leishmaniasis (CL) patients and healthy subjects living in two endemic foci in Nigeria. The study was designed to provide insights into whether cross-species leishmanin, such as that prepared from New World Leishmania could be useful to detect cases of Old World leishmanial infection and to compare the results with LST using L. major-derived leishmanin. The overall LST positivity in individuals from Keana tested with the cross-species leishmanin was 28.7% (27/94), while the positivity rate in the subjects from Kanana tested with the same leishmanin was 54.5% (6/11). Lower positivity values were obtained when L. major (12.5%; 11/88) or L. amazonensis (15.8%; 9/57) was tested as antigen in grossly comparable populations. Moreover, the pooled leishmanin identified most of the subjects (13/14; 92.9%) with active or healed CL, and the maximum reaction sizes were found among positive subjects in this group. No healthy controls (10 total) showed specific DTH response. The LST was useful for assessing the prevalence of subclinical infection and for measuring CL transmission over time. We report for the first time the occurrence of CL in Kanana village of Langtang South local government area of Plateau State

A Laboratory Investigation into the Effects of the Aggregated Related Facts of Critical VMA in Asphalt Paving Mixtures, June 2000

This report summarizes research conducted at Iowa State University on behalf of the Iowa Department of Transportation, focusing on the volumetric state of hot-mix asphalt (HMA) mixtures as they transition from stable to unstable configurations. This has raditionally been addressed during mix design by meeting a minimum voids in the mineral aggregate (VMA) requirement, based solely upon the nominal maximum aggregate size without regard to other significant aggregate-related properties. The goal was to expand the current specification to include additional aggregate properties, e.g., fineness modulus, percent crushed fine and coarse aggregate, and their interactions. The work was accomplished in three phases: a literature review, extensive laboratory testing, and statistical analysis of test results. The literature review focused on the history and development of the current specification, laboratory methods of identifying critical mixtures, and the effects of other aggregate-related factors on critical mixtures. The laboratory testing involved three maximum aggregate sizes (19.0, 12.5, and 9.5 millimeters), three gradations (coarse, fine, and dense), and combinations of natural and manufactured coarse and fine aggregates. Specimens were compacted using the Superpave Gyratory Compactor (SGC), conventionally tested for bulk and maximum theoretical specific gravities and physically tested using the Nottingham Asphalt Tester (NAT) under a repeated load confined configuration to identify the transition state from sound to unsound. The statistical analysis involved using ANOVA and linear regression to examine the effects of identified aggregate factors on critical state transitions in asphalt paving mixtures and to develop predictive equations. The results clearly demonstrate that the volumetric conditions of an HMA mixture at the stable unstable threshold are influenced by a composite measure of the maximum aggregate size and gradation and by aggregate shape and texture. The currently defined VMA criterion, while significant, is seen to be insufficient by itself to correctly differentiate sound from unsound mixtures. Under current specifications, many otherwise sound mixtures are subject to rejection solely on the basis of failing to meet the VMA requirement. Based on the laboratory data and statistical analysis, a new paradigm to volumetric mix design is proposed that explicitly accounts for aggregate factors (gradation, shape, and texture).

Effects of allometry, productivity and lifestyle on rates and limits of body size evolution

Body size affects nearly all aspects of organismal biology, so it is important to understand the constraints and dynamics of body size evolution. Despite empirical work on the macroevolution and macroecology of minimum and maximum size, there is little general quantitative theory on rates and limits of body size evolution. We present a general theory that integrates individual productivity, the lifestyle component of the slow–fast life-history continuum, and the allometric scaling of generation time to predict a clade's evolutionary rate and asymptotic maximum body size, and the shape of macroevolutionary trajectories during diversifying phases of size evolution. We evaluate this theory using data on the evolution of clade maximum body sizes in mammals during the Cenozoic. As predicted, clade evolutionary rates and asymptotic maximum sizes are larger in more productive clades (e.g. baleen whales), which represent the fast end of the slow–fast lifestyle continuum, and smaller in less productive clades (e.g. primates). The allometric scaling exponent for generation time fundamentally alters the shape of evolutionary trajectories, so allometric effects should be accounted for in models of phenotypic evolution and interpretations of macroevolutionary body size patterns. This work highlights the intimate interplay between the macroecological and macroevolutionary dynamics underlying the generation and maintenance of morphological diversity.

Chemical analysis of ash layers in ODP Leg 104 holes

Numerous fresh ash layers comprise about 0.3% by volume of Neogene to Holocene sediments drilled at Leg 104 Sites 642 and 643 (Vøring Plateau, North Atlantic). Median grain sizes of the ashes are about 100 /µm and maximum grain sizes range up to 1200 µm. Rhyolitic pumice shards dominate, with minor bubble wall shards. Basaltic shards are poorly vesicular and blocky or round. Phenocrystic plagioclase, zircon, and clinopyroxene occur in the rhyolitic, plagioclase, and clinopyroxene phenocrysts and basaltic lithics in the basaltic tephra. Quartz, amphibole, clinozoisite, and rutile are interpreted as xenocrysts. All ash layers are well-sorted and represent distal fallout from major explosive eruptions. Most ashes are rhyolitic (high-K and low-K) in composition, some are bimodal (tholeiitic and rhyolitic). Early Miocene tephra is dominantly basaltic. Iceland is inferred to be the likely source region for most ashes. Late Miocene high-K rhyolites may have originated from the K-rich Jan Mayen magmatic province. One Quaternary layer with biotite and alkali feldspar phenocrysts may have been derived from Jan Mayen Island. Four individual Pliocene to Holocene ash layers from Sites 642 and 643 can be correlated fairly well. Upper Miocene layers are tentatively correlated as a sequence between Sites 642 and 643. Average calculated layer frequencies are about three layers/m.y. through the Pliocene and Pleistocene and five to eight layers per m.y. through the middle and late Miocene, suggesting rather continuous volcanic activity in the North Atlantic. Episodic magmatic activity during Neogene epochs in this part of the North Atlantic, as postulated in the literature, cannot be confirmed.

Rutting-induced permeability loss of open graded friction course mixtures

Functionality of an open graded friction course (OGFC) depends on the high interconnected air voids or pores of the OGFC mixture. The authors' previous study indicated that the pores in the OGFC mixture were easily clogged by rutting deformation. Such a deformation-related clogging can cause a significant rutting-induced permeability loss in the OGFC mixture. The objective of this study was to control and reduce the rutting-induced permeability loss of the OGFC based on mixture design and layer thickness. Eight types of the OGFC mixtures with different air void contents, gradations, and nominal maximum aggregate sizes were fabricated in the laboratory. Wheel-tracking rutting tests were conducted on the OGFC slabs to simulate the deformation-related clogging. Permeability tests after different wheel load applications were performed on the rutted OGFC slabs using a falling head permeameter developed in the authors' previous study. The relationships between permeability loss and rutting depth as well as dynamic stability were developed based on the eight OGFC mixtures' test results. The thickness effects of the single-layer and the two-layer OGFC slabs were also discussed in terms of deformation-related clogging and the rutting-induced permeability loss. Results showed that the permeability coefficient decreases linearly with an increasing rutting depth of the OGFC mixtures. Rutting depth was recommended as a design index to control permeability loss of the OGFC mixture rather than the dynamic stability. Permeability loss due to deformation-related clogging can be effectively reduced by using a thicker single-layer OGFC or two-layer OGFC.