Speeding Up Ecological and Evolutionary Computations in R; Essentials of High Performance Computing for Biologists

Computation has become a critical component of research in biology. A risk has emerged that computational and programming challenges may limit research scope, depth, and quality. We review various solutions to common computational efficiency problems in ecological and evolutionary research. Our review pulls together material that is currently scattered across many sources and emphasizes those techniques that are especially effective for typical ecological and environmental problems. We demonstrate how straightforward it can be to write efficient code and implement techniques such as profiling or parallel computing. We supply a newly developed R package (aprof) that helps to identify computational bottlenecks in R code and determine whether optimization can be effective. Our review is complemented by a practical set of examples and detailed Supporting Information material (S1–S3 Texts) that demonstrate large improvements in computational speed (ranging from 10.5 times to 14,000 times faster). By improving computational efficiency, biologists can feasibly solve more complex tasks, ask more ambitious questions, and include more sophisticated analyses in their research.

This is part of the PLOS Computational Biology Education collection.

     

Evidence for a highly asymmetric arrangement of ether- and diacyl-phospholipid subclasses in the plasma membrane of Krebs II ascites cells

(1) Krebs II ascites cells were taken as a model of the neoplastic cells to investigate the transverse distribution of phospholipids in the plasma membrane. The experimental procedure was based on non-lytic degradation of phospholipids in the intact cell by Naja naja phospholipase A₂ and Staphylococcus aureus sphingomyelinase C and on phopholipid analysis of purified plasma membranes. It was shown that the three major phospholipids, i.e., phosphatidylcholine, phosphatidylethanolamine and sphingomyelin, are randomly distributed between the two halves of the membranes, whereas phosphatidylserine remains located in the inner leaflet. (2) The membrane localization of phosphatidylcholine and phosphatidylethanolamine subclasses (diacyl, alkylacyl and alkenylacyl) was also examined, using a new procedure of ether-phospholipid determination. The method involves a selective removal of diacyl species by guinea pig pancreas phospholipase A₁ and of alkenylacyl species by acidolysis. This analysis revealed a 50% increase of ether phospholipids in the plasma membrane as compared to the whole cell (36.5 and 23.1% of total phospholipid, respectively). Furthermore, a strong membrane asymmetry was demonstrated for the three phosphatidylcholine subclasses, since 1-alkyl-2-acyl-sn-glycerol-3-phosphocholine (alkylacyl-GPC) was entirely found in the inner leaflet, whereas both diacyl- and alkenylacyl-GPC displayed an external localization. The same pattern was observed for phosphatidylethanolamine subclasses, except for 1-alkenyl-2-acyl-sn-glycero-3-phosphoethanolamine, which was found randomly distributed. These results are discussed in relation to the process of cell malignant transformation and to the biosynthesis of platelet-activating factor (PAF-acether or 1-alkyl-2-acetyl-GPC).     

Construction of Engineered Bifunctional Enzymes and Their Overproduction in Aspergillus niger for Improved Enzymatic Tools To Degrade Agricultural By-Products

Two chimeric enzymes, FLX and FLXLC, were designed and successfully overproduced in Aspergillus niger. FLX construct is composed of the sequences encoding the feruloyl esterase A (FAEA) fused to the endoxylanase B (XYNB) of A. niger. A C-terminal carbohydrate-binding module (CBM family 1) was grafted to FLX, generating the second hybrid enzyme, FLXLC. Between each partner, a hyperglycosylated linker was included to stabilize the constructs. Hybrid proteins were purified to homogeneity, and molecular masses were estimated to be 72 and 97 kDa for FLX and FLXLC, respectively. Integrity of hybrid enzymes was checked by immunodetection that showed a single form by using antibodies raised against FAEA and polyhistidine tag. Physicochemical properties of each catalytic module of the bifunctional enzymes corresponded to those of the free enzymes. In addition, we verified that FLXLC exhibited an affinity for microcrystalline cellulose (Avicel) with binding parameters corresponding to a K_d of 9.9 × 10⁻⁸ M for the dissociation constant and 0.98 μmol/g Avicel for the binding capacity. Both bifunctional enzymes were investigated for their capacity to release ferulic acid from natural substrates: corn and wheat brans. Compared to free enzymes FAEA and XYNB, a higher synergistic effect was obtained by using FLX and FLXLC for both substrates. Moreover, the release of ferulic acid from corn bran was increased by using FLXLC rather than FLX. This result confirms a positive role of the CBM. In conclusion, these results demonstrated that the fusion of naturally free cell wall hydrolases and an A. niger-derived CBM onto bifunctional enzymes enables the increase of the synergistic effect on the degradation of complex substrates.     

