普通野生稻线粒体蛋白质编码基因密码子使用偏好性的分析

以普通野生稻(Oryza rufipogon Griff.)线粒体基因组为对象,分析其蛋白质编码基因的密码子使用特征及与亚洲栽培稻(O. sativa L.)的差异,探讨其密码子偏性形成的影响因素和进化过程。结果显示:普通野生稻线粒体基因组编码序列第1、第2和第3位碱基的GC含量依次为49.18%、42.67%和40.86%;有效密码子数(Nc)分布于45.32～61.00之间,其密码子偏性较弱; Nc值仅与GC_3呈显著相关,密码子第3位的碱基组成对密码子偏性影响较大;第1向量轴上显示9.91%的差异,其与GC_3s、Nc、密码子偏好指数(CBI)和最优密码子使用频率(Fop)的相关性均达到显著水平;而GC_3和GC₁₂的相关性未达到显著水平。因此,普通野生稻线粒体基因组密码子的使用偏性主要受自然选择压力影响而形成。本研究确定了21个普通野生稻线粒体基因组的最优密码子,大多以A或T结尾,与叶绿体密码子具有趋同进化,但是与核基因组具有不同的偏好性。同义密码子相对使用度(RSCU)、PR2偏倚分析和中性绘图分析显示,普通野生稻线粒体基因功能和其密码子使用密切相关,且线粒体密码子使用在普通野生稻、粳稻(O. sativa L. subsp. japonica Kato)和籼稻(O. sativa L. subsp.indica Kato)内具有同质性。     

Sex-specific expression of a HOX gene associated with rapid morphological evolution

Animal diversity is shaped by the origin and diversification of new morphological structures. Many examples of evolutionary innovations are provided by male-specific traits involved in mating and sexual selection. The origin of new sex-specific characters requires the evolution of new regulatory interactions between sex-determining genes and genes that control spatial patterning and cell differentiation. Here, we show that sex-specific regulation of the HOX gene Sex combs reduced (Scr) is associated with the origin and evolution of the Drosophila sex comb — a novel and rapidly diversifying male-specific organ. In species that primitively lack sex combs, Scr expression shows little spatial modulation, whereas in species that have sex combs, Scr is upregulated in the presumptive sex comb region and is frequently sexually dimorphic. Phylogenetic analysis shows that sex-specific regulation of Scr has been gained and lost multiple times in Drosophila evolution and correlates with convergent origin of similar sex comb morphologies in several independent lineages. Some of these transitions occurred on microevolutionary timescales, indicating that HOX gene expression can evolve with surprising ease. This is the first example of a sex-specific regulation of a HOX gene contributing to the development and evolution of a secondary sexual trait.     

Synthetic biology and metabolic engineering of actinomycetes for natural product discovery

Actinomycetes are one of the most valuable sources of natural products with industrial and medicinal importance. After more than half a century of exploitation, it has become increasingly challenging to find novel natural products with useful properties as the same known compounds are often repeatedly re-discovered when using traditional approaches. Modern genome mining approaches have led to the discovery of new biosynthetic gene clusters, thus indicating that actinomycetes still harbor a huge unexploited potential to produce novel natural products. In recent years, innovative synthetic biology and metabolic engineering tools have greatly accelerated the discovery of new natural products and the engineering of actinomycetes. In the first part of this review, we outline the successful application of metabolic engineering to optimize natural product production, focusing on the use of multi-omics data, genome-scale metabolic models, rational approaches to balance precursor pools, and the engineering of regulatory genes and regulatory elements. In the second part, we summarize the recent advances of synthetic biology for actinomycetal metabolic engineering including cluster assembly, cloning and expression, CRISPR/Cas9 technologies, and chassis strain development for natural product overproduction and discovery. Finally, we describe new advances in reprogramming biosynthetic pathways through polyketide synthase and non-ribosomal peptide synthetase engineering. These new developments are expected to revitalize discovery and development of new natural products with medicinal and other industrial applications.     

Compliance subsampling designs for comparative research: estimation and optimal planning

For studies with treatment noncompliance, analyses have been developed recently to better estimate treatment efficacy. However, the advantage and cost of measuring compliance data have implications on the study design that have not been as systematically explored. In order to estimate better treatment efficacy with lower cost, we propose a new class of compliance subsampling (CSS) designs where, after subjects are assigned treatment, compliance behavior is measured for only subgroups of subjects. The sizes of the subsamples are allowed to relate to the treatment assignment, the assignment probability, the total sample size, the anticipated distributions of outcome and compliance, and the cost parameters of the study. The CSS design methods relate to prior work (i) on two-phase designs in which a covariate is subsampled and (ii) on causal inference because the subsampled postrandomization compliance behavior is not the true covariate of interest. For each CSS design, we develop efficient estimation of treatment efficacy under binary outcome and all-or-none observed compliance. Then we derive a minimal cost CSS design that achieves a required precision for estimating treatment efficacy. We compare the properties of the CSS design to those of conventional protocols in a study of patient choices for medical care at the end of life.     

