A novel bioinformatics pipeline to discover genes related to arbuscular mycorrhizal symbiosis based on their evolutionary conservation pattern among higher plants.

BACKGROUND: Genes involved in arbuscular mycorrhizal (AM) symbiosis have been identified primarily by mutant screens, followed by identification of the mutated genes (forward genetics). In addition, a number of AM-related genes has been identified by their AM-related expression patterns, and their function has subsequently been elucidated by knock-down or knock-out approaches (reverse genetics). However, genes that are members of functionally redundant gene families, or genes that have a vital function and therefore result in lethal mutant phenotypes, are difficult to identify. If such genes are constitutively expressed and therefore escape differential expression analyses, they remain elusive. The goal of this study was to systematically search for AM-related genes with a bioinformatics strategy that is insensitive to these problems. The central element of our approach is based on the fact that many AM-related genes are conserved only among AM-competent species. RESULTS: Our approach involves genome-wide comparisons at the proteome level of AM-competent host species with non-mycorrhizal species. Using a clustering method we first established orthologous/paralogous relationships and subsequently identified protein clusters that contain members only of the AM-competent species. Proteins of these clusters were then analyzed in an extended set of 16 plant species and ranked based on their relatedness among AM-competent monocot and dicot species, relative to non-mycorrhizal species. In addition, we combined the information on the protein-coding sequence with gene expression data and with promoter analysis. As a result we present a list of yet uncharacterized proteins that show a strongly AM-related pattern of sequence conservation, indicating that the respective genes may have been under selection for a function in AM. Among the top candidates are three genes that encode a small family of similar receptor-like kinases that are related to the S-locus receptor kinases involved in sporophytic self-incompatibility. CONCLUSIONS: We present a new systematic strategy of gene discovery based on conservation of the protein-coding sequence that complements classical forward and reverse genetics. This strategy can be applied to diverse other biological phenomena if species with established genome sequences fall into distinguished groups that differ in a defined functional trait of interest.

Global Genomic Diversity of Human Papillomavirus 6 Based on 724 Isolates and 190 Complete Genome Sequences.

Human papillomavirus type 6 (HPV6) is the major etiological agent of anogenital warts and laryngeal papillomas and has been included in both the quadrivalent and nonavalent prophylactic HPV vaccines. This study investigated the global genomic diversity of HPV6, using 724 isolates and 190 complete genomes from six continents, and the association of HPV6 genomic variants with geographical location, anatomical site of infection/disease, and gender. Initially, a 2,800-bp E5a-E5b-L1-LCR fragment was sequenced from 492/530 (92.8%) HPV6-positive samples collected for this study. Among them, 130 exhibited at least one single nucleotide polymorphism (SNP), indel, or amino acid change in the E5a-E5b-L1-LCR fragment and were sequenced in full. A global alignment and maximum likelihood tree of 190 complete HPV6 genomes (130 fully sequenced in this study and 60 obtained from sequence repositories) revealed two variant lineages, A and B, and five B sublineages: B1, B2, B3, B4, and B5. HPV6 (sub)lineage-specific SNPs and a 960-bp representative region for whole-genome-based phylogenetic clustering within the L2 open reading frame were identified. Multivariate logistic regression analysis revealed that lineage B predominated globally. Sublineage B3 was more common in Africa and North and South America, and lineage A was more common in Asia. Sublineages B1 and B3 were associated with anogenital infections, indicating a potential lesion-specific predilection of some HPV6 sublineages. Females had higher odds for infection with sublineage B3 than males. In conclusion, a global HPV6 phylogenetic analysis revealed the existence of two variant lineages and five sublineages, showing some degree of ethnogeographic, gender, and/or disease predilection in their distribution. IMPORTANCE: This study established the largest database of globally circulating HPV6 genomic variants and contributed a total of 130 new, complete HPV6 genome sequences to available sequence repositories. Two HPV6 variant lineages and five sublineages were identified and showed some degree of association with geographical location, anatomical site of infection/disease, and/or gender. We additionally identified several HPV6 lineage- and sublineage-specific SNPs to facilitate the identification of HPV6 variants and determined a representative region within the L2 gene that is suitable for HPV6 whole-genome-based phylogenetic analysis. This study complements and significantly expands the current knowledge of HPV6 genetic diversity and forms a comprehensive basis for future epidemiological, evolutionary, functional, pathogenicity, vaccination, and molecular assay development studies.

Evaluation of the potential use of trabecular bone score to complement bone mineral density in the diagnosis of osteoporosis: a preliminary spine BMD-matched, case-control study.

