Weighted bootstrapping: a correction method for assessing the robustness of phylogenetic trees

Background  

Non-parametric bootstrapping is a widely-used statistical procedure for assessing confidence of model parameters based on the empirical distribution of the observed data [1] and, as such, it has become a common method for assessing tree confidence in phylogenetics [2]. Traditional non-parametric bootstrapping does not weigh each tree inferred from resampled (i.e., pseudo-replicated) sequences. Hence, the quality of these trees is not taken into account when computing bootstrap scores associated with the clades of the original phylogeny. As a consequence, traditionally, the trees with different bootstrap support or those providing a different fit to the corresponding pseudo-replicated sequences (the fit quality can be expressed through the LS, ML or parsimony score) contribute in the same way to the computation of the bootstrap support of the original phylogeny.     

Relationships between photosystem II efficiency and photochemical reflectance index under different levels of illumination: comparison among species grown at high- and low elevations through different seasons

Previously, we found a significant association between photosystem II efficiency (ΦPSII) and photochemical reflectance index (PRI) measured at predawn among different species at different elevations and throughout several seasons. However, this relationship has not been evaluated under varied levels of illumination. Here, we used the Taiwan species Pinus taiwanensis (a conifer distributed at 750–3,000 m a.s.l.), Stranvaesia niitakayamensis (an evergreen tree, 1,700–3,100 m) and two Miscanthus spp. (perennial C₄ Gramineae, coastline–3,200 m) to elucidate the ΦPSII–PRI relationship. We studied six levels of photosynthetic photon flux density (PPFD) (0, 200, 400, 800, 1,200 and 2,000 μmol m⁻² s⁻¹) over several growth seasons at high (2,600 m a.s.l.) and low (800 m a.s.l.) elevation sites. In comparing the same species or genus, ΦPSII and PRI were closely correlated in darkness or under the same level of PPFD, with data obtained from different seasons and elevations pooled for regression analysis. Because both the intercept and slope of the ΦPSII–PRI equation showed a negative curvilinear correlation with PPFD, we could fit an empirical regression model, ΦPSII = c + d·ln(PPFD) + e·[ln(PPFD)]² + f·PRI + g·PRI·ln(PPFD) + h·PRI·[ln(PPFD)]², for multiple regression analysis. Using this model, we found a close correlation between the estimated and measured ΦPSII (r ² = 0.842−0.937, P < 0.001) for all four species examined and for mango (Mangifera indica) measured under both artificial illumination and sunlight (data from Weng et al. 2010). This empirical regression model could simulate both seasonal and diurnal variations of leaf-scale photosynthetic efficiency at high and low elevations.     

Across-species patterns of genetic variation in forest trees of Central Europe

The study focuses on geographical patterns of genetic variation at allozyme loci common for four main tree species of Central Europe (Norway spruce, silver fir, common beech and sessile oak). Moving-window averaging of four indicators of allelic richness and diversity (proportion of polymorphic loci, mean number of alleles per locus, effective number of alleles and expected heterozygosity) with window size of 50 × 50 km was used to identify the patterns. Moreover, local genetic divergence was assessed using the G _ST (Nei, Molecular population genetic and evolution, Amsterdam and Oxford, North-Holland, 1975) and D _j (Gregorius and Roberds, Theor Appl Genet 71:826–834, 1986) statistics for common beech and silver fir, where raw genotype data were available. Spatial patterns of diversity and allelic richness were quite similar. Romanian Carpathians were identified as the most important hotspot of genetic diversity and evolutionary divergence in Central Europe. Implications for genetic conservation are briefly discussed.     

Kinetic models of photosystem II should accommodate the effect of donor side quenching on variable chlorophyll <Emphasis Type="Italic">a</Emphasis> fluorescence in the microseconds time range

Quantitative data on laser flash-induced variable fluorescence in the 100 ns to 1 ms time range (Belyaeva et al. in Photosynth Res 98:105–119, 2008) confirming those of others (Steffen et al. in Biochemistry 40:173–180, 2001, Biochemistry 44:3123–3132, 2005; Belyaeva et al. in Biophysics 51(6):976–990, 2006), need a substantial correction with respect to magnitude of the normalized variable fluorescence associated with single turnover-induced charge separation in RCs of PS II. Their data are conclusive with the involvement of donor side quenching, the release of which occurs with a rate constant in the range of tens of ms⁻¹, and presumed to be associated with reduction of Y_\textz ⁺ Y_{\text{z}}^{ + } by the OEC.     

