Phylogenetic identification of lateral genetic transfer events

Background  

Lateral genetic transfer can lead to disagreements among phylogenetic trees comprising sequences from the same set of taxa. Where topological discordance is thought to have arisen through genetic transfer events, tree comparisons can be used to identify the lineages that may have shared genetic information. An 'edit path' of one or more transfer events can be represented with a series of subtree prune and regraft (SPR) operations, but finding the optimal such set of operations is NP-hard for comparisons between rooted trees, and may be so for unrooted trees as well.     

Implications of considering metal bioavailability in estimates of freshwater ecotoxicity: examination of two case studies

Purpose  

Previous methods of estimating characterization factors (CFs) of metals in life cycle impact assessment (LCIA) models were based on multimedia fate, exposure, and effect models originally developed to address the potential impacts of organic chemicals. When applied to metals, the models neglect the influence of ambient chemistry on metal speciation, bioavailability and toxicity. Gandhi et al. (2010) presented a new method of calculating CFs for freshwater ecotoxicity that addresses these metal-specific issues. In this paper, we compared and assessed the consequences of using the new method versus currently available LCIA models for calculating freshwater ecotoxicity, as applied to two case studies previously examined by Gloria et al. (2006): (1) the production of copper (Cu) pipe and (2) a zinc (Zn) gutter system.     

Methodology for systematic analysis and improvement of manufacturing unit process life cycle inventory (UPLCI) CO2PE! initiative (cooperative effort on process emissions in manufacturing). Part 2: case studies

Purpose  

This report presents two case studies, one for both the screening approach and the in-depth approach, demonstrating the application of the life cycle assessment-oriented methodology for systematic inventory analysis of the machine tool use phase of manufacturing unit processes, which has been developed in the framework of the CO₂PE! collaborative research programme (CO₂PE! 2011) and is described in part 1 of this paper (Kellens et al. 2011).     

A response to information criterion-based clustering with order-restricted candidate profiles in short time-course microarray experiments

Background  

For gene expression data obtained from a time-course microarray experiment, Liu et al. [1] developed a new algorithm for clustering genes with similar expression profiles over time. Performance of their proposal was compared with three other methods including the order-restricted inference based methodology of Peddada et al. [2, 3]. In this note we point out several inaccuracies in Liu et al. [1] and conclude that the order-restricted inference based methodology of Peddada et al. (programmed in the software ORIOGEN) indeed operates at the desired nominal Type 1 error level, an important feature of a statistical decision rule, while being computationally substantially faster than indicated by Liu et al. [1].     

Envelope structure of Synechococcus sp. WH8113, a nonflagellated swimming cyanobacterium

Background  

Many bacteria swim by rotating helical flagellar filaments [1]. Waterbury et al. [15] discovered an exception, strains of the cyanobacterium Synechococcus that swim without flagella or visible changes in shape. Other species of cyanobacteria glide on surfaces [2,7]. The hypothesis that Synechococcus might swim using traveling surface waves [6,13] prompted this investigation.     

Feature selection and classification for microarray data analysis: Evolutionary methods for identifying predictive genes

Background  

In the clinical context, samples assayed by microarray are often classified by cell line or tumour type and it is of interest to discover a set of genes that can be used as class predictors. The leukemia dataset of Golubet al.[1] and the NCI60 dataset of Rosset al.[2] present multiclass classification problems where three tumour types and nine cell lines respectively must be identified. We apply an evolutionary algorithm to identify the near-optimal set of predictive genes that classify the data. We also examine the initial gene selection step whereby the most informative genes are selected from the genes assayed.     

Life cycle assessment of Australian sugarcane products with a focus on cane processing

Purpose  

This work generates attributional life cycle assessment (LCA) results for products produced from Australian sugarcane—raw sugar, molasses, electricity (from bagasse combustion), and ethanol (from molasses). It focuses on cane processing in sugar mills and is a companion to the work presented in (Renouf et al. 2010), where the focus is on cane growing. This work also examines the preferred approach for assigning impacts to the multiple products from cane processing, and the influence that variability in cane growing has on the results.     

Use of weighted reference panels based on empirical estimates of ancestry for capturing untyped variation