The trabecular bone score (TBS) is a new parameter that is determined from gray-level analysis of dual-energy X-ray absorptiometry (DXA) images. It relies on the mean thickness and volume fraction of trabecular bone microarchitecture. This was a preliminary case-control study to evaluate the potential diagnostic value of TBS as a complement to bone mineral density (BMD), by comparing postmenopausal women with and without fractures. The sample consisted of 45 women with osteoporotic fractures (5 hip fractures, 20 vertebral fractures, and 20 other types of fracture) and 155 women without a fracture. Stratification was performed, taking into account each type of fracture (except hip), and women with and without fractures were matched for age and spine BMD. BMD and TBS were measured at the total spine. TBS measured at the total spine revealed a significant difference between the fracture and age- and spine BMD-matched nonfracture group, when considering all types of fractures and vertebral fractures. In these cases, the diagnostic value of the combination of BMD and TBS likely will be higher compared with that of BMD alone. TBS, as evaluated from standard DXA scans directly, potentially complements BMD in the detection of osteoporotic fractures. Prospective studies are necessary to fully evaluate the potential role of TBS as a complementary risk factor for fracture.

MP2RAGE Multiple Sclerosis Magnetic Resonance Imaging at 3 T.

OBJECTIVES: Lesion detection and characterization in multiple sclerosis (MS) are an essential part of its clinical diagnosis and an important research field. In this pilot study, we applied the recently introduced two inversion-contrast magnetization-prepared rapid gradient echo sequence (MP2RAGE) to patients with early-stage MS.¦MATERIALS AND METHODS: The MP2RAGE is a 3-dimensional (3D) magnetization-prepared rapid gradient echo derivative providing homogeneous T1 weighting and simultaneous T1 mapping. The MP2RAGE performance was compared with that of 2 clinical routine sequences (2D fluid-attenuated inversion recovery [FLAIR] and 3D magnetization-prepared rapid gradient echo [MP-RAGE]) and 2 state-of-the art clinical research sequences (the 3D FLAIR-SPACE [sampling perfection with application-optimized contrasts by using different flip-angle evolutions], a fluid-attenuated variable flip-angle fast spin echo technique, and the 3D double-inversion recovery SPACE). A cohort of 10 early-stage female MS patients (age, 31.6 ± 4.7 years; disease duration, 3.8 ± 1.9 years; median expanded disability status scale score, 1.75) and 10 age- and gender-matched controls were enrolled after approval of the local institutional review board was obtained. Multiple sclerosis lesions were identified and assigned to brain locations and tissue types by two experienced physicians in all 5 contrasts. Subsequently, lesions were manually delineated for comparison and statistical analysis of lesion count, volume and quantitative measures.¦RESULTS AND CONCLUSIONS: The results show that the 3D T1-weighted high-resolution MP2RAGE contrast provides a sensitive means for MS lesion assessment. The additional quantitative T1 relaxation time maps obtained with the MP2RAGE provide further potential diagnostic and prognostic information that could help (a) to better discriminate lesion subtypes and (b) to stage and predict the activity and the evolution of MS. Results also indicate that the T2-weighted double-inversion recovery and FLAIR-SPACE contrasts are attractive complements to the MP2RAGE for lesion detection.

Density-based hierarchical clustering of pyro-sequences on a large scale--the case of fungal ITS1.

MOTIVATION: Analysis of millions of pyro-sequences is currently playing a crucial role in the advance of environmental microbiology. Taxonomy-independent, i.e. unsupervised, clustering of these sequences is essential for the definition of Operational Taxonomic Units. For this application, reproducibility and robustness should be the most sought after qualities, but have thus far largely been overlooked. RESULTS: More than 1 million hyper-variable internal transcribed spacer 1 (ITS1) sequences of fungal origin have been analyzed. The ITS1 sequences were first properly extracted from 454 reads using generalized profiles. Then, otupipe, cd-hit-454, ESPRIT-Tree and DBC454, a new algorithm presented here, were used to analyze the sequences. A numerical assay was developed to measure the reproducibility and robustness of these algorithms. DBC454 was the most robust, closely followed by ESPRIT-Tree. DBC454 features density-based hierarchical clustering, which complements the other methods by providing insights into the structure of the data. AVAILABILITY: An executable is freely available for non-commercial users at ftp://ftp.vital-it.ch/tools/dbc454. It is designed to run under MPI on a cluster of 64-bit Linux machines running Red Hat 4.x, or on a multi-core OSX system. CONTACT: dbc454@vital-it.ch or nicolas.guex@isb-sib.ch.

Contemporary approach to neurologic prognostication of coma after cardiac arrest.

