Intracellular reduction of selenite into glutathione peroxidase. Evidence for involvement of NADPH and not glutathione as the reductant

Selenium (Se) in selenite is present in an oxidized state, and must be reduced for it to be incorporated as selenocysteine into selenoenzymes such as glutathione peroxidase (GPx). In vitro, Se, as in selenite, can be reduced utilizing glutathione (GSH) and glutathione reductase (GRed). We determined the effects of decreasing GSH levels, inhibiting GRed activity, and decreasing cellular NADPH on the selenite-dependent rate of GPx synthesis in cultured cells: PC3, CHO, and the E89 glucose-6-phosphate dehydrogenase (G-6-PD)-deficient cell line. A novel statistical analysis method was developed (using Box Cox transformed regression and a bootstrap method) in order to assess the effects of these manipulations singly and in combinations. Buthionine sulfoximine (BSO) was used to decrease GSH levels, 1,3 bis-(2 chloroethyl)-1-nitrosourea (BCNU) was used to inhibit GRed activity and methylene blue (MB) was used to decrease cellular NADPH levels. This statistical method evaluates the effects of BSO, BCNU, MB and selenite alone and in combinations on GPx activity. Decreasing the GSH level (< 5% of control) did not have an effect on the selenite-dependent rate of GPx synthesis in PC3 or CHO cells, but did have a small inhibitory effect on the rate of GPx synthesis in E89 cells. Inhibiting GRed activity was also associated with either no effect (CHO, E89) or a small effect (PC3) on GPx activity. In contrast, decreasing NADPH levels in cells treated with MB was associated with a large decrease in the selenite-dependent rate of GPx synthesis to 36, 34 and 25% of control in PC3, CHO, and E89 cells, respectively. The effects of BSO plus BCNU were not synergistic in any of the cell lines. The effects of BSO plus MB were synergistic in G-6-PD-deficient E89 cells, but not in PC3 or CHO cells. We therefore conclude that under normal culture conditions, NADPH, and not glutathione, is the primary reductant of Se in selenite to forms that are eventually incorporated into GPx. For cells with abnormal ability to generate NADPH, lowering the GSH levels had a small effect on selenite-dependent GPx synthesis. GRed activity is not required for the selenite-dependent synthesis of GPx.     

Rectangular tolerance regions and multivariate normal reference regions in laboratory medicine

Reference intervals are widely used in the interpretation of results of biochemical and physiological tests of patients. When there are multiple biochemical analytes measured from each subject, a multivariate reference region is needed. Because of their greater specificity against false positives, such reference regions are more desirable than separate univariate reference intervals that disregard the cross-correlations between variables. Traditionally, under multivariate normality, reference regions have been constructed as ellipsoidal regions. This approach suffers from a major drawback: it cannot detect component-wise extreme observations. In the present work, procedures are developed to construct rectangular reference regions in the multivariate normal setup. The construction is based on the criteria for tolerance intervals. The problems addressed include the computation of a rectangular tolerance region and simultaneous tolerance intervals. Also addressed is the computation of mixed reference intervals that include both two-sided and one-sided limits, simultaneously. A parametric bootstrap approach is used in the computations, and the accuracy of the proposed methodology is assessed using estimated coverage probabilities. The problem of sample size determination is also addressed, and the results are illustrated using examples that call for the computation of reference regions.     

Estimating index of population trend by re-sampling techniques (jackknife and bootstrap) and its application to the life table study of rice leaf roller, Cnaphalocrocis medinalis (Lepidoptera: Pyralidae)

Taking a published natural population life table of rice leaf roller, Cnaphalocrocis medinalis (Lepidoptera: Pyralidae), as an example, we estimated the population trend index, I, via re‐sampling methods (jackknife and bootstrap), determined its statistical properties and illustrated the application of these methods in determining the control effectiveness of bio‐agents and chemical insecticides. Depending on the simulation outputs, the smoothed distribution pattern of the estimates of I by delete‐1 jackknife is visually distinguishable from the normal density, but the smoothed pattern produced by delete‐d jackknife, and logarithm‐transformed smoothed patterns produced by both empirical and parametric bootstraps, matched well the corresponding normal density. Thus, the estimates of I produced by delete‐1 jackknife were not used to determine the suppressive effect of wasps and insecticides. The 95% percent confidence intervals or the narrowest 95 percentiles and Z‐test criterion were employed to compare the effectiveness of Trichogramma japonicum Ashmead and insecticides (powder, 1.5% mevinphos + 3% alpha‐hexachloro cyclohexane) against the rice leaf roller based on the estimates of I produced by delete‐d jackknife and bootstrap techniques. At α= 0.05 level, there were statistical differences between wasp treatment and control, and between wasp and insecticide treatments, if the normality is ensured, or by the narrowest 95 percentiles. However, there is still no difference between insecticide treatment and control. By Z‐test criterion, wasp treatment is better than control and insecticide treatment with P‐value < 0.01. Insecticide treatment is similar to control with P‐value > 0.2 indicating that 95% confidence intervals procedure is more conservative. Although similar conclusions may be drawn by re‐sampling techniques, such as the delta method, about the suppressive effect of trichogramma and insecticides, the normality of the estimates can be checked and guaranteed, and the correlation among sequential life stages of rice leaf roller is also considered in the estimation. Judged by the P‐values from Z‐test, the delta method is more conservative.     