Many association methods use a subset of genotyped single nucleotide polymorphisms (SNPs) to capture or infer genotypes at other untyped SNPs. We and others previously showed that tag SNPs selected to capture common variation using data from The International HapMap Consortium (Nature 437:1299–1320, 2005), The International HapMap Consortium (Nature 449:851–861, 2007) could also capture variation in populations of similar ancestry to HapMap reference populations (de Bakker et al. in Nat Genet 38:1298–1303, 2006; González-Neira et al. in Genome Res 16:323–330, 2006; Montpetit et al. in PLoS Genet 2:282–290, 2006; Mueller et al. in Am J Hum Genet 76:387–398, 2005). To capture variation in admixed populations or populations less similar to HapMap panels, a “cosmopolitan approach,” in which all samples from HapMap are used as a single reference panel, was proposed. Here we refine this suggestion and show that use of a “weighted reference panel,” constructed based on empirical estimates of ancestry in the target population (relative to available reference panels), is more efficient than the cosmopolitan approach. Weighted reference panels capture, on average, only slightly fewer common variants (minor allele frequency > 5%) than the cosmopolitan approach (mean r ² = 0.977 vs. 0.989, 94.5% variation captured vs. 96.8% at r ² > 0.8), across the five populations of the Multiethnic Cohort, but entail approximately 25% fewer tag SNPs per panel (average 538 vs. 718). These results extend a recent study in two Indian populations (Pemberton et al. in Ann Hum Genet 72:535–546, 2008). Weighted reference panels are potentially useful for both the selection of tag SNPs in diverse populations and perhaps in the design of reference panels for imputation of untyped genotypes in genome-wide association studies in admixed populations. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Genome-wide meta-analysis for rheumatoid arthritis

Meta-analysis is being increasingly used as a tool for integrating data from different studies of complex phenotypes, because the power of any one study to identify causal loci is limited. We applied a novel meta-analytical approach (Loesgen et al. in Genet Epidemiol 21(Suppl 1):S142–S147, 2001) in compiling results from four studies of rheumatoid arthritis in Caucasians including two studies from NARAC (Jawaheer et al. in Am J Hum Genet 68:927–936, 2001; Jawaheer et al. in Arthritis Rheum 48:906–916, 2003), one study from the UK (MacKay et al. in Arthritis Rheum 46:632–639, 2001) and one from France (Cornelis et al. in Proc Natl Acad Sci USA 95:10746–10750, 1998). For each study, we obtained NPL scores by performing interval mapping (2 cM intervals) using GeneHunter2 (Kruglyak et al. in Am J Hum Genet 58:1347–1363, 1996; Markianos et al. in Am J Hum Genet 68:963–977, 2001). The marker maps differed among the three consortium groups, therefore, the marker maps were aligned after the interval mapping was completed and the NPL scores that were within 1 cM of each other were combined using the method of Loesgen et al. (Genet Epidemiol 21(Suppl 1):S142–S147, 2001) by calculating the weighted average of the NPL score. This approach avoids some problems in analysis encountered by using GeneHunter2 when some markers in the sample are not genotyped. This procedure provided marginal evidence (P<0.05) of linkage on chromosome 1, 2, 5 and 18, strong evidence (P<0.01) on chromosomes 8 and 16, and overwhelming evidence in the HLA region of chromosome 6.     

Consistency of Topological Moves Based on the Balanced Minimum Evolution Principle of Phylogenetic Inference

Many phylogenetic algorithms search the space of possible trees using topological rearrangements and some optimality criterion. FastME is such an approach that uses the {em balanced minimum evolution (BME)} principle, which computer studies have demonstrated to have high accuracy. FastME includes two variants: {em balanced subtree prune and regraft (BSPR)} and {em balanced nearest neighbor interchange (BNNI)}. These algorithms take as input a distance matrix and a putative phylogenetic tree. The tree is modified using SPR or NNI operations, respectively, to reduce the BME length relative to the distance matrix, until a tree with (locally) shortest BME length is found. Following computer simulations, it has been conjectured that BSPR and BNNI are consistent, i.e. for an input distance that is a tree-metric, they converge to the corresponding tree. We prove that the BSPR algorithm is consistent. Moreover, even if the input contains small errors relative to a tree-metric, we show that the BSPR algorithm still returns the corresponding tree. Whether BNNI is consistent remains open.     

The human Ago2 MC region does not contain an eIF4E-like mRNA cap binding motif

Background  

Argonaute (Ago) proteins interact with small regulatory RNAs to mediate gene regulatory pathways. A recent report by Kiriakidou et al. [1] describes an MC sequence region identified in Ago2 that displays similarity to the cap-binding motif in translation initiation factor 4E (eIF4E). In a cap-bound eIF4E structure, two important aromatic residues of the motif stack on either side of a 7-methylguanosine 5'-triphosphate (m⁷Gppp) base. The corresponding Ago2 aromatic residues (F450 and F505) were hypothesized to perform the same cap-binding function. However, the detected similarity between the MC sequence and the eIF4E cap-binding motif was questionable.     

Rhythmic Oscillations in Quantitative EEG Measured During a Continuous Performance Task

