Bootstrapping phylogenies inferred from rearrangement data

ABSTRACT: Large-scale sequencing of genomes has enabled the inference of phylogenies based on the evolution of genomic architecture, under such events as rearrangements, duplications, and losses. Many evolutionary models and associated algorithms have been designed over the last few years and have found use in comparative genomics and phylogenetic inference. However, the assessment of phylogenies built from such data has not been properly addressed to date. The standard method used in sequence-based phylogenetic inference is the bootstrap, but it relies on a large number of homologous characters that can be resampled; yet in the case of rearrangements, the entire genome is a single character. Alternatives such as the jackknife suffer from the same problem, while likelihood tests cannot be applied in the absence of well established probabilistic models. We present a new approach to the assessment of distance-based phylogenetic inference from whole-genome data; our approach combines features of the jackknife and the bootstrap and remains nonparametric. For each feature of our method, we give an equivalent feature in the sequence-based framework; we also present the results of extensive experimental testing, in both sequence-based and genome-based frameworks. Through the feature-by-feature comparison and the experimental results, we show that our bootstrapping approach is on par with the classic phylogenetic bootstrap used in sequence-based reconstruction, and we establish the clear superiority of the classic bootstrap for sequence data and of our corresponding new approach for rearrangement data over proposed variants. Finally, we test our approach on a small dataset of mammalian genomes, verifying that the support values match current thinking about the respective branches. Our method is the first to provide a standard of assessment to match that of the classic phylogenetic bootstrap for aligned sequences. Its support values follow a similar scale and its receiver-operating characteristics are nearly identical, indicating that it provides similar levels of sensitivity and specificity. Thus our assessment method makes it possible to conduct phylogenetic analyses on whole genomes with the same degree of confidence as for analyses on aligned sequences. Extensions to search-based inference methods such as maximum parsimony and maximum likelihood are possible, but remain to be thoroughly tested.     

Contribution of sucrose synthase, ADP-glucose pyrophosphorylase and starch synthase to starch synthesis in developing pea seeds

Using genetic variability existing amongst nine pea genotypes (Pisum sativum L.), the biochemical basis of sink strength in developing pea seeds was investigated. Sink strength was considered to be reflected by the rate of starch synthesis (RSS) in the embryo, and sink activity in the seed was reflected by the relative rate of starch synthesis (RRSS). These rates were compared to the activities of three enzymes of the starch biosynthetic pathway [sucrose synthase (Sus), ADP-glucose pyrophosphorylase and starch synthase] at three developmental stages during seed filling (25, 50 and 75% of the dry seed weight). Complete sets of data collected during seed filling for the nine genotypes showed that, for all enzyme activities (expressed on a protein basis), only Sus in the embryo and seed coat was linearly and significantly correlated to RRSS. The contribution of the three enzyme activities to the variability in RSS and RRSS was evaluated by multiple regression analysis for the first two developmental stages. Only Sus activity in the embryo could explain, at least in part, the significant variability observed for both the RSS and the RRSS at each developmental stage. We conclude that Sus activity is a reliable marker of sink activity in developing pea seeds.     

Sea urchin Hox genes: insights into the ancestral Hox cluster

We describe the Hox cluster in the radially symmetric sea urchin and compare our findings to what is known from clusters in bilaterally symmetric animals. Several Hox genes from the direct-developing sea urchin Heliocidaris erythrogramma are described. CHEF gel analysis shows that the Hox genes are clustered on a < or = 300 kilobase (kb) fragment of DNA, and only a single cluster is present, as in lower chordates and other nonvertebrate metazoans. Phylogenetic analyses of sea urchin, amphioxus, Drosophila, and selected vertebrate Hox genes confirm that the H. erythrogramma genes, and others previously cloned from other sea urchins, belong to anterior, central, and posterior groups. Despite their radial body plan and lack of cephalization, echinoderms retain at least one of the anterior group Hox genes, an orthologue of Hox3. The structure of the echinoderm Hox cluster suggests that the ancestral deuterostome had a Hox cluster more similar to the current chordate cluster than was expected Sea urchins have at least three Abd-B type genes, suggesting that Abd-B expansion began before the radiation of deuterostomes.      

Smoking cessation and bronchial epithelial remodelling in COPD: a cross-sectional study

Background

Chronic Obstructive Pulmonary Disease (COPD) is associated with bronchial epithelial changes, including squamous cell metaplasia and goblet cell hyperplasia. These features are partially attributed to activation of the epidermal growth factor receptor (EGFR). Whereas smoking cessation reduces respiratory symptoms and lung function decline in COPD, inflammation persists. We determined epithelial proliferation and composition in bronchial biopsies from current and ex-smokers with COPD, and its relation to duration of smoking cessation.

Methods

114 COPD patients were studied cross-sectionally: 99 males/15 females, age 62 ± 8 years, median 42 pack-years, no corticosteroids, current (n = 72) or ex-smokers (n = 42, median cessation duration 3.5 years), postbronchodilator FEV₁ 63 ± 9% predicted. Squamous cell metaplasia (%), goblet cell (PAS/Alcian Blue⁺) area (%), proliferating (Ki-67⁺) cell numbers (/mm basement membrane), and EGFR expression (%) were measured in intact epithelium of bronchial biopsies.