The use of parsimony network analysis for the formal delineation of phylogenetic species of yeasts: Candida apicola, Candida azyma, and Candida parazyma sp. nov., cosmopolitan yeasts associated with floricolous insects

Parsimony network analysis of rDNA sequences was used to delimit phylogenetic species of yeasts in an objective, formal manner. Many strains assigned to Candida apicola (Starmerella clade), when compared to the type, fell outside the inclusion limits proposed by Kurtzman and Robnett (1998) based on a pair-wise comparison of the large subunit rRNA gene D1/D2 domains. However, when these sequences were analyzed jointly with ITS rDNA sequences by parsimony network analysis, 28 of the 30 strains formed a cohesive set. Two strains, MUCL 45721 and CBS 4353, were excluded from the species, but there was no evident justification to subdivide the rest. A similar analysis of 81 isolates originally assigned to Candida azyma (Wickerhamiella clade) yielded dramatically different results, giving rise to six independent networks corresponding to Candida azyma sensu stricto (18 strains), Candida azymoides (2 strains), a pair of isolates from Australian hibiscus flowers, a single isolate from the same substrate, a single isolate from Malaysian bertam palm nectar, and 57 isolates that are assigned to the new species Candida parazyma (type = UWOPS 91-652.1^T = CBS 11563^T = NRRL Y-48669^T). The strains retained in C. azyma sensu stricto differed from one another by up to four substitutions in their D1/D2 sequences, but their polymorphism at the level of the ITS was considerable and suggested a history of divergence resulting from dispersal. Strains of C. parazyma fell into seven variant haplotypes based on sequences of the rDNA ITS and D1/D2 regions. The most abundant haplotype occurred across the global range of the species. Others were either endemic to Belize, Costa Rica, Rarotonga, or Tennessee, suggestive of vicariance, or occurred across remote localities, offering partial support to the notion of rapid dispersal.     

Phylogenetic position of the monotypic Des Murs’ Wiretail ( Sylviorthorhynchus desmursii , Aves: Furnariidae) based on mitochondrial and nuclear DNA

Parallelisms and paraphyletic assemblages are common among ovenbirds. Molecular markers are therefore the best approach when studying the evolutionary relationships among the members of this unparalleled diversified family. We obtained nucleotide sequence data from mitochondrial (cytochrome b) and nuclear genes (myoglobin and glyceraldehyde-3-phosphodehydrogenase) and used these to deduce the phylogenetic position of a monotypic genus endemic to the austral temperate rainforests of southern South America, the Des Murs’ Wiretail (Sylviorthorhynchus desmursii Des Murs, 1847, Aves: Furnariidae). Phylogenetic analyses based on maximum parsimony, maximum likelihood and Bayesian inference all converged into a congruent topology, with a basal position of Des Murs’ Wiretail within Synallaxinae together with Tit–Spinetails (Leptasthenura). Our data reject the hypothesis of a phylogenetic relationship between Des Murs’ Wiretail and thistletails (Schizoeaca) which exhibit parallelisms in morphology, tail structure and nest architecture. Using a molecular clock based on the myoglobin intron 2 gene, we estimated a divergence time of Des Murs’ Wiretail from Tit-Spinetails of 14–15 Myr, which is associated with the appearance of sclerophyllous forest elements in Chile at the Middle–Upper Miocene. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Phylogeny of the sabertoothed felids (Carnivora: Felidae: Machairodontinae)