The objective of the present investigation was to determine if cyclic variations in human performance recorded during a 30 min continuous performance task would parallel cyclic variations in right-hemisphere beta-wave activity. A fast fourier transformation was performed on the quantitative electroencephalogram (qEEG) and the performance record of each participant (N = 62), producing an individual periodogram for each outcome measure. An average periodogram was then produced for both qEEG and performance by combining (averaging) the amplitudes associated with each periodicity in the 62 original periodograms. Periodicities ranging from 1.00 to 2.00 min and from 4.70 to 5.70 min with amplitudes greater than would be expected due to chance were retained (Smith et al. 2003). The results of the present investigation validate the existence of cyclic variations in human performance that have been identified previously (Smith et al. 2003) and extend those findings by implicating right-hemisphere mediated arousal in the process (Arruda et al. 1996, 1999, 2007). Significant cyclic variations in left-hemisphere beta-wave activity were not observed. Taken together, the findings of the present investigation support a model of sustained attention that predicts cyclic changes in human performance that are the result of cyclic changes in right-hemisphere arousal.

Life cycle assessment of magnetic and electronic ballast for 36-W fluorescent lamp

Background, aim, and scope   

Ballast is a device in a fluorescent lamp that supports the production of light. In this study, the environmental impacts of two types of Malaysian ballast, magnetic ballast and electronic ballast, were identified and compared using the life cycle assessment approach through the ISO 14040 (2005) series.     

Environmental impacts of conventional plastic and bio-based carrier bags

Background, aim, and scope  

Worldwide, the production of biodegradable and compostable plastics has steadily grown. In Part 1 (Khoo et al. 2010), life cycle assessment (LCA) was applied to compare the production stages of a bio-based bag (made from polyhydroxyalkanoate or bio-plastic (PHA)) with polyethylene plastic bag. The scope of the study is within the context of Singapore and does not include other types of conventional or bio-based polymers (e.g., polylactic acid (PLA), thermoplastics, high-density polyethylene (HDPE), EPS, etc). This article (part 2) proposes to investigate the end-of-life options of both bags.     

Isometric gene tree reconciliation revisited

Background

Isometric gene tree reconciliation is a gene tree/species tree reconciliation problem where both the gene tree and the species tree include branch lengths, and these branch lengths must be respected by the reconciliation. The problem was introduced by Ma et al. in 2008 in the context of reconstructing evolutionary histories of genomes in the infinite sites model.

Results

In this paper, we show that the original algorithm by Ma et al. is incorrect, and we propose a modified algorithm that addresses the problems that we discovered. We have also improved the running time from \(O(N^2)\) to \(O(N\log N)\), where N is the total number of nodes in the two input trees. Finally, we examine two new variants of the problem: reconciliation of two unrooted trees and scaling of branch lengths of the gene tree during reconciliation of two rooted trees.

Conclusions

We provide several new algorithms for isometric reconciliation of trees. Some questions in this area remain open; most importantly extensions of the problem allowing for imprecise estimates of branch lengths.

     

A Computational Analysis of Localized Ca<Superscript>2+</Superscript>-Dynamics Generated by Heterogeneous Release Sites

We investigate the role of heterogeneous expression of IP₃R and RyR in generating diverse elementary Ca²⁺ signals. It has been shown empirically (Wojcikiewicz and Luo in Mol. Pharmacol. 53(4):656–662, 1998; Newton et al. in J. Biol. Chem. 269(46):28613–28619, 1994; Smedt et al. in Biochem. J. 322(Pt. 2):575–583, 1997) that tissues express various proportions of IP₃ and RyR isoforms and this expression is dynamically regulated (Parrington et al. in Dev. Biol. 203(2):451–461, 1998; Fissore et al. in Biol. Reprod. 60(1):49–57, 1999; Tovey et al. in J. Cell Sci. 114(Pt. 22):3979–3989, 2001). Although many previous theoretical studies have investigated the dynamics of localized calcium release sites (Swillens et al. in Proc. Natl. Acad. Sci. U.S.A. 96(24):13750–13755, 1999; Shuai and Jung in Proc. Natl. Acad. Sci. U.S.A. 100(2):506–510, 2003a; Shuai and Jung in Phys. Rev. E, Stat. Nonlinear Soft Matter Phys. 67(3 Pt. 1):031905, 2003b; Thul and Falcke in Biophys. J. 86(5):2660–2673, 2004; DeRemigio and Smith in Cell Calcium 38(2):73–86, 2005; Nguyen et al. in Bull. Math. Biol. 67(3):393–432, 2005), so far all such studies focused on release sites consisting of identical channel types. We have extended an existing mathematical model (Nguyen et al. in Bull. Math. Biol. 67(3):393–432, 2005) to release sites with two (or more) receptor types, each with its distinct channel kinetics. Mathematically, the release site is represented by a transition probability matrix for a collection of nonidentical stochastically gating channels coupled through a shared Ca²⁺ domain. We demonstrate that under certain conditions a previously defined mean-field approximation of the coupling strength does not accurately reproduce the release site dynamics. We develop a novel approximation and establish that its performance in these instances is superior. We use this mathematical framework to study the effect of heterogeneity in the Ca²⁺-regulation of two colocalized channel types on the release site dynamics. We consider release sites consisting of channels with both Ca²⁺-activation and inactivation (“four-state channels”) and channels with Ca²⁺-activation only (“two-state channels”) and show that for the appropriate parameter values, synchronous channel openings within a release site with any proportion of two-state to four-state channels are possible, however, the larger the proportion of two-state channels, the more sensitive the dynamics are to the exact spatial positioning of the channels and the distance between channels. Specifically, the clustering of even a small number of two-state channels interferes with puff/spark termination and increases puff durations or leads to a tonic response.     