Assessing group differences in biodiversity by simultaneously testing a user‐defined selection of diversity indices

Comparing diversities between groups is a task biologists are frequently faced with, for example in ecological field trials or when dealing with metagenomics data. However, researchers often waver about which measure of diversity to choose as there is a multitude of approaches available. As Jost (2008, Molecular Ecology, 17 , 4015) has pointed out, widely used measures such as the Shannon or Simpson index have undesirable properties which make them hard to compare and interpret. Many of the problems associated with the use of these ‘raw’ indices can be corrected by transforming them into ‘true’ diversity measures. We introduce a technique that allows the comparison of two or more groups of observations and simultaneously tests a user‐defined selection of a number of ‘true’ diversity measures. This procedure yields multiplicity‐adjusted P‐values according to the method of Westfall and Young (1993, Resampling‐Based Multiple Testing: Examples and Methods for p‐Value Adjustment, 49, 941), which ensures that the rate of false positives (type I error) does not rise when the number of groups and/or diversity indices is extended. Software is available in the R package ‘simboot’.     

Molecular cloning and functional analyses of glutathione peroxidase homologous genes from Chlorella sp. NJ-18

The extreme variability of the mitochondrial (mito) genomes of bivalves makes it difficult to understand their evolutionary dynamics, given that species from different families do not share comparable features. We compared the mitogenomes from four Paphia clams (three of them were firstly sequenced) and found that mitogenome reorganization among the four congeneric species is not random but follows phylogenetic trends. Start/stop codon variations are species-correlated rather than gene-correlated, and bear useful phylogenetic information. Unique start/stop codon usage in P. euglypta and A+T content in P. amabilis indicates that these mitogenome-level characters, usually considered to be conservative features in other lineages, may not be phylogenetically evolved, but may have evolved via species-specific mitogenomic maintenance mechanisms. Variable divergence of two trnM genes in different lineages may demonstrate differences in mechanisms by which paralogous trnM genes are maintained. Sequence alignment analysis indicates that the VNTRs in the four mitogenomes have a common origin. The rationale of the subgenus Neotapes Kuroda and Habe, 1971 was supported by evidence from morphological characters, mitogenomic features, as well as phylogenetic analyses using cox1 and rrnS genes. The data suggest that the taxonomic basis of the subgenus should be “smooth surface” but not “undulated lines,” and P. textile should be classified to the Neotapes subgenus.     

An improved statistical method for detecting heterotachy in nucleotide sequences

The principle of heterotachy states that the substitution rate of sites in a gene can change through time. In this article, we propose a powerful statistical test to detect sites that evolve according to the process of heterotachy. We apply this test to an alignment of 1289 eukaryotic rRNA molecules to 1) determine how widespread the phenomenon of heterotachy is in ribosomal RNA, 2) to test whether these heterotachous sites are nonrandomly distributed, that is, linked to secondary structure features of ribosomal RNA, and 3) to determine the impact of heterotachous sites on the bootstrap support of monophyletic groupings. Our study revealed that with 21 monophyletic taxa, approximately two-thirds of the sites in the considered set of sequences is heterotachous. Although the detected heterotachous sites do not appear bound to specific structural features of the small subunit rRNA, their presence is shown to have a large beneficial influence on the bootstrap support of monophyletic groups. Using extensive testing, we show that this may not be due to heterotachy itself but merely due to the increased substitution rate at the detected heterotachous sites.     

Accounting for molecular stochasticity in systematic revisions: Species limits and phylogeny of Paroaria

Different frameworks have been proposed for using molecular data in systematic revisions, but there is ongoing debate on their applicability, merits and shortcomings. In this paper we examine the fit between morphological and molecular data in the systematic revision of Paroaria, a group of conspicuous songbirds endemic to South America. We delimited species based on examination of >600 specimens, and developed distance-gap, and distance- and character-based coalescent simulations to test species limits with molecular data. The morphological and molecular data collected were then analyzed using parsimony, maximum likelihood, and Bayesian phylogenetics. The simulations were better at evaluating the new species limits than using genetic distances. Species diversity within Paroaria had been underestimated by 60%, and the revised genus comprises eight species. Phylogenetic analyses consistently recovered a congruent topology for the most recently derived species in the genus, but the most basal divergences were not resolved with these data. The systematic and phylogenetic hypotheses developed here are relevant to both setting conservation priorities and understanding the biogeography of South America.