Production of a chimeric enzyme tool associating the Trichoderma reesei swollenin with the Aspergillus niger feruloyl esterase A for release of ferulic acid

The main goals of this work were to produce the fusion protein of the Trichoderma reesei swollenin I (SWOI) and Aspergillus niger feruloyl esterase A (FAEA) and to study the effect of the physical association of the fusion partners on the efficiency of the enzyme. The fusion protein was produced up to 25 mg l⁻¹ in the T. reesei strains Rut-C30 and CL847. In parallel, FAEA alone was produced for use as a control protein in application tests. Recombinant FAEA and SWOI–FAEA were purified to homogeneity and characterized. The biochemical and kinetic characteristics of the two recombinant proteins were found to be similar to those of native FAEA, except for the temperature stability and specific activity of the SWOI–FAEA. Finally, the SWOI–FAEA protein was tested for release of ferulic acid from wheat bran. A period of 24 h of enzymatic hydrolysis with the SWOI–FAEA improved the efficiency of ferulic acid release by 50% compared with the results obtained using the free FAEA and SWOI. Ferulic acid is used as an antioxidant and flavor precursor in the food and pharmaceutical industries. This is the first report of a potential application of the SWOI protein fused with an enzyme of industrial interest.     

Black tea,green tea,and tea polyphenols

Previous studies have shown that tea consumption can impair trace element metabolism, particularly iron status, and increase the risk of anemia in humans and animals. More recently, however, evidence has been accumulating to show that, in animals, consumption of green tea or its polyphenols is associated with a reduction of the incidence and severity of a variety of experimentally induced cancers. In this study we have monitored the growth, trace element status, including hematological parameters of weanling rats given either (1) water, (2) 1% black tea, (3) 1% green, tea, or (4) 0.2% crude green tea extract as their sole drinking fluid while consuming diets containing either adequate or low amounts of iron. With the exception of manaanese, none of the trace elements studied (iron, copper, zinc, and manganese) or the hematological indices measured were affected by the type of beverage supplied, even though the polyphenol extract was shown to chelate metals in vitro and all the animals fed the low iron diet were shown to be anemic. There appeared to be an effect of black and green teas on manganese balance in, both the first and last weeks of the study. A lower level of brain managanese was associated with green tea consumption, and a higher level of this element in the kidneys of animals fed black tea. The results demonstrate that both black and green teas and a green tea polyphenol extract do not represent a risk to animals consuming the beverages as their sole fluid intake with respect to iron availability, although the interactions with manganese deserve further study.     

Fractionation of subcellular membranes of the secretory pathway from the peroxidase-producing white-rot fungus Phanerochaete chrysosporium

Abstract Mycelia from the basidiomycete Phanerochaete chrysosporium , producing lignin and manganese peroxidases, were homogenized and fractionated on a sucrose gradient. The main subcellular fungal membrane fractions were successfully separated. Lipid composition analyses of the isolated membranes as well as associated marker enzymes distribution gave evidence to similarities with membranes originating from plants. Lignin and manganese peroxidases were investigated by immunodetection in subcellular fractions. Our results show that lignin and manganese peroxidases are mainly associated with Golgi apparatus vesicles and, to a lesser extent, with endoplasmic reticulum and light density vesicles, but not with plasma membranes.     

Candida Endocarditis Treated with 5-Fluorocytosine

Three patients with endocarditis due to Candida species were treated with 5-fluorocytosine. One had endocarditis imposed on senile aortic sclerosis and survived and two had infections on prosthetic valves and died. In the latter two patients chemotherapy failed to eradicate the infection and was associated with marrow depression and possibly hepatic necrosis.     

Insights into Exo- and Endoglucanase Activities of Family 6 Glycoside Hydrolases from Podospora anserina

The ascomycete Podospora anserina is a coprophilous fungus that grows at late stages on droppings of herbivores. Its genome encodes a large diversity of carbohydrate-active enzymes. Among them, four genes encode glycoside hydrolases from family 6 (GH6), the members of which comprise putative endoglucanases and exoglucanases, some of them exerting important functions for biomass degradation in fungi. Therefore, this family was selected for functional analysis. Three of the enzymes, P. anserina Cel6A (PaCel6A), PaCel6B, and PaCel6C, were functionally expressed in the yeast Pichia pastoris. All three GH6 enzymes hydrolyzed crystalline and amorphous cellulose but were inactive on hydroxyethyl cellulose, mannan, galactomannan, xyloglucan, arabinoxylan, arabinan, xylan, and pectin. PaCel6A had a catalytic efficiency on cellotetraose comparable to that of Trichoderma reesei Cel6A (TrCel6A), but PaCel6B and PaCel6C were clearly less efficient. PaCel6A was the enzyme with the highest stability at 45°C, while PaCel6C was the least stable enzyme, losing more than 50% of its activity after incubation at temperatures above 30°C for 24 h. In contrast to TrCel6A, all three studied P. anserina GH6 cellulases were stable over a wide range of pHs and conserved high activity at pH values of up to 9. Each enzyme displayed a distinct substrate and product profile, highlighting different modes of action, with PaCel6A being the enzyme most similar to TrCel6A. PaCel6B was the only enzyme with higher specific activity on carboxymethylcellulose (CMC) than on Avicel and showed lower processivity than the others. Structural modeling predicts an open catalytic cleft, suggesting that PaCel6B is an endoglucanase.     

