Robust inference of positive selection from recombining coding sequences

MOTIVATION: Accurate detection of positive Darwinian selection can provide important insights to researchers investigating the evolution of pathogens. However, many pathogens (particularly viruses) undergo frequent recombination and the phylogenetic methods commonly applied to detect positive selection have been shown to give misleading results when applied to recombining sequences. We propose a method that makes maximum likelihood inference of positive selection robust to the presence of recombination. This is achieved by allowing tree topologies and branch lengths to change across detected recombination breakpoints. Further improvements are obtained by allowing synonymous substitution rates to vary across sites. RESULTS: Using simulation we show that, even for extreme cases where recombination causes standard methods to reach false positive rates >90%, the proposed method decreases the false positive rate to acceptable levels while retaining high power. We applied the method to two HIV-1 datasets for which we have previously found that inference of positive selection is invalid owing to high rates of recombination. In one of these (env gene) we still detected positive selection using the proposed method, while in the other (gag gene) we found no significant evidence of positive selection. AVAILABILITY: A HyPhy batch language implementation of the proposed methods and the HIV-1 datasets analysed are available at http://www.cbio.uct.ac.za/pub_support/bioinf06. The HyPhy package is available at http://www.hyphy.org, and it is planned that the proposed methods will be included in the next distribution. RDP2 is available at http://darwin.uvigo.es/rdp/rdp.html     

CodonTest: Modeling Amino Acid Substitution Preferences in Coding Sequences

Codon models of evolution have facilitated the interpretation of selective forces operating on genomes. These models, however, assume a single rate of non-synonymous substitution irrespective of the nature of amino acids being exchanged. Recent developments have shown that models which allow for amino acid pairs to have independent rates of substitution offer improved fit over single rate models. However, these approaches have been limited by the necessity for large alignments in their estimation. An alternative approach is to assume that substitution rates between amino acid pairs can be subdivided into rate classes, dependent on the information content of the alignment. However, given the combinatorially large number of such models, an efficient model search strategy is needed. Here we develop a Genetic Algorithm (GA) method for the estimation of such models. A GA is used to assign amino acid substitution pairs to a series of rate classes, where is estimated from the alignment. Other parameters of the phylogenetic Markov model, including substitution rates, character frequencies and branch lengths are estimated using standard maximum likelihood optimization procedures. We apply the GA to empirical alignments and show improved model fit over existing models of codon evolution. Our results suggest that current models are poor approximations of protein evolution and thus gene and organism specific multi-rate models that incorporate amino acid substitution biases are preferred. We further anticipate that the clustering of amino acid substitution rates into classes will be biologically informative, such that genes with similar functions exhibit similar clustering, and hence this clustering will be useful for the evolutionary fingerprinting of genes.     

Decreased External Skeletal Robustness in Schoolchildren – A Global Trend? Ten Year Comparison of Russian and German Data

Objectives

Obesity and a reduced physical activity are global developments. Physical activity affects the external skeletal robustness which decreased in German children. It was assumed that the negative trend of decreased external skeletal robustness can be found in other countries. Therefore anthropometric data of Russian and German children from the years 2000 and 2010 were compared.

Methods

Russian (2000/2010 n = 1023/268) and German (2000/2010 n = 2103/1750) children aged 6–10 years were investigated. Height, BMI and external skeletal robustness (Frame-Index) were examined and compared for the years and the countries. Statistical analysis was performed by Mann-Whitney-Test.

Results

Comparison 2010 and 2000: In Russian children BMI was significantly higher; boys were significantly taller and exhibited a decreased Frame-Index (p = .002) in 2010. German boys showed significantly higher BMI in 2010. In both sexes Frame-Index (p = .001) was reduced in 2010. Comparison Russian and German children in 2000: BMI, height and Frame-Index were different between Russian and German children. German children were significantly taller but exhibited a lower Frame-Index (p<.001). Even German girls showed a significantly higher BMI. Comparison Russian and German children in 2010: BMI and Frame-Index were different. Russian children displayed a higher Frame-Index (p<.001) compared with Germans.

Conclusions

In Russian children BMI has increased in recent years. Frame-Index is still higher in Russian children compared with Germans however in Russian boys Frame-Index is reduced. This trend and the physical activity should be observed in the future.     

Derivation and large-scale expansion of multipotent astroglial neural progenitors from adult human brain

The isolation and expansion of human neural cell types has become increasingly relevant in restorative neurobiology. Although embryonic and fetal tissue are frequently envisaged as providing sufficiently primordial cells for such applications, the developmental plasticity of endogenous adult neural cells remains largely unclear. To examine the developmental potential of adult human brain cells, we applied conditions favoring the growth of neural stem cells to multiple cortical regions, resulting in the identification and selection of a population of adult human neural progenitors (AHNPs). These nestin(+) progenitors may be derived from multiple forebrain regions, are maintainable in adherent conditions, co-express multiple glial and immature markers, and are highly expandable, allowing a single progenitor to theoretically form sufficient cells for approximately 4x10(7) adult brains. AHNPs longitudinally maintain the ability to generate both glial and neuronal cell types in vivo and in vitro, and are amenable to genetic modification and transplantation. These findings suggest an unprecedented degree of inducible plasticity is retained by cells of the adult central nervous system.     

瓜类枯萎病研究进展(综述)

瓜类枯萎病是由尖孢镰刀菌Fusarium oxysporum侵染引起,在全球大部分瓜类产区均有发生。该病已成为影响我国瓜类品质和产量的重要病害。本文从瓜类作物枯萎病的防治方法、致病机理和分子生物学研究等方面综述其研究现状,并对今后的研究方向进行展望。     

Detecting phylogenetic signal in mutualistic interaction networks using a Markov process model