Coma after cardiac arrest (CA) is an important cause of admission to the ICU. Prognosis of post-CA coma has significantly improved over the past decade, particularly because of aggressive postresuscitation care and the use of therapeutic targeted temperature management (TTM). TTM and sedatives used to maintain controlled cooling might delay neurologic reflexes and reduce the accuracy of clinical examination. In the early ICU phase, patients' good recovery may often be indistinguishable (based on neurologic examination alone) from patients who eventually will have a poor prognosis. Prognostication of post-CA coma, therefore, has evolved toward a multimodal approach that combines neurologic examination with EEG and evoked potentials. Blood biomarkers (eg, neuron-specific enolase [NSE] and soluble 100-β protein) are useful complements for coma prognostication; however, results vary among commercial laboratory assays, and applying one single cutoff level (eg, > 33 μg/L for NSE) for poor prognostication is not recommended. Neuroimaging, mainly diffusion MRI, is emerging as a promising tool for prognostication, but its precise role needs further study before it can be widely used. This multimodal approach might reduce false-positive rates of poor prognosis, thereby providing optimal prognostication of comatose CA survivors. The aim of this review is to summarize studies and the principal tools presently available for outcome prediction and to describe a practical approach to the multimodal prognostication of coma after CA, with a particular focus on neuromonitoring tools. We also propose an algorithm for the optimal use of such multimodal tools during the early ICU phase of post-CA coma.

Context-dependent dual role of SKI8 homologs in mRNA synthesis and turnover.

Eukaryotic mRNA transcription and turnover is controlled by an enzymatic machinery that includes RNA polymerase II and the 3' to 5' exosome. The activity of these protein complexes is modulated by additional factors, such as the nuclear RNA polymerase II-associated factor 1 (Paf1c) and the cytoplasmic Superkiller (SKI) complex, respectively. Their components are conserved across uni- as well as multi-cellular organisms, including yeast, Arabidopsis, and humans. Among them, SKI8 displays multiple facets on top of its cytoplasmic role in the SKI complex. For instance, nuclear yeast ScSKI8 has an additional function in meiotic recombination, whereas nuclear human hSKI8 (unlike ScSKI8) associates with Paf1c. The Arabidopsis SKI8 homolog VERNALIZATION INDEPENDENT 3 (VIP3) has been found in Paf1c as well; however, whether it also has a role in the SKI complex remains obscure so far. We found that transgenic VIP3-GFP, which complements a novel vip3 mutant allele, localizes to both nucleus and cytoplasm. Consistently, biochemical analyses suggest that VIP3-GFP associates with the SKI complex. A role of VIP3 in the turnover of nuclear encoded mRNAs is supported by random-primed RNA sequencing of wild-type and vip3 seedlings, which indicates mRNA stabilization in vip3. Another SKI subunit homolog mutant, ski2, displays a dwarf phenotype similar to vip3. However, unlike vip3, it displays neither early flowering nor flower development phenotypes, suggesting that the latter reflect VIP3's role in Paf1c. Surprisingly then, transgenic ScSKI8 rescued all aspects of the vip3 phenotype, suggesting that the dual role of SKI8 depends on species-specific cellular context.

Improved predictive mapping of indoor radon concentrations using ensemble regression trees based on automatic clustering of geological units.

PURPOSE: According to estimations around 230 people die as a result of radon exposure in Switzerland. This public health concern makes reliable indoor radon prediction and mapping methods necessary in order to improve risk communication to the public. The aim of this study was to develop an automated method to classify lithological units according to their radon characteristics and to develop mapping and predictive tools in order to improve local radon prediction. METHOD: About 240 000 indoor radon concentration (IRC) measurements in about 150 000 buildings were available for our analysis. The automated classification of lithological units was based on k-medoids clustering via pair-wise Kolmogorov distances between IRC distributions of lithological units. For IRC mapping and prediction we used random forests and Bayesian additive regression trees (BART). RESULTS: The automated classification groups lithological units well in terms of their IRC characteristics. Especially the IRC differences in metamorphic rocks like gneiss are well revealed by this method. The maps produced by random forests soundly represent the regional difference of IRCs in Switzerland and improve the spatial detail compared to existing approaches. We could explain 33% of the variations in IRC data with random forests. Additionally, the influence of a variable evaluated by random forests shows that building characteristics are less important predictors for IRCs than spatial/geological influences. BART could explain 29% of IRC variability and produced maps that indicate the prediction uncertainty. CONCLUSION: Ensemble regression trees are a powerful tool to model and understand the multidimensional influences on IRCs. Automatic clustering of lithological units complements this method by facilitating the interpretation of radon properties of rock types. This study provides an important element for radon risk communication. Future approaches should consider taking into account further variables like soil gas radon measurements as well as more detailed geological information.

Experimental Cournot oligopoly and inequity aversion

This paper explores the role of inequity aversion as an explanation for observed behavior in experimental Cournot oligopolies. We show that inequity aversion can change the nature of the strategic interaction: quantities are strategic substitutes for sufficiently asymmetric output levels but strategic complements otherwise. We find that inequity aversion can explain why: (i) some experiments result in higher than Cournot-Nash production levels while others result in lower, (ii) collusion often occurs with only two players whereas with three or more players market outcomes are very close to Cournot-Nash, and (iii) players often achieve equal profits in asymmetric Cournot oligopoly.