Genomic GC level, optimal growth temperature, and genome size in prokaryotes

Two years ago, we showed that positive correlations between optimal growth temperature (T(opt)) and genome GC are observed in 15 out of the 20 families of prokaryotes we analyzed, thus indicating that "T(opt) is one of the factors that influence genomic GC in prokaryotes". Our results were disputed, but these criticisms were demonstrated to be mistaken and based on misconceptions. In a recent report, Wang et al. [H.C. Wang, E. Susko, A.J. Roger, On the correlation between genomic G+C content and optimal growth temperature in prokaryotes: data quality and confounding factors, Biochem. Biophys. Res. Commun. 342 (2006) 681-684] criticize our results by stating that "all previous simple correlation analyses of GC versus temperature have ignored the fact that genomic GC content is influenced by multiple factors including both intrinsic mutational bias and extrinsic environmental factors". This statement, besides being erroneous, is surprising because it applies in fact not to ours but to the authors' article. Here, we rebut the points raised by Wang et al. and review some issues that have been a matter of debate, regarding the influence of environmental factors upon GC content in prokaryotes. Furthermore, we demonstrate that the relationship that exists between genome size and GC level is valid for aerobic, facultative, and microaerophilic species, but not for anaerobic prokaryotes.     

Molecular systematics of the Hyaenidae: relationships of a relictual lineage resolved by a molecular supermatrix

The four extant species of hyenas (Hyaenidae; Carnivora) form a morphologically and ecologically heterogeneous group of feliform carnivorans that are remnants of a formerly diverse group of mammalian predators. They include the aardwolf (Proteles cristatus), a termite-feeding specialist, and three species with a craniodental morphology adapted to cracking the bones of prey and/or carcasses, the spotted hyena (Crocuta crocuta), brown hyena (Parahyaena brunnea), and striped hyena (Hyaena hyaena). Hyenas have been the subject of a number of systematic studies during the last two centuries, due in large part to the extensive fossil record of the group, with nearly 70 described fossil species. Morphological studies incorporating both fossil and living taxa have yielded different conclusions regarding the evolutionary relationships among living hyenas. We used a molecular supermatrix comprised of seven nuclear gene segments and the complete mitochondrial cytochrome b gene to evaluate phylogenetic relationships among the four extant hyaenid species. We also obtained sequence data from representative species of all the main families of the Feliformia (Felidae, Herpestidae, and Viverridae) to estimate the sister group of the Hyaenidae. Maximum parsimony and maximum likelihood analyses of the supermatrix recovered identical topologies. Furthermore, Bayesian phylogenetic analyses of the supermatrix, with among-site rate variation among data partitions parameterized in three different ways, also yielded the same topology. For each phylogeny reconstruction method, all but two nodes received 100% bootstrap or 1.00 posterior probability nodal support. Within the monophyletic Hyaenidae, Parahyaena and Hyaena were joined together, with Crocuta as the sister to this clade, and Proteles forming the most basal lineage. A clade containing two species of mongoose (core Herpestidae) plus Cryptoprocta ferox (currently classified in Viverridae) was resolved as the sister group of Hyaenidae. The pattern of relationships among the three bone-cracking hyaenids (Crocuta, Hyaena, and Parahyaena) is incongruent with recent cladistic assessments based on morphology and suggests the need to reevaluate some of the morphological characters that have been traditionally used to evaluate relationships among hyenas. Divergence time estimates based on a Bayesian relaxed molecular clock indicates that hyaenids diverged from their feliform sister group 29.2 MYA, in the Middle Oligocene. Molecular clock estimates also suggest that the origin of the aardwolf is much more recent (10.6 MYA) than that implied by a cladistic analysis of morphology ( approximately 20 MYA) and suggests that the aardwolf is possibly derived from a bone and meat eating lineage of hyaenids that were present in the Late Miocene. [Hyaenidae; phylogeny; cytochrome b; nuclear gene segments; Proteles; Crocuta; Hyaena; Parahyaena.].     

Computational tools for isotopically instationary 13C labeling experiments under metabolic steady state conditions

(13)C metabolic flux analysis (MFA) has become an important and powerful tool for the quantitative analysis of metabolic networks in the framework of metabolic engineering. Isotopically instationary (13)C MFA under metabolic stationary conditions is a promising refinement of classical stationary MFA. It accounts for the experimental requirements of non-steady-state cultures as well as for the shortening of the experimental duration. This contribution extends all computational methods developed for classical stationary (13)C MFA to the instationary situation by using high-performance computing methods. The developed tools allow for the simulation of instationary carbon labeling experiments (CLEs), sensitivity calculation with respect to unknown parameters, fitting of the model to the measured data, statistical identifiability analysis and an optimal experimental design facility. To explore the potential of the new approach all these tools are applied to the central metabolism of Escherichia coli. The achieved results are compared to the outcome of the stationary counterpart, especially focusing on statistical properties. This demonstrates the specific strengths of the instationary method. A new ranking method is proposed making both an a priori and an a posteriori design of the sampling times available. It will be shown that although still not all fluxes are identifiable, the quality of flux estimates can be strongly improved in the instationary case. Moreover, statements about the size of some immeasurable pool sizes can be made.     

Effects of sex control and twinning on economic optimization of culling cows in Japanese Black cow-calf production systems

