Testing small study effects in multivariate meta-analysis

Small study effects occur when smaller studies show different, often larger, treatment effects than large ones, which may threaten the validity of systematic reviews and meta-analyses. The most well-known reasons for small study effects include publication bias, outcome reporting bias, and clinical heterogeneity. Methods to account for small study effects in univariate meta-analysis have been extensively studied. However, detecting small study effects in a multivariate meta-analysis setting remains an untouched research area. One of the complications is that different types of selection processes can be involved in the reporting of multivariate outcomes. For example, some studies may be completely unpublished while others may selectively report multiple outcomes. In this paper, we propose a score test as an overall test of small study effects in multivariate meta-analysis. Two detailed case studies are given to demonstrate the advantage of the proposed test over various naive applications of univariate tests in practice. Through simulation studies, the proposed test is found to retain nominal Type I error rates with considerable power in moderate sample size settings. Finally, we also evaluate the concordance between the proposed tests with the naive application of univariate tests by evaluating 44 systematic reviews with multiple outcomes from the Cochrane Database.     

The surface structure of a complex substrate revealed by enzyme kinetics and Freundlich constants for α-amylase interaction with the surface of starch

Background

Starch is a main source of carbohydrate in human diets, but differences are observed in postprandial glycaemia following ingestion of different foods containing identical starch contents. Such differences reflect variations in rates at which different starches are digested in the intestine. In seeking explanations for these differences, we have studied the interaction of α-amylase with starch granules. Understanding this key step in digestion should help with a molecular understanding for observed differences in starch digestion rates.

Methods

For enzymes acting upon solid substrates, a Freundlich equation relates reaction rate to enzyme adsorption at the surface. The Freundlich exponent (n) equals 2/3 for a liquid-smooth surface interface, 1/3 for adsorption to exposed edges of ordered structures and 1.0 for solution–solution interfaces. The topography of a number of different starch granules, revealed by Freundlich exponents, was compared with structural data obtained by differential scanning calorimetry and Fourier transform infrared spectroscopy with attenuated total internal reflectance (FTIR-ATR).

Results

Enzyme binding rate and FTIR-ATR peak ratio were directly proportional to n and Δ_gelH was inversely related to n. Amylase binds fastest to solubilised starch and to granules possessing smooth surfaces at the solid–liquid interface and slowest to granules possessing ordered crystalline surfaces.

Conclusions

Freundlich exponents provide information about surface blocklet structures of starch that supplements knowledge obtained from physical methods.

General Significance

Nanoscale structures at the surface of starch granules influence hydrolysis by α-amylase. This can be important in understanding how dietary starch is digested with relevance to diabetes, cardiovascular health and cancer.     

Evolutionary patterns of the mitochondrial genome in the Moorish gecko, Tarentola mauritanica

A previous study on the evolutionary patterns of Tarentola mauritanica demonstrated that low levels of mitochondrial diversity observed in the European populations relative to nuclear markers were consistent with a selective sweep hypothesis. In order to unravel the mitochondrial evolutionary history in this European population and two other lineages of T. mauritanica (Iberian and North African clades), variation within 22 nearly complete mitogenomes was analyzed. Surprisingly, each clade seems to have a distinct evolutionary history; with both the European and Iberian clades presenting a decrease of polymorphism, which in the former is consistent with departure from neutrality of the mtDNA (positive or background selection), but in the latter seems to be the result of a bottleneck after a population expansion. The pattern exhibited by the North African clade seems to be a consequence of adaptation to certain mtDNA variants by positive selection.     

Functional and evolutionary analysis of DXL1, a non-essential gene encoding a 1-deoxy-D-xylulose 5-phosphate synthase like protein in Arabidopsis thaliana