Flexible taxonomic assignment of ambiguous sequencing reads

Background  

To characterize the diversity of bacterial populations in metagenomic studies, sequencing reads need to be accurately assigned to taxonomic units in a given reference taxonomy. Reads that cannot be reliably assigned to a unique leaf in the taxonomy (ambiguous reads) are typically assigned to the lowest common ancestor of the set of species that match it. This introduces a potentially severe error in the estimation of bacteria present in the sample due to false positives, since all species in the subtree rooted at the ancestor are implicitly assigned to the read even though many of them may not match it.     

Phylogeny of Prokaryotes and Chloroplasts Revealed by a Simple Composition Approach on All Protein Sequences from Complete Genomes Without Sequence Alignment

The complete genomes of living organisms have provided much information on their phylogenetic relationships. Similarly, the complete genomes of chloroplasts have helped to resolve the evolution of this organelle in photosynthetic eukaryotes. In this paper we propose an alternative method of phylogenetic analysis using compositional statistics for all protein sequences from complete genomes. This new method is conceptually simpler than and computationally as fast as the one proposed by Qi et al. (2004b) and Chu et al. (2004). The same data sets used in Qi et al. (2004b) and Chu et al. (2004) are analyzed using the new method. Our distance-based phylogenic tree of the 109 prokaryotes and eukaryotes agrees with the biologists tree of life based on 16S rRNA comparison in a predominant majority of basic branching and most lower taxa. Our phylogenetic analysis also shows that the chloroplast genomes are separated to two major clades corresponding to chlorophytes s.l. and rhodophytes s.l. The interrelationships among the chloroplasts are largely in agreement with the current understanding on chloroplast evolution.Reviewing Editor: Dr. John Oakeshott     

Basic properties and information theory of Audic-Claverie statistic for analyzing cDNA arrays

Background  

The Audic-Claverie method [1] has been and still continues to be a popular approach for detection of differentially expressed genes in the SAGE framework. The method is based on the assumption that under the null hypothesis tag counts of the same gene in two libraries come from the same but unknown Poisson distribution. The problem is that each SAGE library represents only a single measurement. We ask: Given that the tag count samples from SAGE libraries are extremely limited, how useful actually is the Audic-Claverie methodology? We rigorously analyze the A-C statistic that forms a backbone of the methodology and represents our knowledge of the underlying tag generating process based on one observation.     

The synaptonemal complex and meiotic recombination in humans: new approaches to old questions

Meiotic prophase serves as an arena for the interplay of two important cellular activities, meiotic recombination and synapsis of homologous chromosomes. Synapsis is mediated by the synaptonemal complex (SC), originally characterized as a structure linked to pairing of meiotic chromosomes (Moses (1958) J Biophys Biochem Cytol 4:633–638). In 1975, the first electron micrographs of human pachytene stage SCs were presented (Moses et al. (1975) Science 187:363–365) and over the next 15 years the importance of the SC to normal meiotic progression in human males and females was established (Jhanwar and Chaganti (1980) Hum Genet 54:405–408; Pathak and Elder (1980) Hum Genet 54:171–175; Solari (1980) Chromosoma 81:315–337; Speed (1984) Hum Genet 66:176–180; Wallace and Hulten (1985) Ann Hum Genet 49(Pt 3):215–226). Further, these studies made it clear that abnormalities in the assembly or maintenance of the SC were an important contributor to human infertility (Chaganti et al. (1980) Am J Hum Genet 32:833–848; Vidal et al. (1982) Hum Genet 60:301–304; Bojko (1983) Carlsberg Res Commun 48:285–305; Bojko (1985) Carlsberg Res Commun 50:43–72; Templado et al. (1984) Hum Genet 67:162–165; Navarro et al. (1986) Hum Reprod 1:523–527; Garcia et al. (1989) Hum Genet 2:147–53). However, the utility of these early studies was limited by lack of information on the structural composition of the SC and the identity of other SC-associated proteins. Fortunately, studies of the past 15 years have gone a long way toward remedying this problem. In this minireview, we highlight the most important of these advances as they pertain to human meiosis, focusing on temporal aspects of SC assembly, the relationship between the SC and meiotic recombination, and the contribution of SC abnormalities to human infertility.The synaptonemal complex–50 years