In recent years, advances in our understanding of feline relationships have cast light on their evolutionary history. In contrast, there have been no phylogenetic analyses on machairodont felids, making it difficult to develop an evolutionary hypothesis based on the recent surge of studies on their craniomandibular morphology and functional anatomy. In this paper, I provide the first phylogenetic hypothesis of machairodont relationships based on 50 craniomandibular and dental characters from a wide range of sabercats spanning more 11 Myr. Exact searches produced 19 most‐parsimonious trees, and a strict consensus was well resolved. The Machairodontinae comprise a number of basal taxa (Promegantereon, Machairodus, Nimravides, Dinofelis, Metailurus) and a well‐supported clade of primarily Plio‐Pleistocene taxa (Megantereon, Smilodon, Amphimachairodus, Homotherium, Xenosmilus) for which the name Eumachairodontia taxon novum is proposed. Previous phenetic grouping of machairodont taxa into three distinct groups, the Smilodontini, Homotherini and Metailurini, was not supported by cladistic parsimony analysis, and forcing monophyly of these groups was significantly incompatible with character distribution. Machairodonts as a clade are not characterized by saberteeth, i.e. hypertrophied, blade‐like upper canines, but by small lower canines, as well as small M¹; and large P³ parastyle. True saberteeth arose later and are a synapomorphy of the Eumachairodontia.     

Axillary shoot proliferation and tuberization of <Emphasis Type="Italic">Dioscorea fordii</Emphasis> Prain et Burk

Dioscorea fordii Prain et Burk., a member of the family Dioscoreaceae, is an important tuberous food in China owing to its edible and medicinal functions. However, the lack of healthy planting material has restricted its production. The present study was an attempt to develop a protocol for in vitro propagation of D. fordii. In an initial two by two (2²) factorial experiment, the effects of culture system and activated charcoal have been observed for axillary shoot proliferation and tuberization of in vitro plantlets. Shoot length, frequency of proliferation, fresh weight and dry weight of shoots, frequency of tuberization and mean number of tubers per plantlet (NTPs) were significantly (P ≤ 0.05) greater in liquid medium compared to semi-solid media. Secondly, an orthogonal experimental design [L₉ (3⁴)] was used to investigate the effects of the ratio of naphthalene acetic acid (NAA) to 6-benzyladenine (BA) (NAA/BA), paclobutrazol, methyl jasmonate and sucrose concentration on tuberization. Sucrose concentration had the record effect on frequency of tuberization, NTPs and frequency of proliferation. The preferred medium for axillary shoot proliferation and tuberization of D. fordii found to be MS basal medium (Murashige and Skoog in Physiol Plant 15:473–479, 1962) supplemented with 1.0 mg l⁻¹ BA, 0.1 mg l⁻¹ NAA, 30 g l⁻¹ sucrose and 1.5 g l⁻¹ AC in liquid culture.     

Distinguishing HIV-1 drug resistance, accessory, and viral fitness mutations using conditional selection pressure analysis of treated versus untreated patient samples

Background  

HIV can evolve drug resistance rapidly in response to new drug treatments, often through a combination of multiple mutations [1–3]. It would be useful to develop automated analyses of HIV sequence polymorphism that are able to predict drug resistance mutations, and to distinguish different types of functional roles among such mutations, for example, those that directly cause drug resistance, versus those that play an accessory role. Detecting functional interactions between mutations is essential for this classification. We have adapted a well-known measure of evolutionary selection pressure (K _a /K _s ) and developed a conditional K _a /K _s approach to detect important interactions.     

Short-term coherency between gross primary production and community respiration in an algal-dominated reef flat

Rates of net community carbon production (mmol C m⁻² h⁻¹) were measured continuously in an algal-dominated reef flat community on the Kaneohe Bay barrier reef, Hawaii, for 12 days at the end of October 2006. The weather became increasingly cloudy during the last 5 days of measurements, resulting in a sevenfold decline in daily incident light (28–4 Ein m⁻² d⁻¹). In response, gross primary production (P) for the reef flat community also decreased sevenfold, varying linearly with light (r ² = 0.92, n = 12). Community respiration (R) decreased fivefold over this same period and was highly correlated with changes in P (r ² = 0.84, n = 12). We reason that this short-term coherence between P and R indicates that most of the carbon fixed during this period was rapidly metabolized via plant respiration. We further conclude that the dominance of autotrophic respiration under general conditions of nutrient-limited growth can explain much of the balance between P and R that is commonly observed in shallow reef communities.     

Bootstrap-based Support of HGT Inferred by Maximum Parsimony

Background  

Maximum parsimony is one of the most commonly used criteria for reconstructing phylogenetic trees. Recently, Nakhleh and co-workers extended this criterion to enable reconstruction of phylogenetic networks, and demonstrated its application to detecting reticulate evolutionary relationships. However, one of the major problems with this extension has been that it favors more complex evolutionary relationships over simpler ones, thus having the potential for overestimating the amount of reticulation in the data. An ad hoc solution to this problem that has been used entails inspecting the improvement in the parsimony length as more reticulation events are added to the model, and stopping when the improvement is below a certain threshold.     