Ecological interaction networks, such as those describing the mutualistic interactions between plants and their pollinators or between plants and their frugivores, exhibit non‐random structural properties that cannot be explained by simple models of network formation. One factor affecting the formation and eventual structure of such a network is its evolutionary history. We argue that this, in many cases, is closely linked to the evolutionary histories of the species involved in the interactions. Indeed, empirical studies of interaction networks along with the phylogenies of the interacting species have demonstrated significant associations between phylogeny and network structure. To date, however, no generative model explaining the way in which the evolution of individual species affects the evolution of interaction networks has been proposed. We present a model describing the evolution of pairwise interactions as a branching Markov process, drawing on phylogenetic models of molecular evolution. Using knowledge of the phylogenies of the interacting species, our model yielded a significantly better fit to 21% of a set of plant–pollinator and plant–frugivore mutualistic networks. This highlights the importance, in a substantial minority of cases, of inheritance of interaction patterns without excluding the potential role of ecological novelties in forming the current network architecture. We suggest that our model can be used as a null model for controlling evolutionary signals when evaluating the role of other factors in shaping the emergence of ecological networks.     

Molecular genetics of complex I-deficient Chinese hamster cell lines

The work from our laboratory on complex I-deficient Chinese hamster cell mutants is reviewed. Several complementation groups with a complete defect have been identified. Three of these are due to X-linked mutations, and the mutated genes for two have been identified. We describe null mutants in the genes for the subunits MWFE (gene: NDUFA1) and ESSS. They represent small integral membrane proteins localized in the Ialpha (Igamma) and Ibeta subcomplexes, respectively [J. Hirst, J. Carroll, I.M. Fearnley, R.J. Shannon, J.E. Walker. The nuclear encoded subunits of complex I from bovine heart mitochondria. Biochim. Biophys. Acta 1604 (7-10-2003) 135-150.]. Both are absolutely essential for assembly and activity of complex I. Epitope-tagged versions of these proteins can be expressed from a poly-cistronic vector to complement the mutants, or to be co-expressed with the endogenous proteins in other hamster cell lines (mutant or wild type), or human cells. Structure-function analyses can be performed with proteins altered by site-directed mutagenesis. A cell line has been constructed in which the MWFE subunit is conditionally expressed, opening a window on the kinetics of assembly of complex I. Its targeting, import into mitochondria, and orientation in the inner membrane have also been investigated. The two proteins have recently been shown to be the targets for a cAMP-dependent kinase [R. Chen, I.M. Fearnley, S.Y. Peak_Chew, J.E. Walker. The phosphorylation of subunits of complex I from bovine heart mitochondria. J. Biol. Chem. xx (2004) xx-xx.]. The epitope-tagged proteins can be cross-linked with other complex I subunits.     

Functional network integration of embryonic stem cell-derived astrocytes in hippocampal slice cultures

Embryonic stem (ES) cells provide attractive prospects for neural transplantation. So far, grafting strategies in the CNS have focused mainly on neuronal replacement. Employing a slice culture model, we found that ES cell-derived glial precursors (ESGPs) possess a remarkable capacity to integrate into the host glial network. Following deposition on the surface of hippocampal slices, ESGPs actively migrate into the recipient tissue and establish extensive cell-cell contacts with recipient glia. Gap junction-mediated coupling between donor and host astrocytes permits widespread delivery of dye from single donor cells. During maturation, engrafted donor cells display morphological, immunochemical and electrophysiological properties that are characteristic of differentiating native glia. Our findings provide the first evidence of functional integration of grafted astrocytes, and depict glial network integration as a potential route for widespread transcellular delivery of small molecules to the CNS.     

The tumor suppressor cybL,a component of the respiratory chain,mediates apoptosis induction

A genetic screen was established to clone apoptosis-inducing genes in a high-throughput format. It led to the isolation of several proapoptotic genes whose proteins are localized to mitochondria. One of the isolated genes is cytochrome bL (cybL also known as SDHC, CII-3, or QPs-1), a component of the respiratory chain complex II. It was further investigated because both cybL and another component of complex II, cybS, have recently been identified as tumor suppressor proteins, some of which act by controlling apoptosis. Our studies reveal that cell death induction by cybL expression is concomitant with a transient inhibition of complex II and the generation of reactive oxygen species. Importantly, cells that are constitutively deficient in cybL are resistant to a variety of proapoptotic cytostatic drugs and to the effects of the Fas receptor. Our results therefore identify complex II as a sensor for apoptosis induction and could explain the unexpected observation that complex II is inactivated in tumors.     

Species-specific and mutant MWFE proteins. Their effect on the assembly of a functional mammalian mitochondrial complex I

The MWFE protein (70 amino acids) is highly conserved in evolution, but the human protein (80% identical to hamster) does not complement a null mutation in Chinese hamster cells. We have identified a small protein segment where significant differences exist between rodents and primates, illustrating very specifically the need for compatibility of the nuclear and mitochondrial genomes in the assembly of complex I. The segment between amino acids 39 and 46 appears to be critical for species-specific compatibility. Amino acid substitutions in this region were tested that caused a reduction of activity of the hamster protein or converted the inactive human protein into a partially active one. Such mutations could be useful in making mice with partial complex I activity as models for mitochondrial diseases. Their potential as dominant negative mutants was explored. More deleterious mutations in the NDUFA1 gene were also characterized. A conservative substitution, R50K, or a short C-terminal deletion makes the protein completely inactive. In the absence of MWFE, no high molecular weight complex was detectable by Blue Native-gel electrophoresis. The MWFE protein itself is unstable in the absence of assembled mitochondrially encoded integral membrane proteins of complex I.