Results

Ex-smokers with COPD had significantly less epithelial squamous cell metaplasia, proliferating cell numbers, and a trend towards reduced goblet cell area than current smokers with COPD (p = 0.025, p = 0.001, p = 0.081, respectively), but no significant difference in EGFR expression. Epithelial features were not different between short-term quitters (<3.5 years) and current smokers. Long-term quitters (≥3.5 years) had less goblet cell area than both current smokers and short-term quitters (medians: 7.9% vs. 14.4%, p = 0.005; 7.9% vs. 13.5%, p = 0.008; respectively), and less proliferating cell numbers than current smokers (2.8% vs. 18.6%, p < 0.001).

Conclusion

Ex-smokers with COPD had less bronchial epithelial remodelling than current smokers, which was only observed after long-term smoking cessation (>3.5 years).

Trial registration

NCT00158847     

Chemopreventive and renal protective effects for docosahexaenoic acid (DHA): implications of CRP and lipid peroxides

Background

The fish oil-derived ω-3 fatty acids, like docosahexanoic (DHA), claim a plethora of health benefits. We currently evaluated the antitumor effects of DHA, alone or in combination with cisplatin (CP) in the EAC solid tumor mice model, and monitored concomitant changes in serum levels of C-reactive protein (CRP), lipid peroxidation (measured as malondialdehyde; MDA) and leukocytic count (LC). Further, we verified the capacity of DHA to ameliorate the lethal, CP-induced nephrotoxicity in rats and the molecular mechanisms involved therein.

Results

EAC-bearing mice exhibited markedly elevated LC (2-fold), CRP (11-fold) and MDA levels (2.7-fold). DHA (125, 250 mg/kg) elicited significant, dose-dependent reductions in tumor size (38%, 79%; respectively), as well as in LC, CRP and MDA levels. These effects for CP were appreciably lower than those of DHA (250 mg/kg). Interestingly, DHA (125 mg/kg) markedly enhanced the chemopreventive effects of CP and boosted its ability to reduce serum CRP and MDA levels. Correlation studies revealed a high degree of positive association between tumor growth and each of CRP (r = 0.85) and leukocytosis (r = 0.89), thus attesting to a diagnostic/prognostic role for CRP. On the other hand, a single CP dose (10 mg/kg) induced nephrotoxicity in rats that was evidenced by proteinuria, deterioration of glomerular filtration rate (GFR, -4-fold), a rise in serum creatinine/urea levels (2–5-fold) after 4 days, and globally-induced animal fatalities after 7 days. Kidney-homogenates from CP-treated rats displayed significantly elevated MDA- and TNF-α-, but reduced GSH-, levels. Rats treated with DHA (250 mg/kg, but not 125 mg/kg) survived the lethal effects of CP, and showed a significant recovery of GFR; while their homogenates had markedly-reduced MDA- and TNF-α-, but -increased GSH-levels. Significant association was detected between creatinine level and those of MDA (r = 0.81), TNF-α ) r = 0.92) and GSH (r = -0.82); implying causal relationships.

Conclusion

DHA elicited prominent chemopreventive effects on its own, and appreciably augmented those of CP as well. The extent of tumor progression in various mouse groups was highly reflected by CRP levels (thus implying a diagnostic/prognostic role for CRP). Further, this study is the first to reveal that DHA can obliterate the lethal CP-induced nephrotoxicity and renal tissue injury. At the molecular level, DHA appears to act by reducing leukocytosis, systemic inflammation, and oxidative stress.     

Horizon-Specific Bacterial Community Composition of German Grassland Soils,as Revealed by Pyrosequencing-Based Analysis of 16S rRNA Genes