Purification of the photosynthetic reaction center from <Emphasis Type="Italic">Heliobacterium modesticaldum</Emphasis>

We have developed a purification protocol for photoactive reaction centers (HbRC) from Heliobacterium modesticaldum. HbRCs were purified from solubilized membranes in two sequential chromatographic steps, resulting in the isolation of a fraction containing a single polypeptide, which was identified as PshA by LC–MS/MS of tryptic peptides. All polypeptides reported earlier as unknown proteins (in Heinnickel et al., Biochemistry 45:6756–6764, 2006; Romberger et al., Photosynth Res 104:293–303, 2010) are now identified by mass spectrometry to be the membrane-bound cytochrome c ₅₅₃ and four different ABC-type transporters. The purified PshA homodimer binds the following pigments: 20 bacteriochlorophyll (BChl) g, two BChl g′, two 8¹-OH-Chl a _F, and one 4,4′-diaponeurosporene. It lacks the PshB polypeptide binding the F_A and F_B [4Fe–4S] clusters. It is active in charge separation and exhibits a trapping time of 23 ps, as judged by time-resolved fluorescence studies. The charge recombination rate of the P₈₀₀ ⁺F_X⁻ state is 10–15 ms, as seen before. The purified HbRC core was able to reduce cyanobacterial flavodoxin in the light, exhibiting a K _M of 10 μM and a k _cat of 9.5 s⁻¹ under near-saturating light. There are ~1.6 menaquinones per HbRC in the purified complex. Illumination of frozen HbRC in the presence of dithionite can cause creation of a radical at g = 2.0046, but this is not a semiquinone. Furthermore, we show that high-purity HbRCs are very stable in anoxic conditions and even remain active in the presence of oxygen under low light.     

Mitochondria from Dipodascus (Endomyces) magnusii and Yarrowia lipolytica yeasts did not undergo a Ca2+-dependent permeability transition even under anaerobic conditions

In this study we used tightly-coupled mitochondria from Yarrowia lipolytica and Dipodascus (Endomyces) magnusii yeasts. The two yeast strains are good alternatives to Saccharomyces cerevisiae, being aerobes containing well-structured mitochondria (thus ensuring less structural limitation to observe their appreciable swelling) and fully competent respiratory chain with three invariantly functioning energy conservation points, including Complex I, that can be involved in induction of the canonical Ca²⁺/P_i-dependent mitochondrial permeability transition (mPTP pore) with an increased open probability when electron flux increases (Fontaine et al. J Biol Chem 273:25734–25740, 1998; Bernardi et al. FEBS J 273:2077–2099, 2006). High-amplitude swelling and collapse of the membrane potential were used as parameters for demonstrating pore opening. Previously (Kovaleva et al. J Bioenerg Biomembr 41:239–249, 2009; Kovaleva et al. Biochemistry (Moscow) 75:297–303, 2010) we have shown that mitochondria from Y. lipolytica and D. magnusii were very resistant to the Ca²⁺ overload combined with varying concentrations of P_i, palmitic acid, SH-reagents, carboxyatractyloside (an inhibitor of ADP/ATP translocator), as well as depletion of intramitochondrial adenine nucleotide pools, deenergization of mitochondria, and shifting to acidic pH values in the presence of high [P_i]. Here we subjected yeast mitochondria to other conditions known to induce an mPTP in animal and plant mitochondria, namely to Ca²⁺ overload under hypoxic conditions (anaerobiosis). We were unable to observe Ca²⁺-induced high permeability of the inner membrane of D. magnusii and Y. lipolytica yeast mitochondria under anaerobic conditions, thus suggesting that an mPTP-like pore, if it ever occurs in yeast mitochondria, is not coupled with the Ca²⁺ uptake. The results provide the first demonstration of ATP-dependent energization of yeast mitochondria under conditions of anaerobiosis.     

TGF-Mediated Dynamics in a System of Many Coupled Nephrons