Highly efficient production of laccase by the basidiomycete Pycnoporus cinnabarinus

An efficient transformation and expression system was developed for the industrially relevant basidiomycete Pycnoporus cinnabarinus. This was used to transform a laccase-deficient monokaryotic strain with the homologous lac1 laccase gene placed under the regulation of its own promoter or that of the SC3 hydrophobin gene or the glyceraldehyde-3-phosphate dehydrogenase (GPD) gene of Schizophyllum commune. SC3-driven expression resulted in a maximal laccase activity of 107 nkat ml(-1) in liquid shaken cultures. This value was about 1.4 and 1.6 times higher in the cases of the GPD and lac1 promoters, respectively. lac1-driven expression strongly increased when 25 g of ethanol liter(-1) was added to the medium. Accordingly, laccase activity increased to 1,223 nkat ml(-1). These findings agree with the fact that ethanol induces laccase gene expression in some fungi. Remarkably, lac1 mRNA accumulation and laccase activity also strongly increased in the presence of 25 g of ethanol liter(-1) when lac1 was expressed behind the SC3 or GPD promoter. In the latter case, a maximal laccase activity of 1,393 nkat ml(-1) (i.e., 360 mg liter(-1)) was obtained. Laccase production was further increased in transformants expressing lac1 behind its own promoter or that of GPD by growth in the presence of 40 g of ethanol liter(-1). In this case, maximal activities were 3,900 and 4,660 nkat ml(-1), respectively, corresponding to 1 and 1.2 g of laccase per liter and thus representing the highest laccase activities reported for recombinant fungal strains. These results suggest that P. cinnabarinus may be a host of choice for the production of other proteins as well.     

Roles of cytoplasmic osmolytes, water, and crowding in the response of Escherichia coli to osmotic stress: biophysical basis of osmoprotection by glycine betaine

To better understand the biophysical basis of osmoprotection by glycine betaine (GB) and the roles of cytoplasmic osmolytes, water, and macromolecular crowding in the growth of osmotically stressed Escherichia coli, we have determined growth rates and amounts of GB, K(+), trehalose, biopolymers, and water in the cytoplasm of E. coli K-12 grown over a wide range of high external osmolalities (1.02-2.17 Osm) in MOPS-buffered minimal medium (MBM) containing 1 mM betaine (MBM+GB). As osmolality increases, we observe that the amount of cytoplasmic GB increases, the amounts of K(+) (the other major cytoplasmic solute) and of biopolymers remain relatively constant, and the growth rate and the amount of cytoplasmic water decrease strongly, so concentrations of biopolymers and all solutes increase with increasing osmolality. We observe the same correlation between the growth rate and the amount of cytoplasmic water for cells grown in MBM+GB as in MBM, supporting our proposal that the amount of cytoplasmic water is a primary determinant of the growth rate of osmotically stressed cells. We also observe the same correlation between cytoplasmic concentrations of biopolymers and K(+) for cells grown in MBM and MBM+GB, consistent with our hypothesis of compensation between the anticipated large perturbing effects on cytoplasmic protein-DNA interactions of increases in cytoplasmic concentrations of K(+) and biopolymers (crowding) with increasing osmolality. For growth conditions where the amount of cytoplasmic water is relatively large, we find that cytoplasmic osmolality is adequately predicted by assuming that contributions of individual solutes to osmolality are additive and using in vitro osmotic data on osmolytes and a local bulk domain model for cytoplasmic water. At moderate growth osmolalities (up to 1 Osm), we conclude that GB is an efficient osmoprotectant because it is almost as excluded from the biopolymer surface in the cytoplasm as it is from native protein surface in vitro. At very high growth osmolalities where cells contain little cytoplasmic water, predicted cytoplasmic osmolalities greatly exceed observed osmolalities, and the efficiency of GB as an osmolality booster decreases as the amount of cytoplasmic water decreases.