The diversity of bacteria in soil is enormous, and soil bacterial communities can vary greatly in structure. Here, we employed a pyrosequencing-based analysis of the V2-V3 16S rRNA gene region to characterize the overall and horizon-specific (A and B horizons) bacterial community compositions in nine grassland soils, which covered three different land use types. The entire data set comprised 752,838 sequences, 600,544 of which could be classified below the domain level. The average number of sequences per horizon was 41,824. The dominant taxonomic groups present in all samples and horizons were the Acidobacteria, Betaproteobacteria, Actinobacteria, Gammaproteobacteria, Alphaproteobacteria, Deltaproteobacteria, Chloroflexi, Firmicutes, and Bacteroidetes. Despite these overarching dominant taxa, the abundance, diversity, and composition of bacterial communities were horizon specific. In almost all cases, the estimated bacterial diversity (H′) was higher in the A horizons than in the corresponding B horizons. In addition, the H′ was positively correlated with the organic carbon content, the total nitrogen content, and the C-to-N ratio, which decreased with soil depth. It appeared that lower land use intensity results in higher bacterial diversity. The majority of sequences affiliated with the Actinobacteria, Bacteroidetes, Cyanobacteria, Fibrobacteres, Firmicutes, Spirochaetes, Verrucomicrobia, Alphaproteobacteria, Betaproteobacteria, and Gammaproteobacteria were derived from A horizons, whereas the majority of the sequences related to Acidobacteria, Chloroflexi, Gemmatimonadetes, Nitrospira, TM7, and WS3 originated from B horizons. The distribution of some bacterial phylogenetic groups and subgroups in the different horizons correlated with soil properties such as organic carbon content, total nitrogen content, or microbial biomass.Soil is probably the most complex microbial environment on Earth with respect to species richness and community size. The microbial richness in soils exceeds that of other environments (44) and is higher by orders of magnitude than the biodiversity of plants and animals. Cultivated soil or grassland soil contains an estimated 2 × 10⁹ prokaryotic cells per gram (12). Soil microbial communities are an important factor of agriculturally managed systems, as they are responsible for most nutrient transformations in soil and influence the above-ground plant diversity and productivity (53).To analyze the bacterial community in soils, most approaches target the 16S rRNA gene by PCR amplification and subsequent analysis employing sequencing of clone libraries (10, 24), denaturing gradient gel electrophoresis (DGGE) (38), or terminal restriction fragment length polymorphism (T-RFLP) (17, 52). Most of these approaches provided limited insights into the structure of soil bacterial communities, as the survey sizes and the number of compared sampling sites were small with respect to the enormous bacterial diversity present in different soil samples. For example, the reported clone libraries vary considerably in size, but small sample sizes (500 or fewer 16S rRNA gene sequences) are usually analyzed and employed for the theoretical estimation of species richness (39). This provides snapshots of the predominant bacterial community members, but phylogenetic groups that are present in a low abundance and which may possess important ecosystem functions are not assessed (47). In addition, it has been shown that rich sampling (several thousands of clones) of complex bacterial communities is required to perform robust measurements and estimations of community diversity parameters (37). Thus, the detection bias accompanying analyses of small sample sizes can lead to invalidated assumptions. Genetic profiling techniques such as DGGE and T-RFLP have high-throughput capability. These approaches allow researchers to unravel differences in community structure but are limited for assessing diversity (23, 40). To deeply survey the diversity and the composition of the bacterial communities within different soil samples, large-scale pyrosequencing of partial 16S rRNA genes has been employed recently. Previous pyrosequencing-based studies of soil (1, 30, 34, 43) have generated large data sets, which comprised 39,707 (30) to 152,359 (34) 16S rRNA partial gene sequences. Those studies provided comprehensive insights into the biogeography of bacterial soil communities and taxa that were present in a low abundance. However, all those studies focused on the analysis of microbial communities present in topsoil. The subsoil is also known to harbor an important part of the soil microbial biomass (18). It has been shown that the microbial population in the shallow subsurface is impacted by agricultural production to a similar extent as that in topsoil (5).In this study, we performed large-scale pyrosequencing-based analyses of 16S rRNA genes to assess the bacterial community composition in topsoil and the corresponding subsoil of nine different grassland sites in the Hainich region (Thuringia, Germany). To provide a high level of coverage at the species level (97% genetic distance) and minimize detection bias, we exceeded the above-described numbers of analyzed 16S rRNA gene sequences (752,838 in this study). To examine the impact of land use on bacterial diversity and community composition, the selected grassland sites covered a range of three different land use types, including samples from unfertilized pastures grazed by cattle, fertilized mown pastures grazed by cattle, and fertilized meadows. In many recent studies, surveys were focused on comprehensive analyses of a single soil or a few soil samples (1, 14, 37, 43). This allowed the determination of overall bacterial species richness and community composition, but the assessment of spatial patterns and environmental factors that drive these patterns is hampered by the limited number of examined soils. To assess spatial distribution and the impact of soil edaphic factors and land use on community structure, we used triplicate samples of each land use type from different locations. In addition, composite samples derived from five soil cores after the separation of soil horizons were employed.     

感染球孢白僵菌后红火蚁体内蛋白质含量的变化

为了解球孢白僵菌Beauveria bassiana成功感染红火蚁Solenopsis invicta的生理生化机制, 本研究在室内条件下测定了经球孢白僵菌感染后红火蚁工蚁、蛹和幼蚁体内蛋白质含量的变化。结果表明: 红火蚁各虫态感染球孢白僵菌后体内蛋白质含量变化存在差异。在接种后12-72 h, 工蚁体内蛋白质含量在24 h时最高, 为47.12 mg/g, 显著低于对照的56.40 mg/g(P＜0.05), 此后呈下降趋势, 72 h时蛋白质含量下降至25.16 mg/g。幼蚁在接种36 h时体内蛋白质含量最高, 为24.13 mg/g, 此后呈下降趋势, 72 h时蛋白质含量为8.95 mg/g, 显著低于对照的24.80 mg/g。蛹感染球孢白僵菌后体内蛋白质含量一直呈下降趋势, 从12 h的36.57 mg/g降到72 h的10.42 mg/g, 与对照的38.61 mg/g相比差异显著(P＜0.05)。结果显示球孢白僵菌感染能显著降低红火蚁各虫态体内蛋白质含量。