This paper presents a mathematical model of a system of many coupled nephrons branching from a common cortical radial artery, and accompanying analysis of that system. This modeling effort is a first step in understanding how coupling magnifies the tendency of nephrons to oscillate owing to tubuloglomerular feedback. Central to the present work is the single nephron integral model (as in Pitman et al., The IMA Volumes in Mathematics and Its Applications, vol. 129, pp. 345–364, 2002 and in Zaritski, Ph.D. Dissertation, 1999) which is a simplification of the single nephron PDE model of Layton et al. (Am. J. Physiol. 261, F904–F919, 1991). A second principal idea used in the present model is a coupling of model nephrons, generalizing the work of Pitman et al. (Bull. Math. Biol. 66, 1463–1492, 2004) who proposed a model of two coupled nephrons. In this study, we couple nephrons through a nearest neighbor interaction. Speaking generally, our results suggest that a series of similar nephrons coupled to their nearest neighbors are more prone to be found in an oscillatory mode, relative to a single nephron with the same properties. More specifically, we show analytically that, for N coupled identical nephrons, the region supporting oscillatory solutions in the time delay–gain parameter plane increases with N. Numerical simulations suggest that, if N nephrons have gains and time delays that do not differ by much, the system is, again, more prone to oscillate, relative to a single nephron, and the oscillations tend to be approximately synchronous and in-phase. We examine the effect of parameters on bifurcation. We also examine alternative models of coupling; this analysis allows us to conclude that the increased propensity of coupled nephrons to oscillate is a robust finding, true for several models of nephron interaction.     

<Emphasis Type="Italic">In vitro</Emphasis> propagation of <Emphasis Type="Italic">Chimonanthus praecox</Emphasis> (L.), a winter flowering ornamental shrub

A procedure for the micropropagation of Chimonanthus praecox (L) Link, wintersweet, has been developed using buds from adult trees excised in spring. Shoot cultures established on Murashige and Skoog (1962) medium supplemented with 0.5 mg l⁻¹ 6-benzyladenine (BAP) and 0.1 mg l⁻¹ indole-3-butyric acid (IBA) were difficult to maintain in vitro through extended periods of time due to browning of the medium, shoot and leaf necrosis, and hyperhydricity. A treatment combining the use of 0.1% w/v activated charcoal and addition of a double phase agar-solidified/liquid medium improved propagation, enabling a successful in vitro propagation scheme to be developed. Optimal shoot multiplication occurred on medium containing 0.5 mg l⁻¹ BAP, and rooting on medium with 2.0 mg l⁻¹ IBA for 7 d, followed by transfer to hormone-free medium. Rooted plantlets were easily acclimated in a glasshouse and replanted and cultured outdoors.     

Development of a fermentation process based on a defined medium for the production of pregallidermin, a nontoxic precursor of the lantibiotic gallidermin

In this work, a defined medium was developed and optimized for the mutant strain Staphylococcus gallinarum ΔP, which produces pregallidermin (PGDM), a nontoxic precursor of the lantibiotic gallidermin (GDM). The availability of a defined medium is a prerequisite for a rational process development and the investigation of medium effects on final product concentration, yield, and volumetric productivity. We identified four vitamins and three metal ions as essential for growth and PGDM production with S. gallinarum ΔP. The strain was capable of growing without any added amino acids, but the addition of proline had a strong growth-stimulatory effect. The concentrations of all essential compounds were balanced in a continuous culture using a medium-shift technique. Based on this balanced medium, a fed-batch process was developed in which S. gallinarum ΔP was grown up to a biomass concentration of 67 g l⁻¹ and produced 1.95 g l⁻¹ PGDM, equivalent to 0.57 mM. In the fermentation broth, we identified other GDM precursors in addition to those with a 12 or 14-amino-acid-long leader peptide that had been observed previously. Including those precursors with shorter leader sequences, the final concentration would correspond to 0.69 mM. In molar terms, this represents a roughly fourfold or fivefold increase, respectively, over established, complex medium-based gallidermin production processes (Kempf et al. 2000). With the same medium and feed protocol, the maximum concentration of mature GDM produced by wild-type S. gallinarum Tü 3928 was only 0.08 mM.     