The synthesis of 1-deoxy-D-xylulose 5-phosphate (DXP), catalyzed by the enzyme DXP synthase (DXS), represents a key regulatory step of the 2-C-methyl-D-erythritol 4-phosphate (MEP) pathway for isoprenoid biosynthesis. In plants DXS is encoded by small multigene families that can be classified into, at least, three specialized subfamilies. Arabidopsis thaliana contains three genes encoding proteins with similarity to DXS, including the well-known DXS1/CLA1 gene, which clusters within subfamily I. The remaining proteins, initially named DXS2 and DXS3, have not yet been characterized. Here we report the expression and functional analysis of A. thaliana DXS2. Unexpectedly, the expression of DXS2 failed to rescue Escherichia coli and A. thaliana mutants defective in DXS activity. Coherently, we found that DXS activity was negligible in vitro, being renamed as DXL1 following recent nomenclature recommendation. DXL1 is targeted to plastids as DXS1, but shows a distinct expression pattern. The phenotypic analysis of a DXL1 defective mutant revealed that the function of the encoded protein is not essential for growth and development. Evolutionary analyses indicated that DXL1 emerged from DXS1 through a recent duplication apparently specific of the Brassicaceae lineage. Divergent selective constraints would have affected a significant fraction of sites after diversification of the paralogues. Furthermore, amino acids subjected to divergent selection and likely critical for functional divergence through the acquisition of a novel, although not yet known, biochemical function, were identified. Our results provide with the first evidences of functional specialization at both the regulatory and biochemical level within the plant DXS family.     

Hemocyanin gene family evolution in spiders (Araneae), with implications for phylogenetic relationships and divergence times in the infraorder Mygalomorphae

Hemocyanins are multimeric copper-containing hemolymph proteins involved in oxygen binding and transport in all major arthropod lineages. Most arachnids have seven primary subunits (encoded by paralogous genes a–g), which combine to form a 24-mer (4 × 6) quaternary structure. Within some spider lineages, however, hemocyanin evolution has been a dynamic process with extensive paralog duplication and loss. We have obtained hemocyanin gene sequences from numerous representatives of the spider infraorders Mygalomorphae and Araneomorphae in order to infer the evolution of the hemocyanin gene family and estimate spider relationships using these conserved loci. Our hemocyanin gene tree is largely consistent with the previous hypotheses of paralog relationships based on immunological studies, but reveals some discrepancies in which paralog types have been lost or duplicated in specific spider lineages. Analyses of concatenated hemocyanin sequences resolved deep nodes in the spider phylogeny and recovered a number of clades that are supported by other molecular studies, particularly for mygalomorph taxa. The concatenated data set is also used to estimate dates of higher-level spider divergences and suggests that the diversification of extant mygalomorphs preceded that of extant araneomorphs. Spiders are diverse in behavior and respiratory morphology, and our results are beneficial for comparative analyses of spider respiration. Lastly, the conserved hemocyanin sequences allow for the inference of spider relationships and ancient divergence dates.     

Analysis of genetic variability and relationships among Mentha L. using the limonene synthase gene, LS

The genus Mentha comprises a group of aromatic plants with worldwide distribution. Because of frequent interspecific hybridization, the genetic relationships within the genus are not clearly understood. Limonene synthase, which catalyses the first committed step in the essential oil monoterpene biosynthetic pathway, is considered to be a possible rate limiting enzyme. With the homology-based cloning method, primers were designed according to cDNA sequence to amplify full-length DNA sequences in 13 Mentha samples from five species, using Perilla as an outgroup. Analyses of gene structure, length variation, GC-content, Ts/Tv ratio and evolutionary diversity were carried out. Consensus phylogenetic trees were obtained using maximum likelihood, neighbor-joining, and maximum parsimony, respectively, based on the full-length genomic DNA sequences, complete ORF coding sequences and predicted amino acid sequences. The results presented here based on the sequence of MhLS provide the first credibly supported genetic relationships for Mentha, which enables a basis for further mint taxonomy, cultivation and breeding.     

Myeloperoxidase G-463A polymorphism and susceptibility to coronary artery disease: A meta-analysis