The effects of sex control and twinning techniques on determination of optimal culling parity of cows in beef cow-calf production systems were deterministically analyzed using a herd model simulation. The model simulated the annualized net revenue as an economic indicator during the whole life cycle of a cow. Biological factors (survivability, growth, reproduction, and feed requirements) and economic factors (returns from sales of live calves and cows' carcasses and production costs) were included in the model. Some biological and economic parameters relating to these factors were altered from a base condition in order to adapt the production systems with sex control and twinning techniques. Based on the model, early culling was optimal for all production systems when biological efficiency was used as an indicator of production; however, later culling was optimal for single production, but slightly earlier culling was optimal for twin production, when annualized net revenue was evaluated. The introduction of sex control did not greatly affect the determination of the optimal culling parity of cows. When production included the sex control, female sexing increased biological efficiency, whereas male sexing increased annualized net revenue. In the present beef cow-calf production circumstances in Japan, introduction of sex control did not have economically appreciable effects, but twinning was economically beneficial. For production involving sex control, improvement in the conception rate per mating and/or reduction of technical cost were required for this technology to be profitable.     

The optimal burst of mutation to create a phenotype

Mutagenesis is commonly applied to genes and genomes to create novel variants with desired properties. This paper calculates the level of mutagenesis that maximizes the appearance of favorable mutants, assuming that the mutagenesis is applied in a single episode. The downside of mutagenesis is that a substantial fraction of mutations will destroy gene/genome function. The upside of mutagenesis is the production of beneficial mutations, but the desired phenotype may require that 1, 2 or more beneficial mutations be present simultaneously (the phenotype dimensionality). The optimum level of mutagenesis is sensitive to both properties. In the simplest model, the mutation optimum occurs when number of lethal equivalents per genome equals the phenotype dimensionality, a result first derived by Mundry and Gierer [1958. Production of mutations in tobacco mosaic virus by chemical treatment of its nucleic acid in vitro. Z. Vererbungsl. 89 (4), 614-630]. This level of mutation is shown to be an upper bound for the optimum in various extensions of the model, and the recovery of mutants is also reasonably tolerant to deviations from the optimum.     

Expiration rate drives human airway design

Analyses of human airway architecture based on calculations of airflow resistance or energy dissipation suggest that the branching pattern is not optimized for minimizing energy loss by flow dissipation during respiration. Airway flow dissipates only a few percent of the total body work during normal breathing, so branching patterns deviate from minimum energy loss to also optimize other physiological needs. Studies of airway performance often record some measure of expiration, such as FEV1 (Forced Expiratory Volume in 1 s), because airway constriction during expiration limits the rate of rapid respiration. We posit that lung structure is optimized for the rate of expiration as well as minimum energy loss. By increasing the daughter-to-parent airway diameter ratio (h) from 0.794 (corresponding to the energy minimum for symmetrically branching airways) to 0.85 (the observed value in humans) luminal pressures at airway generations 4-15 were substantially increased during exercise (a 4.5 and 15 cmH₂O increase during moderate and heavy exercise, respectively). Values of h somewhat larger than 0.794 help airways remain open during expiration by increasing both viscous pressure drop and convective acceleration pressure drop. Asymmetric bifurcations also exhibit higher proximal airway pressures than symmetric ones, but the improvement was not large.     

Mutation and control of the human immunodeficiency virus

We examine the dynamics of infection by the human immunodeficiency virus (HIV), as well as therapies that minimize viral load, restore adaptive immunity, and use minimal dosage of anti-HIV drugs. Virtual therapies for wild-type infections are demonstrated; however, the HIV infection is never cured, requiring continued treatment to keep the condition in remission. With high viral turnover and mutation rates, drug-resistant strains of HIV evolve quickly. The ability of optimal therapy to contain drug-resistant strains is shown to depend upon the relative fitness of mutant strains.     

Explaining the optimality of U-shaped age-specific mortality

Mortality is U-shaped with age for many species, declining from birth to sexual maturity, then rising in adulthood, sometimes with postreproductive survival. We show analytically why the optimal life history of a species with determinate growth is likely to have this shape. An organism allocates energy among somatic growth, fertility and maintenance/survival at each age. Adults may transfer energy to juveniles, who can then use more energy than they produce. Optimal juvenile mortality declines from birth to maturity, either to protect the increasingly valuable cumulative investments by adults in juveniles or to exploit the compounding effects of early investment in somatic growth, since early growth raises subsequent energy production, which in turn supports further growth. Optimal adult mortality rises after maturity as expected future reproduction declines as in Hamilton, but intergenerational transfers lead to postreproductive survival as in Lee. Here the Hamilton and transfer effects are divided by probabilities of survival in contrast to the fitness impact measures, which are relevant for mutation-selection balance. If energetic efficiency rises strongly with adult experience, then adult mortality could initially be flat or declining.     

Allee Effect and Control of Lake System Invasion

We consider the model of invasion prevention in a system of lakes that are connected via traffic of recreational boats. It is shown that in presence of an Allee effect, the general optimal control problem can be reduced to a significantly simpler stationary optimization problem of optimal invasion stopping. We consider possible values of model parameters for zebra mussels. The general N-lake control problem has to be solved numerically, and we show a number of typical features of solutions: distribution of control efforts in space and optimal stopping configurations related with the clusters in lake connection structure.     

Optimal fed-batch operation of recombinant cells subject to plasmid instability and death

Optimal feed rate strategy is studied for fed-batch culture of recombinant cells with plasmid instability and with different death rates for the plasmid-free cells (PFC) and plasmid-bearing cells (PBC). Most of the fed-batch fermentation is known to have first-order singularity and therefore a single singular arc. However, this study shows that a singular arc with second-order singularity and therefore two distinct singular arcs are possible for a recombinant cell process if PFC and PBC are subjected to death, and their specific growth rates are proportional to each other. Two types of singular arcs are elucidated and analyzed. The optimal policies over the singular arcs are theoretically explored as these findings reveal qualitative information on the singular arc, which is critically important in providing the optimal initial conditions in numerical computation of optimal feed rate profile.