Evaluation of green fluorescent protein as a reporter gene and phosphinothricin as the selective agent for achieving a higher recovery of transformants in cucumber (<Emphasis Type="Italic">Cucumis sativus</Emphasis> L. cv. Poinsett76) via <Emphasis Type="Italic">Agrobacterium tumefaciens</Emphasis>

Efficient Agrobacterium tumefaciens-mediated transformation and a higher recovery of transformed plants of cucumber cv. Poinsett76 were achieved via direct organogenesis from cotyledon explants. Stable transformants were obtained by inoculating explants with A. tumefaciens strains EHA105 or LBA4404, both harboring the binary vector pME508, which contains the neomycin phosphotransferase II (nptII) and phosphinothricin resistance genes (bar) conferring resistance to kanamycin and PPT, respectively, as selectable markers and the sgfp-tyg gene for the green fluorescent protein (GFP) as a visual marker driven by the constitutive CaMV35S promoter in the presence of acetosyringone (50 μM). Transformed shoots were obtained on MS Murashige and Skoog (Plant Physiol. 15: 473–497, 1962) medium supplemented with 1 mg L⁻¹ benzyladenine (BA), 20 mg L⁻¹ l-glutamine and 2 mg L⁻¹ phosphinothricin (PPT) or 100 mg L⁻¹ kanamycin. The regenerated shoots were examined in vivo using a hand-held long wave UV lamp for GFP expression. The GFP screening helped identify escapes and chimeric shoots at regular intervals to increase the growth of transformed shoots on cotyledon explants. Elongation and rooting of putative transformants were achieved on PPT (2 mg L⁻¹) containing MS media with 0.5 mg L⁻¹ gibberellic acid (GA₃) and 0.6 mg L⁻¹ indole butyric acid (IBA), respectively. PCR and Southern analyses confirmed the integration of the sgfp gene into the genome of T₀ and the progenies. T₁ segregation of transgenic progeny exhibited Mendelian inheritance of the transgene. The use of EHA105 resulted in 21% transformation efficiency compared to 8.5% when LBA4404 was used. This higher rate was greatly facilitated by PPT selection coupled with effective screening of transformants for GFP expression, thus making the protocol highly useful for the recovery of a higher number of transgenic cucumber plants.     

Variation in the determinants of power of chemically skinned type I rat soleus muscle fibres

We explored to which extent maximal velocity of shortening (Vmax), force per cross-sectional area (specific tension, Po) and curvature of the force–velocity relationship (a/Po in the Hill equation) contribute to differences in peak power of single, chemically skinned rat type I fibres. Force–velocity relationships were determined from isotonic contractions of 94 maximally activated fibres. Peak power (±SD) was 3.50 ± 1.64 W L⁻¹. There was a tenfold range of peak power and five-, six- and fourfold ranges for Po, Vmax and a/Po, respectively. None of the differences between fibres was explicable by differences in myosin heavy or light chain composition. The inverse relationship between a/Po and Vmax suggests a similar underlying cause. Fitting the data to the Huxley (Progr Biophys Biophys Chem 7:255–318, 1957) cross-bridge model showed that the rate constant g ₂ and the sum of the rate constants (f + g ₁) co-varied, both being low in the slowest fibres. Approximately 16% of the variation in Po could be explained by variation in the proportion of attached cycling cross-bridges (f/(f + g ₁)), but the origin of most of the variance in Po remains unknown.     

The expression of nr0b1, the earliest gene in zebrafish tooth development,is a marker for human tooth and ameloblastoma formation

Zebrafish teeth develop on pharyngeal jaws in the 5th branchial arch, but early tooth development is remarkably similar to mammals (Borday-Birraux et al., Evol Dev 8:130, 2006). Recently, eve1 has been shown to be associated with the primary tooth (4V¹) and early ameloblast development, the enamel organ precursor (Laurenti et al., Dev Dyn 230:727, 2004). dax1 is initially expressed in the 5th branchial arch in zebrafish at approximately 26 h postfertilization (hpf) and colocalizes with eve1 expression at ~48 hpf. Embryos injected with dax1 morpholino show downregulation of eve1 expression. Based on the zebrafish observations, we demonstrated novel DAX1 expression in normal human dental, benign ameloblastoma, and malignant ameloblastoma tissues. The association of NR0B1 and its protein product DAX1 with primary tooth development and ameloblastoma tumorigenesis is an association not previously described.