Likelihood-based experimental design for estimation of ED50

In selecting the best dosage choice for the estimation of ED50, it is natural to try to minimize the length of the confidence intervals. In this presentation, the dose allocation that minimizes the length of the likelihood-based confidence intervals is presented and compared with alternative allocations that have been proposed based on the length of different types of confidence intervals, such as those based on the asymptotic variance or on Fieller's Theorem. Effective strategies to deal with the parameter dependence of these allocations are explored. A series of experiments to evaluate the effect of small doses per fraction on the radiation tolerance of the rat cervical spinal cord provide the motivation and an illustration for the proposed procedures.     

Optimal HIV treatment by maximising immune response

We present an optimal control model of drug treatment of the human immunodeficiency virus (HIV). Our model is based upon ordinary differential equations that describe the interaction between HIV and the specific immune response as measured by levels of natural killer cells. We establish stability results for the model. We approach the treatment problem by posing it as an optimal control problem in which we maximise the benefit based on levels of healthy CD4+ T cells and immune response cells, less the systemic cost of chemotherapy. We completely characterise the optimal control and compute a numerical solution of the optimality system via analytic continuation.Research supported by the Natural Science and Engineering Research Council (NSERC) and the Mathematics of Information Technology and Complex Systems (MITACS) of Canada     

Optimal choice between feedforward and feedback control in gene expression to cope with unpredictable danger

All organisms face risks of unpredictable danger caused by harmful physical environments, pathogens, parasites or predators. Organisms may have several alternative ways of coping with such dangers. These differ in cost, effectiveness and activation time. We study the conditions under which it is optimal to use different alternatives for damage control. As an example we consider a microbe (such as E. coli), which may experience heat shocks that cause denaturation of proteins in the cell. To restore the denatured proteins the organism produces heat-shock proteins (HSP). There are two different pathways for production of HSP. Some HSP are produced immediately after a heat shock (feedforward control), but additional HSP may be produced thereafter, stimulated by the presence of denatured proteins (feedback control). Feedforward is based solely on heat-shock intensity without accurate information on the resulting amount of denatured proteins. We examine the optimal combination of the two pathways that minimizes the sum of the damage caused by the presence of untreated denatured proteins and the production cost of HSP. The optimal response depends on the time delay for feedback control, the effectiveness of HSP in processing denatured proteins, the production cost of HSP, the severity of damage by denatured proteins and the probability distribution of the abundance of denatured protein conditional on heat-shock intensity. We find that feedforward control should always be used. Additional HSP may be produced by feedback control when the abundance of denatured protein is large whilst no feedback control should be used when it is small. All the HSP are produced by feedforward control when the maximum is close to the mean of denatured protein abundance conditional on the heat-shock intensity.     

Mathematical analysis of a model describing evolution of an asexual population in a changing environment

We investigate a mathematical model for an asexual population with non-overlapping (discrete) generations, that exists in a changing environment. Sexual populations are also briefly discussed at the end of the paper. It is assumed that selection occurs on the value of a single polygenic trait, which is controlled by a finite number of loci with discrete-effect alleles. The environmental change results in a moving fitness optimum, causing the trait to be subject to a combination of stabilising and directional selection.This model is different from that investigated by Waxman and Peck [Genetics 153 (1999) 1041] where overlapping generations and continuous effect alleles were considered. In this paper, we consider non-overlapping generations and discrete effect alleles. However in [Genetics 153 (1999) 1041] and the present work, there is the same pattern of environmental change, namely a constant rate of change of the optimum.From [Genetics 153 (1999) 1041], no rigorous theoretical conclusion can be drawn about the form of the solutions as t grows large. Numerical work carried out in [Genetics 153 (1999) 1041] suggests that the solution is a lagged travelling wave solution, but no mathematical proof exists for the continuous model. Only partial results, regarding existence of travelling wave solutions and perturbed solutions, have been established (see [Nonlin. Anal. 53 (2003) 683; An integral equation describing an asexual population in a changing environment, Preprint]). For the discrete case of this paper, under the assumption that the ratio between the unit of genotypic value and the speed of environment change is a rational number, we are able to give rigorous proof of the following conclusion: the population follows the environmental change with a small lag behind, moreover, the lag is represented using a calculable quantity.