Published data on the association between the myeloperoxidase (MPO) G-463A polymorphism and coronary artery disease (CAD) are inconclusive. To derive a more precise estimation of the relationship, a meta-analysis on this topic was performed. PubMed, EMBASE and Chinese national knowledge infrastructure were searched for studies regarding the association between the MPO G-463A polymorphism and CAD. A logistic regression analysis was used to estimate the genetic effect and the possible genetic model of action. Summary odds ratios (ORs) with their corresponding 95% confidence intervals (CIs) were calculated. There was strong evidence for an association between the MPO G-463A polymorphism and CAD. The genetic model of action was most likely to be co-dominant. Overall, the data showed that AA and GA genotypes were significantly associated with reduced risk of CAD (AA vs. GG: OR = 0.37, 95% CI = 0.17–0.78; GA vs. GG: OR = 0.73, 95% CI = 0.57–0.92). In subgroup analyses by study population and sources of controls, statistically significant results were observed in the Chinese population (AA vs. GG: OR = 0.21, 95% CI = 0.10–0.43; GA vs. GG: OR = 0.57, 95% CI =0.44–0.74) and in hospital-based control studies (AA vs. GG: OR = 0.20, 95% CI = 0.10–0.39; GA vs. GG: OR = 0.61, 95% CI = 0.48–0.77). This meta-analysis suggests that the MPO G-463A variant genotypes may be associated with decreased risk of CAD. However, given the limited number of studies and the potential biases, the influence of this polymorphism on CAD risk needs further investigation.     

Phylogenetics of Escallonia (Escalloniaceae) based on plastid DNA sequence data

Escallonia (Escalloniaceae) is a New World genus of c. 39 species distributed mainly in the South American highlands. Plastid DNA sequence data from the intergenic spacers trnS‐trnG and 3′ trnV‐ndhC and the ndhF gene for 32 species were used to examine the relationships among species and related genera and to analyse the relationship between phylogeny and the geographical distribution of the species. Maximum parsimony and Bayesian inference were employed to analyse the data. The sister relationship of Escallonia to Forgesia and Valdivia was corroborated. We recovered five strongly supported clades that are geographically structured, suggesting that the evolutionary history of the genus may be linked to historical processes, including the uplift of mountainous systems in South America. © 2013 The Linnean Society of London, Botanical Journal of the Linnean Society, 2013, 173 , 442–451.     

Responses of epidermal cell turgor pressure and photosynthetic activity of leaves of the atmospheric epiphyte Tillandsia usneoides (Bromeliaceae) after exposure to high humidity

It has been well-established that many epiphytic bromeliads of the atmospheric-type morphology, i.e., with leaf surfaces completely covered by large, overlapping, multicellular trichomes, are capable of absorbing water vapor from the atmosphere when air humidity increases. It is much less clear, however, whether this absorption of water vapor can hydrate the living cells of the leaves and, as a consequence, enhance physiological processes in such cells. The goal of this research was to determine if the absorption of atmospheric water vapor by the atmospheric epiphyte Tillandsia usneoides results in an increase in turgor pressure in leaf epidermal cells that subtend the large trichomes, and, by using chlorophyll fluorescence techniques, to determine if the absorption of atmospheric water vapor by leaves of this epiphyte results in increased photosynthetic activity. Results of measurements on living cells of attached leaves of this epiphytic bromeliad, using a pressure probe and of whole-shoot fluorescence imaging analyses clearly illustrated that the turgor pressure of leaf epidermal cells did not increase, and the photosynthetic activity of leaves did not increase, following exposure of the leaves to high humidity air. These results experimentally demonstrate, for the first time, that the absorption of water vapor following increases in atmospheric humidity in atmospheric epiphytic bromeliads is mostly likely a physical phenomenon resulting from hydration of non-living leaf structures, e.g., trichomes, and has no physiological significance for the plant's living tissues.     

Development of thermodynamic optimum searching (TOS) to improve the prediction accuracy of flux balance analysis

Flux balance analysis (FBA) has been widely used in calculating steady‐state flux distributions that provide important information for metabolic engineering. Several thermodynamics‐based methods, for example, quantitative assignment of reaction directionality and energy balance analysis have been developed to improve the prediction accuracy of FBA. However, these methods can only generate a thermodynamically feasible range, rather than the most thermodynamically favorable solution. We therefore developed a novel optimization method termed as thermodynamic optimum searching (TOS) to calculate the thermodynamically optimal solution, based on the second law of thermodynamics, the minimum magnitude of the Gibbs free energy change and the maximum entropy production principle (MEPP). Then, TOS was applied to five physiological conditions of Escherichia coli to evaluate its effectiveness. The resulting prediction accuracy was found significantly improved (10.7–48.5%) by comparing with the ¹³C‐fluxome data, indicating that TOS can be considered an advanced calculation and prediction tool in metabolic engineering. Biotechnol. Bioeng. 2013; 110: 914–923. © 2012 Wiley Periodicals, Inc.